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SUMMARY 


This  report  presents  the  mathematical  formulation  and  verification  of  a 
numerical  model  that  simulates  wave  transformation  and  longshore  current  over 
a  bar  and  trough  profile.  The  model  is  intended  for  application  on  projects 
involving  dredged  material  berms  and  other  situations  where  a  bar  and  trough 
topography  prevents  use  of  simple  analytic  or  numerical  solutions  that  are 
restricted  to  a  uniformly  sloping  beach.  The  model  is  called  NMLONG,  an 
acronym  that  stands  for  Numerical  Model  of  the  LONGshore  current,  and  it  can 
be  conveniently  run  on  a  desk- top  (personal)  computer. 

The  model  incorporates  all  known  features  of  the  wave  and  longshore 
current  system  that  appear  in  research- type  engineering  models  run  on  main¬ 
frame  computer  systems.  These  features  include  wave  and  wind  driving,  wave 
breaking  and  reformation  over  multiple  bar  and  trough  profiles,  and  lateral 
mixing.  The  model  also  allows  choice  of  linear  or  quadratic  bottom  friction 
and  regular  or  random  wave  heights  as  options.  The  main  restriction  of  the 
model  is  longshore  uniformity  of  the  waves  and  beach  topography. 

The  model  was  verified  with  several  field  and  laboratory  data  sets.  The 
wave  calculation  reproduced  cross -shore  heights  of  laboratory  monochromatic 
waves  on  both  plane-sloping  and  bar  and  trough  profiles,  and  it  also  repro¬ 
duced  field  measurements  of  wave  height  frequency  of  occurrence  at  different 
depths  in  the  surf  zone .  The  longshore  current  model  reproduced  measured 
currents  on  plane -sloping  laboratory  beaches  and  currents  measured  in  the 
field  on  irregular  bottoms.  Sensitivity  analyses  examined  linear  and  nonlin¬ 
ear  bottom  friction  formulations,  strength  of  lateral  mixing,  effect  of  wind 
on  setup/setdown  and  the  current,  and  wave-current  interaction. 

Several  time-saving  algorithms  were  developed  that  enable  NMLONG  to  run 
effectively  on  a  personal  computer.  The  model  also  contains  an  interface  that 
allows  data  entry  and  graphical  display  of  results  similar  to  a  commercial 
software  product.  Range  checking  of  input  data  and  default  values  of  princi¬ 
pal  model  parameters  is  also  provided.  The  graphical  output  consists  of  the 
beach  profile  and  the  wave  height  and  longshore  current  velocity  along  the 
profile.  Output  data  files  include  the  wave  height,  wave  direction,  longshore 
current  velocity,  and  water  surface  elevation  (wave-  and  wind- induced  setup 
and  setdown)  across  the  beach  profile . 
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NMLONG :  NUMERICAL  MODEL  FOR  SIMULATING  THE  LONGSHORE  CURRENT 
MODEL  DEVELOPMENT  AND  TESTS 

PART  I :  INTRODUCTION 

Problem  Statement  and  Objectives 

1.  The  longshore  current  is  a  major  and  pervasive  flow  in  the  near¬ 
shore,  and  estimates  of  its  velocity  are  needed  in  most  coastal  engineering 
projects  involving  sediment  transport,  such  as  navigation  channel  maintenance 
and  shore  protection.  The  distribution  of  the  longshore  current  across  the 
nearshore  zone  depends  on  many  factors  that  are  described  in  detail  in  this 
report.  This  distribution  is  complex  and  typically  has  one  or  more  peaks,  for 
example,  one  peak  just  shoreward  of  the  wave  breaker  line  and  another  on  the 
foreshore . 

2.  One  of  the  major  purposes  for  calculating  the  longshore  current  is 
to  estimate  the  cross -shore  distribution  of  the  longshore  sand  transport  rate. 
For  example,  the  amount  of  sand  transported  around  the  tip  of  a  jetty  and, 
conversely,  the  amount  impounded  by  the  jetty  correspond  with  the  profile  of 
the  longshore  current.  If  average  and  seasonal  trends  in  the  current  and  sand 
transport  patterns  at  the  jetty  can  be  established,  channel  dredging  and  sand 
bypassing  schedules  and  resources  can  be  more  accurately  planned.  In  general, 
beach  evolution  is  controlled  by  spatial  and  temporal  changes  in  the  longshore 
sand  transport  rate,  which  are  directly  related  to  changes  in  the  longshore 
current  velocity. 

3 .  On  many  beaches ,  longshore  bars  and  troughs  appear  at  intervals 
offshore  and  follow  the  coast  at  almost  fixed  distances  from  the  shore.  If  a 
coast  has  a  rich  sand  supply,  a  large  bar  will  tend  to  form  at  the  mouths  of 
inlets  and  harbor  entrances.  These  bars  are  created  in  part  by  the  sand 
transported  along  the  coast  by  the  longshore  current  and  have  been  hypothe¬ 
sized  to  act  as  conduits  of  sand  around  inlets  and  entrances.  Longshore  sand 
transport  is  thus  a  major  factor  in  navigation  channel  infilling  and  bypassing 
of  sand  around  inlets  and  entrances.  An  estimate  of  the  transport  rate 
requires  knowledge  of  the  cross-shore  distribution  of  the  longshore  current. 
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4.  In  recent  years,  the  concept  of  placing  dredged  material  in  the 
nearshore  in  the  form  of  shore-parallel  longshore  bars  or  "underwater  berms" 
has  been  advanced  as  a  beneficial  use  of  dredged  material  for  shore  protection 
and  also  as  a  means  to  reduce  dredging  costs  (Murden  1988,  McLellan  1990). 

Such  berms  may  act  to  break  erosive  storm  waves  and,  possibly,  serve  as  a 
source  of  beach  nourishment.  It  is  necessary  to  know  the  properties  of  the 
longshore  current  profile  over  these  bar  and  trough  features  both  to  under¬ 
stand  the  nearshore  hydrodynamic  environment  and  to  estimate  the  longshore 
sand  transport  rate  and  evolution  of  these  bars. 

5.  The  above  overview  of  common  coastal  engineering  situations  involv¬ 
ing  longshore  sand  transport  leads  to  an  appreciation  that  realistic  calcula¬ 
tion  of  the  longshore  current  can  be  made  only  under  the  general  conditions  of 
arbitrary  bottom  topography  (necessarily  restricted  to  be  uniform  alongshore 
in  the  present  study,  as  is  discussed  in  the  following  two  chapters)  and  arbi¬ 
trary  waves,  wind,  and  water  level.  Both  waves  and  wind  are  driving  forces  of 
the  longshore  current,  and  both  can  change  the  water  level  near  the  shore.  On 
a  barred  bottom  topography,  waves  may  break,  reform,  and  break  again,  implying 
the  existence  of  a  multipeaked  distribution  of  the  longshore  current  across 
the  surf  zone  and  a  corresponding  multipeaked  distribution  of  the  longshore 
sand  transport  rate.  The  independent  actions  of  wind  and  waves  can  either 
positively  or  negatively  reinforce  in  generating  the  current,  again  creating  a 
complex  current  profile. 

6.  The  main  objective  of  the  present  study  is  development  of  a  numer¬ 
ical  model  of  the  longshore  current  applicable  to  nearshore  dredged  material 
placement  operations  involving  creation  of  berms  and  coastal  engineering 
projects  involving  longshore  currents  and  longshore  sand  transport.  An 
additional  objective  is  to  bring  the  developed  technology  to  the  desk  of  the 
engineer  or  planner  in  the  form  of  a  convenient  computer  program  that  can  be 
easily  run  on  a  desk- top  or  personal  computer.  By  use  of  modern  mathematical 
and  numerical  procedures,  the  resultant  program  has  the  desired  generality  and 
reliability;  yet  because  of  its  structured  interface  and  graphical  output,  it 
is  more  convenient  than  following  a  handbook  procedure . 
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Model  Design  Requirements 


7.  The  model  was  developed  to  satisfy  requirements  dictated  by  known 
and  anticipated  applications.  The  following  is  a  list  of  those  requirements: 

a.  The  model  should  calculate  the  longshore  current  profile  over  an 
arbitrary  sea  bottom  of  plane  and  parallel  contours.  Such  a 
profile  could  consist  of  multiple  longshore  bars  and  troughs. 

b.  The  model  solution  scheme  should  represent  the  physical  phenom¬ 
ena  in  the  most  accurate  way  consistent  with  engineering  data, 
uses  of  the  model ,  and  accuracy  of  the  governing  equations . 

c.  The  model  should  require  only  readily  available  or  estimated 
input  data. 

d.  Both  wave  and  wind  forcing  should  be  represented. 

e.  The  change  in  mean  water  level,  which  enters  in  the  longshore 
current  calculation,  should  be  calculated. 

f.  The  model  should  be  robust,  that  is,  be  numerically  stable  and 
provide  reliable  answers  for  a  wide  range  of  conditions  without 
requiring  significant  changes  in  empirical  factors  entering  the 
model . 

g.  The  model  should  be  sufficiently  flexible  to  allow  ready  modifi¬ 
cation  and  incorporation  of  new  developments. 

h.  The  model  should  be  verified  with  data  encompassing  a  wide  range 
of  wave,  wind,  current,  and  topographic  conditions. 

jL.  The  model  should  be  sufficiently  efficient  to  run  on  a  desk- top 
computer . 

8.  These  requirements  necessarily  led  to  development  of  a  numerical 
modeling  system,  called  NMLONG  for  Numerical  Model  of  the  LONGshore  current, 
which  consists  of  a  wave  transformation  model  and  a  longshore  current  model. 
The  basic  structure  of  a  wave  transformation  model  meeting  the  requirements 
was  available  from  the  literature  and  required  minor  refinements  and  exten¬ 
sions,  whereas  the  longshore  current  model  was  developed  directly  from  the 
governing  equations,  guided  by  previous  work. 

9.  In  the  course  of  model  development,  a  number  of  fundamental  topics 
were  investigated.  These  include  the  significance  of  the  wave  and  current 
interaction,  comparison  of  the  effect  of  linear  and  quadratic  bottom  friction 
laws,  and  relative  strength  of  the  wind  and  waves  in  generating  a  longshore 
current. 
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Report  Contents 


10.  Part  I  gives  the  motivation  and  objectives  of  the  study.  An 
orientation  to  the  topic  of  longshore  current  calculations  and  review  of 
pertinent  studies  are  presented  in  Part  II.  Parts  III  and  IV  describe  wave 
model  development  and  verification,  respectively,  and  Parts  V  and  VI  similarly 
pertain  to  the  longshore  current  model.  A  summary  of  main  results  is  given  in 
Part  VII.  Discussion  of  the  numerical  method  used  to  efficiently  perform  the 
time  average  of  the  nonlinear  bottom  friction  stress  is  given  in  Appendix  A, 
and  data  sets  used  in  model  verification  are  listed  in  Appendix  B. 

11.  The  present  report  describes  model  development  and  tests.  Report  2 
in  the  NMLONG  series  will  be  a  user's  guide  and  description  of  the  implementa¬ 
tion  of  the  model  on  a  desk- top  computer. 
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LITERATURE  REVIEW 

verview 

ascription  of  the  longshore  current,  which 
'.rculation,  was  independently  and  simul- 
Longuet-Higgins  (1970a,  b)  and  Thornton 
flux  or  radiation  stress  concepts  in- 
rt  (1962,  1963,  1964).  Because  of  the 
in  coastal  engineering  and  science,  over 
.  mathematical  studies  and,  to  a  lesser 
tave  been  conducted  to  improve  under - 
e  considerable  literature  on  the 
st  classical"  in  a  certain  sense;  yet 

iiade  of  the  literature  on  longshore 
o  and  Coleman  1982,  Gourlay  1982,  Komar 
is  so  important  that  new  analytic 

.  - or  the  governing  equations  are  published  almost 

yearly.  Considerable  mathematical  ingenuity  and  clever  specification  of 
idealized  conditions  have  gone  into  development  of  these  solutions.  Of  all 
the  analytic  solutions,  the  model  of  Longuet-Higgins  (1970a,  b)  is  remarkably 
insightful  yet  simple.  This  model  has  served  as  the  starting  point  for  a 
large  number  of  solutions  that  have  removed  one  or  more  restrictions  from  the 
original  model  or  added  new  capabilities.  As  analytic  solutions  are  made  more 
general,  however,  their  evaluation  becomes  mathematically  and  computationally 
complex.  The  authors  believe  that  numerical  models  such  as  described  in  this 
report  are  more  computationally  efficient  than  the  complex  analytic  or  "exact- 
solutions  and  are  easier  to  modify  to  incorporate  improvements.  In  contrast, 
a  particular  analytic  solution  has  very  little  flexibility  to  admit  change 
owing  to  the  conditions  imposed  to  obtain  that  solution. 

14.  The  alternative  approach  to  analytic  models  of  the  longshore 
current  is  to  solve  the  governing  equations  numerically.  Indeed,  in  one  of 
the  original  formulations,  Thornton  (1970)  developed  a  numerical  model  and 
solved  the  governing  differential  equations  under  much  less  restrictive 
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conditions  than  are  possible  with  analytic  solutions.  Since  that  time, 
numerous  numerical  models  of  the  longshore  current  have  been  published. 

15.  This  literature  review  is  divided  by  calculation  approaches,  dis¬ 
cussing  analytic  and  numerical  models  separately.  It  begins  by  introducing 
important  terminology,  concepts,  and  field  observations.  The  review  is  not 
comprehensive  in  that  articles  are  discussed  to  compile  results  important  for 
understanding  the  background  of  the  problem  and  assessing  the  capabilities  of 
previous  models  and  the  present  model. 

Terminology 

16.  In  this  section  selected  terms  that  appear  in  the  literature  on  the 
longshore  current  are  introduced. 

Velocity 

17.  The  word  "velocity"  in  a  strict  sense  denotes  a  vector  quantity 
that  has  a  magnitude  (speed)  and  a  direction.  In  the  case  of  the  longshore 
current,  the  flow  is  along  the  coordinate  axis  oriented  parallel  to  the  coast, 
and  the  flow  direction  can  be  only  positive  or  negative  along  this  axis. 
Therefore,  the  phrase  "longshore  current  velocity"  can  be  used  to  denote  speed 
with  the  understanding  that  the  sign  gives  the  direction,  since  this  current 
is  only  one  component  of  the  velocity  of  the  general  two-dimensional  (2-D) 
horizontal  circulation  pattern.  By  the  same  convention,  vector  symbols  are 
not  used  in  writing  equations  containing  the  longshore  current  velocity. 
Profile  and  distribution 

18.  Another  phrase  that  commonly  appears  in  the  literature  is  "long¬ 
shore  current  profile,"  which  refers  to  the  distribution  of  the  depth- averaged 
current  from  the  mean  water  shoreline  to  a  point  seaward  of  the  breaking  wave 
zone  where  the  wave -induced  longshore  current  becomes  negligible  (typically 
one  or  two  times  the  width  of  the  surf  zone).  In  contrast,  referring  to  the 
distribution  of  the  wave  height  across  the  surf  zone  as  the  "wave  profile"  is 
not  appropriate,  because  that  term  normally  refers  to  the  form  of  an  indivi¬ 
dual  wave.  However,  the  word  "distribution"  in  reference  to  wave  height  can 
cause  confusion  because  it  may  also  describe  the  statistical  distribution  of 
wave  height  at  a  single  point,  for  example,  the  probability  density  function 
of  wave  height.  In  this  report,  both  uses  of  the  word  "distribution"  are 
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necessary  in  discussion  of  the  wave  height,  but  in  referring  to  its  spatial 
distribution,  the  term  "wave  height  distribution  across  the  shore"  or  similar 
phrase  will  be  employed. 

Steady  state  or  mean 

19 .  The  model  developed  in  this  report  and  in  most  studies  on  the 
longshore  current  pertains  to  the  steady- state  solution  of  the  governing 
equations,  i.e.,  time-averaged  or  mean  quantities  as  opposed  to  instantaneous 
quantities.  Therefore,  unless  otherwise  stated,  reference  to  the  longshore 
current  pertains  to  the  steady  or  mean  longshore  current.  The  temporal  or 
transient  behavior  of  the  longshore  current  has  been  little  studied,  although, 
for  example,  accelerations  and  decelerations  in  the  water  are  expected  to  be 
important  in  mobilizing  and  moving  sand.  All  modeling  studies  that  calculate 
the  temporal  behavior  of  the  longshore  current  have  merely  used  the  time- 
dependent  term  in  the  governing  equation  to  step  in  time  to  a  steady- state 
solution  without  detailed  examination  of  the  transient  behavior  of  the 
current.  In  engineering  approaches  to  calculating  the  longshore  current, 
temporal  changes  in  the  current  have  been  simulated  by  moving  from  one  steady- 
state  solution  to  another  in  discrete  steps.  This  approach  is  used  here,  but 
in  the  future  the  temporal  behavior  of  the  current  will  need  to  be  dealt  with 
more  rigorously. 


Field  Observations 

20.  In  this  section  major  concepts  and  processes  pertaining  to  long¬ 
shore  currents  observed  in  the  field  are  discussed  to  understand  the  various 
phenomena  that  a  mathematical  model  must  describe  and  potential  limitations 
that  may  be  manifested  in  such  an  approach. 

Nearshore  current  system 

21.  The  term  "nearshore  currents"  was  introduced  by  Shepard  and  Inman 
(1950)  to  distinguish  currents  generated  directly  by  incident  waves  from 

coastal  currents  associated  with  the  tide  and  large-scale  oceanic  circula¬ 
tion  systems .  In  this  report  the  local  wind  is  also  taken  to  be  a  forcing 
mechanism  of  the  nearshore  current.  The  nearshore  current  system  consists  of 
the  longshore  current  that  flows  parallel  to  the  trend  of  the  shore,  currents 
that  flow  strongly  offshore  in  narrow  bands  called  "rip  currents"  (Shepard, 
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Emery,  and  La  Fond  1941),  and  weaker  cross-shore  currents  that  flow  onshore 
and  offshore.  The  boundaries  of  the  nearshore  circulation  system  are  dynamic: 
the  shoreward  limit  is  the  highest  runup  in  a  time -dependent  description  or 
the  mean  shoreline  in  the  case  of  a  time -average  description,  and  the  seaward 
limit  is  in  the  vicinity  of  the  seaward-most  tip  of  rip  current  heads,  which 
extend  the  length  of  about  two  or  three  widths  of  the  surf  zone.  The  stron¬ 
gest  nearshore  currents  are  located  in  the  surf  zone,  which  is  the  region  of 
intensive  depth-limited  wave  breaking  and  maximum  longshore  current  speed,  and 
in  the  throats  of  rip  currents  with  their  strong  seaward  flow.  In  the  present 
situation,  lateral  boundaries  are  assumed  to  be  located  sufficiently  distant 
so  as  not  to  influence  the  wave  and  current  field  in  the  area  of  interest. 

22.  On  the  basis  of  extensive  field  observations,  Harris  (1969)  clas¬ 
sified  nearshore  current  patterns  primarily  by  the  angle  of  wave  incidence  as: 

a.  Symmetrical  and  cellular,  in  which  rip  currents  divide  the 
nearshore  into  clear  cells  delineated  by  rip  currents,  tending 
to  form  under  near-normal  wave  incidence. 

b.  Alongshore  system,  in  which  the  current  mainly  flows  parallel 
to  the  shoreline,  and  tending  to  form  under  large  oblique 
angles  of  wave  incidence. 

c.  Asymmetrical  and  cellular,  an  intermediate  pattern  between  the 
symmetrical  and  cellular  pattern  and  the  alongshore  system, 
tending  to  form  under  moderate  angles  of  wave  incidence. 

23.  From  his  synoptic  field  measurements  of  waves,  currents,  and  bottom 
bathymetry,  Sonu  (1972)  concluded  that  the  shape  of  the  bottom  was  an  impor¬ 
tant  factor  in  addition  to  the  angle  of  wave  incidence  in  controlling  the 
pattern  of  the  nearshore  circulation.  In  a  classification  scheme  similar  to 
that  of  Harris  (1969) ,  Sonu  defined  a  cellular  pattern,  meandering  pattern, 
and  longshore  pattern.  In  a  less  strict  sense,  the  longshore  component  of  an 
asymmetrical  and  cellular  (or  meandering  current)  system  may  be  approximated 
as  an  alongshore  system  if  an  average  is  taken  over  sufficiently  long  time. 
Longshore  uniformity 

24.  Conditions  most  favorable  for  the  existence  of  a  longshore  current 
are  large  angles  of  wave  incidence  and  relatively  uniform  bottom  contours 
alongshore.  If  the  contours  are  not  uniform,  the  nearshore  current  will 
necessarily  have  components  directed  onshore  or  offshore.  Similarly,  if  the 
wave  field  is  not  uniform  alongshore,  the  longshore  current  may  reverse  or 
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weaken  at  points,  and  the  nearshore  current  pattern  may  tend  to  form  rip 
currents  with  longshore  currents  serving  as  "feeder  currents"  between  the  rips 
(Bowen  1969).  Although  it  is  not  necessary  to  invoke  longshore  uniformity  in 
development  of  a  longshore  current  model  (e.g.,  Komar  1975,  Gourlay  1976, 

Keely  and  Bowen  1977),  it  is  most  realistic  and  practical  to  do  so.  The  most 
important  exception  in  common  engineering  applications  where  longshore 
uniformity  is  not  a  good  approximation  involves  diffraction  at  headlands, 
detached  breakwaters,  and  long  jetties,  which  creates  longshore  gradients  in 
wave  height  and  wave  angle.  These  gradients  can  produce  currents  of  magnitude 
comparable  to  that  of  oblique  wave  incidence  (Gourlay  1976). 

"Unsteady"  steady  current 

25.  Although  the  concept  of  a  steady  longshore  current  is  a  convenient 
construct,  it  is  only  an  approximate  description  of  the  longshore  current 
existing  in  the  field.  By  using  tethered  floats  configured  to  reside  either 
near  the  water  surface  or  near  the  bottom,  Inman  and  Quinn  (1952)  measured 
longshore  currents  on  a  Pacific  Ocean  beach  and  found  that  the  standard 
deviation  in  the  current  could  exceed  its  mean  value.  In  nature,  the  forcing 
conditions,  principally  the  waves  but  also  the  wind,  are  not  stationary,  and 
the  longshore  current  shows  significant  variations  in  time  (Inman  and  Quinn 
1952,  Wood  and  Meadows  1975,  Meadows  1976)  and  along  the  shore  (Keely  and 
Bowen  1977)  .  Guza  and  Thornton  (1978)  discuss  temporal  and  spatial  varia¬ 
bility  in  their  longshore  current  field  measurements,  review  the  literature  on 
the  phenomenon,  and  suggest  that  large-scale  eddy  motions  and  edge  waves  may 
be  responsible. 

Long-period  waves 

26.  Surface  gravity  waves  that  are  visibly  observed  to  arrive  at  the 
coast  may  be  classified  as  short-period  waves  with  a  typical  period  on  the 
order  of  10  sec.  It  is  these  waves  that  are  taken  to  be  one  of  the  two 
driving  mechanisms  of  the  longshore  current  in  this  report,  the  other  driving 
mechanism  being  the  wind.  Other  types  of  wave  motion  in  the  nearshore  not 
readily  visible  to  the  eye  have  periods  longer  than  the  incident  gravity 
waves.  The  so-called  infragravity  waves,  wave  with  periods  longer  (for 
example,  100  sec)  than  gravity  waves,  are  associated  with  wave  groups  or  surf 
beat,  edge  waves,  and  standing  waves.  These  waves  may  modulate  the  longshore 
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current  (Shepard  and  Inman  1950,  Inman  and  Quinn  1952)  but  are  not  expected  to 
greatly  change  the  magnitude  of  the  mean  longshore  current. 

27.  Recently,  a  new  type  of  long  wave,  called  a  shear  wave  or  far 
infragravity  wave  with  a  period  in  the  range  of  100  to  1,000  sec  has  been  dis¬ 
covered  (Bowen  and  Holman  1990;  Oltman-Shay,  Howd,  and  Birkemeier  1990)  that 
may  at  times  contribute  to  the  total  longshore  current  speed  which  is  on  the 
same  order  as  the  current  produced  by  oblique  (gravity)  wave  incidence.  Shear 
waves  originate  from  instabilities  in  the  shear  in  the  profile  of  the  steady 
(mean)  wave- induced  longshore  current  discussed  in  this  report.  At  present, 
knowledge  of  the  formation  and  persistence  of  shear  currents  is  only  in  its 
infancy.  Since  these  waves  owe  their  existence  to  the  presence  of  a  longshore 
current,  accurate  prediction  of  the  profile  of  a  mean  current  takes  on  even 
greater  importance . 

Paradox  of  lateral  mixing 

28.  As  in  other  types  of  hydrodynamic  and  circulation  studies,  the  term 
"lateral  mixing"  refers  to  the  exchange  of  fluid  momentum  in  the  horizontal 
plane,  in  the  present  context  perpendicular  to  the  shore.  Lateral  exchange  of 
momentum  mixes  water  particles  having  different  longshore  velocities,  thereby 
smearing  the  current  profile.  Significant  smoothing  of  the  longshore  current 
profile  in  the  surf  zone  is  expected  because  of  the  intense  turbulence  created 
by  breaking  waves,  which  is  advected  back  and  forth  by  oscillatory  wave 
motion.  In  mathematical  modeling  of  longshore  current  profiles  generated  in 
carefully  performed  laboratory  experiments  with  incident  waves  of  constant 
height,  period,  and  direction,  it  is  necessary  to  include  lateral  mixing 
inside  the  surf  zone  as  well  as  outside  where  the  tail  of  the  profile  extends 
substantially  seaward.  On  the  other  hand,  some  investigators  (e.g.,  Battjes 
1972,  Thornton  and  Guza  1986)  have  suggested  that  in  the  field  the  combination 
of  waves  of  different  heights  and  periods  from  different  directions  will 
produce  a  smooth  current  profile,  greatly  reducing  the  lateral  mixing  effect, 
and  perhaps  eliminating  the  necessity  to  include  it  in  models.  This  obser¬ 
vation  leads  to  the  "lateral  mixing  paradox,"  the  fact  that  in  most  modeling 
studies  approximately  the  same  strength  of  lateral  mixing  has  been  required  to 
describe  longshore  currents  in  the  laboratory  generated  by  regular  waves  as  in 
the  field  with  currents  calculated  by  regular  or  as,  in  the  present  study, 
random  waves . 
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29.  A  possible  resolution  of  this  paradox  may  lie  in  the  recent 
discovery  of  the  shear  waves  described  in  the  previous  section  (Bowen  and 
Holman  1990;  Oltman-Shay,  Howd,  and  Birkemeier  1990).  Bowen  and  Holman  note 
that  the  greater  cross -shore  length  scale  of  a  shear  wave  (order  of  50  m)  as 
compared  with  the  depth- limited  orbital  wave  motion  of  a  gravity  wave  in  the 
surf  zone  (order  of  5  m)  suggests  that  shear  waves  may  be  a  dominant  mixing 
mechanism,  spreading  the  wave -induced  turbulence  across  the  nearshore.  Bowen 
and  Holman  also  state  that  "wave  tank  experiments  (of  the  longshore  current) 
may  in  fact  be  contaminated  with  shear  wave  motions"  that  cannot  be  detected 
by  normal  means.  The  presence  of  shear  waves  is  a  plausible  reason  for  the 
high  tails  in  the  longshore  current  profile  seaward  of  the  breaker  line  found 
in  some  laboratory  experiments.  Other  reasons  for  such  high  tails  are 
circulation  in  a  closed  basin  and  direct  wave -induced  lateral  mixing  combined 
with  the  oscillating  position  of  the  breakers  that  occurs  even  in  experiments 
with  regular  waves . 

Effect  of  wind 

30.  The  wind  contributes  in  at  least  three  ways  as  a  forcing  function 
in  nearshore  hydrodynamics : 

a.  It  creates  longshore  and  cross -shore  currents. 

b.  It  changes  the  mean  water  level. 

c.  It  changes  the  characteristics  of  the  breaking  waves,  which 
feed  back  to  the  current  and  water  surface  changes . 

These  three  processes  are  discussed  in  the  following  paragraphs. 

31.  Wind  and  the  longshore  current.  Winds  blowing  over  the  sea  far 
from  a  project  site  create  a  portion  of  the  waves  that  arrive  at  the  site,  and 
these  waves  may  generate  a  longshore  current.  Besides  generating  wavesthat 
break  and  drive  currents,  the  local  wind  generates  currents  directly  by 
exerting  a  shear  on  the  water  surface.  Nummedal  and  Finley  (1978)  performed 
multiple  regression  analysis  using  both  linear  and  nonlinear  combinations  of 
observed  wave  and  wind  parameters  in  attempted  correlations  with  the  longshore 
current  velocity.  The  data  set  encompassed  1  year  of  visually  obtained  wind, 
wave,  and  current  measurements  made  at  Debidue  Island,  South  Carolina,  facing 
the  Atlantic  Ocean.  They  found  that  the  longshore  component  of  the  wind  speed 
(laS6ed  by  2  to  3  hr  between  reversals)  was  the  dominant  parameter  describing 
the  longshore  current  velocity  in  the  surf  zone.  Because  the  local  wind  also 
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generates  waves,  the  longshore  component  of  the  current  contains  a  contribu¬ 
tion  from  obliquely  incident  waves,  which  would  tend  to  raise  the  correlation. 
The  demonstration  by  Nummedal  and  Finley  of  the  significance  of  wind  on  the 
longshore  current  does  not  appear  to  have  been  adequately  appreciated  in  the 
coastal  engineering  and  science  communities. 

32.  Everts*  described  an  experience  of  SCUBA  diving  in  approximately 
13-m  depth  off  the  Coastal  Engineering  Research  Center's  (CERC's)  Field 
Research  Facility  (FRF)  at  Duck,  North  Carolina,  while  a  fast-moving  North¬ 
easter  approached  with  a  strong  longshore  wind  component.  In  a  matter  of 
minutes  after  the  wind  began  (order  of  10  min) ,  the  current  went  from  nondete- 
ctable  to  such  a  great  speed  that  the  divers  at  the  bottom  could  not  maintain 
position  against  it.  This  observation  indicates  that  the  lag  between  long¬ 
shore  current  generation  by  wind  and  changes  in  the  wind  velocity  is  very 
short  for  strong  winds  blowing  over  the  shallow  water  of  the  nearshore. 

33.  Hubertz  (1986)  reports  measurements  of  wind,  waves,  and  longshore 
current  made  at  the  FRF.  A  direct  relationship  between  the  mean  longshore 
current  and  wind  speed  was  found,  and  Hubertz  recommended  that  meteorological 
factors  be  included  in  longshore  current  estimation  techniques.  Hubertz 
(1984)  also  performed  numerical  simulations  of  the  wind-  and  wave-driven 
longshore  current  and  found  that  inclusion  of  wind  forcing  significantly 
improved  agreement  with  the  measurements  over  simulations  made  using  only  wave 
forcing. 

34.  Wind  and  water  surface  elevation.  The  contribution  to  the  cross¬ 
shore  momentum  equation  governing  the  change  in  mean  water  level  produced  by 
an  imposed  wind  is  well  known  (e.g.,  Bretschneider  1966;  DeVries  and  Amorocho 
1980) .  It  is  included  in  the  present  model  development  together  with  wave- 
induced  changes  in  mean  water  level.  The  change  in  water  surface  elevation 
will  change  the  width  of  the  surf  zone  and,  therefore,  the  profile  of  the 
longshore  current.  High-quality  correlations  between  the  mean  water  surface 
elevation  and  wind  velocity  for  the  nearshore  zone  of  an  open  coast  are  not 
known  to  the  authors,  but  numerical  surge  models  incorporating  wind  forcing 
can  well  produce  water  marks  occurring  during  storm  inundation.  An  example  of 
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time -dependent  change  in  water  surface  elevation  was  given  by  DeVries  and 
Amorocho  (1980)  for  both  ends  of  a  9 . 7-km-long,  37-m-wide,  and  8.4-m-deep 
section  of  the  California  Aqueduct,  which  they  successfully  simulated  by  using 
a  numerical  model.  Their  study  developed  a  predictive  equation  for  the  wind 
drag  coefficient  (Amorocho  and  DeVries  1980)  that  has  been  used  in  coastal 
applications,  but  the  expression  may  not  be  suitable  at  higher  wind  speeds  due 
to  the  restricted  fetch  in  their  experiments.  Amorocho  and  DeVries  report 
that  setup  in  the  aqueduct  was  negligible  if  the  wind  speed  at  the  2-m 
elevation  fell  below  approximately  4  m/sec. 

35.  Wind  and  wave  breaking.  The  local  wind  can  modify  the  characteris¬ 
tics  of  depth-limited  breaking  waves,  such  as  breaker  height  and  type  (as 
plunging  or  spilling).  The  only  quantitative  work  on  this  subject  known  to 
the  authors  is  that  of  Douglas  and  Weggel  (1988).  For  a  certain  type  of  wave 
condition  in  their  physical  model  experiment,  Douglas  and  Weggel  found  that  a 
wind  blowing  offshore  decreased  the  surf  zone  width  by  43  percent  due  primari¬ 
ly  to  a  lessening  of  breaking  wave  depth.  Thus,  incipient  depth- limited 
breaking  waves  were  able  to  travel  farther  through  the  "support"  given  by  the 
offshore  wind.  This  phenomenon  has  not  gone  unnoticed  by  surfers  who  pursue 
longer  rides  during  days  with  an  offshore -directed  breeze.  Implications  of 
the  effect  of  the  cross -shore  component  of  the  wind  on  wave  breaking  and  the 
longshore  current  are  profound,  since  changes  in  breaking  wave  type  changes 
the  wave  energy  dissipation  and  wave  height  distribution  across  the  shore, 
thereby  changing  the  longshore  current  profile.  Additionally,  a  change  in 
surf  zone  width  due  to  the  wind  correspondingly  changes  the  width  of  the 
current  profile  and  its  location  in  the  nearshore.  The  effect  of  the  wind  on 
wave  breaker  type  and  other  breaking  wave  characteristics  is  not  included  in 
the  present  model. 


Analytic  Solutions 

36.  Starting  from  the  first  analytic  solutions  for  the  longshore 
current  profile  based  on  the  momentum  equations  with  wave  radiation  stresses, 
numerous  approximations  have  been  introduced  to  improve  these  solutions  and 
extend  them  to  more  general  situations.  Despite  the  mathematical  skill  and 
physical  insight  that  have  gone  into  development  of  these  models,  the  result 
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is  always  an  exact  mathematical  solution  of  an  incomplete  formulation  of  the 
problem,  as  is  seen  in  Table  1,  which  summarizes  the  capabilities  of  represen¬ 
tative  analytic  models.  It  is  clear  that  the  analytic  models  satisfy  only  a 
few  of  the  requirements  of  an  engineering  model  as  enumerated  in  Part  I. 

37.  With  proliferation  of  inexpensive  and  powerful  desk-top  computers, 
the  main  function  of  analytic  models  now  lies  in  promoting  understanding  of 
the  physics  of  the  phenomenon  through  compact  display  of  the  relationships 
between  the  various  governing  parameters  and  how  they  control  the  current 
profile.  The  more  elaborate  analytic  solutions  are  complex  and  require 
evaluation  of  computation- intensive  special  functions.  Also,  matching  condi¬ 
tions  are  required  since  solutions  must  be  obtained  separately  seaward  and 
shoreward  of  the  wave  break  point,  the  result  of  an  assumed  discontinuity  in 
the  wave  driving  force  at  the  break  point.  Numerical  solutions,  described  in 
the  next  section,  automatically  compute  across  the  breaker  line. 

38.  In  the  following,  selected  analytic  solutions  based  on  momentum 
flux  concepts  are  discussed  to  introduce  the  important  phenomena  and  to 
present  model  capabilities.  The  discussion  is  summarized  in  Table  1,  which 
compares  the  capabilities  of  the  analytic  models  ranging  from  the  oldest  to 
the  more  recent.  A  similar  table  will  be  given  for  numerical  models  in  the 
next  section.  Familiarity  with  basic  surf  zone  hydrodynamics  is  assumed,  and 
Part  III  can  be  consulted  for  a  review  of  the  governing  equations. 

39.  Discussion  is  omitted  of  closed- form  solutions  for  only  the  mean 
longshore  current  velocity,  which  are  obtained  by  averaging  the  current 
profile  across  the  surf  zone.  Such  calculations  have  been  presented  by 
Longue t- Higgins  (1970b)  for  oblique  wave  incidence,  and  Komar  (1975)  and 
Gourlay  (1976,  1982)  for  oblique  wave  incidence  and  longshore  variation  in 
wave  height  (and  in  wave  angle,  Keely  and  Bowen  1977). 

40.  To  obtain  closed- form  or  analytic  solutions,  truncated  power  series 
and  perturbation  techniques  are  universally  employed,  principally  to  linearize 
the  form  of  the  quadratic  bottom  friction  stress,  which  is  the  main  retarding 
force  of  the  current.  These  procedures  introduce  mathematical  limitations  on 
the  range  of  validity  of  the  solution,  since  the  depth  and  expansion  parame¬ 
ters  should  formally  remain  small.  Shallow-water  approximations  to  linear- 
wave  theory  must  also  be  used  to  eliminate  hyperbolic  functions  appearing  in 
the  equations.  Thus,  analytic  solutions  are  not  only  exact  solutions  of  a 
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Table  1 


Capabilities  of  Selected  Analytic  Longshore  Current  Models 


Model  by  Author* 


Capability 

B 

L-H 

L&D 

K&S 

D 

T&G 

S&K 

M&H 

Concave  beach 

-- 

-- 

-- 

Yes 

Barred  beach 

-- 

-- 

-- 

-- 

-- 

Large  wave  angles 

-- 

Yes 

Yes 

Yes 

Wave  refraction 

-- 

Yes 

Yes 

Yes 

Yes 

Nonuniform  wave 
height  decay 

-- 

-- 

-- 

-- 

Yes 

Linear  bottom 
friction 
(weak  current) 

Yes 

Yes 

Yes 

Yes 

Yes 

Yes 

Yes 

Yes 

Nonlinear  bottom 
friction  (strong 
current) 

Yes* ** 

Lateral  mixing 

Yes 

Yes 

Yes 

-- 

Yes 

Yes 

Wave -current 
interaction 

Yes 

Wind  driving 

-- 

-- 

-- 

-- 

-- 

-- 

Random  wave  height 

_ 

Yes 

_  _ 

-  - 

Random  wave 
direction 


*  Key  to  authors  and  models:  B  =  Bowen  (1969);  L-H  =  Longuet- 
Higgins  (1970a,  b) ;  L&D  =  Liu  and  Dalrymple  (1978);  K&S  =  Kraus 
and  Sasaki  (1979);  D  =  Dalrymple  (1980);  T&G  =  Thornton  and  Guza 
(1986);  S6cK  =  Smith  and  Kraus  (1987);  M&H  =  McDougal  and  Hudspeth 
(1989). 

**  Weak  current  and  strong  current  solutions  are  different  models. 
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limited-scope  problem,  but  also  a  mathematically  approximate  description  of 
the  included  processes. 

41.  With  reference  to  Table  1  and  for  comparison  of  the  capabilities  of 
both  analytic  and  numerical  models,  limitations  found  in  most  analytic  models 
may  be  summarized  as : 

a.  Small  incident  wave  angle  (approximately  less  than  20  deg  at 
the  breaker  line) . 

b.  Linearized  bottom  friction  stress. 

c.  Regular  waves  (fixed  wave  height,  direction,  and  period). 

d.  Plane-sloping  beach. 

e.  Wave  height  proportional  to  the  water  depth  in  the  surf  zone 
(for  broken  waves) . 

f.  No  refraction  in  the  surf  zone. 

g.  No  wind  forcing. 

Use  of  shallow-water  wave  theory,  although  reducing  computation  time,  is  not 
convenient  for  general  engineering  applications  because  the  surf  width  and 
depth  will  vary.  In  the  following  discussion,  improvements  removing  one  or 
more  restrictions  will  be  noted. 

42.  The  earliest  analytic  models  (Bowen  1969;  Lonquet-Higgins  1970a,  b) 
demonstrated  that  the  momentum-based  formulation  led  to  qualitatively  and 
quantitatively  correct  shapes  and  magnitudes  of  the  wave -induced  longshore 
current  profile.  The  leading-order  processes  were  identified  as  bottom 
friction  and  the  excess  momentum  flux  or  radiation  stresses  of  the  waves,  with 
lateral  mixing  acting  primarily  to  adjust  the  shape  of  the  profile.  The 
formulation  of  the  bottom  friction  force  was  validated  with  laboratory  data, 
and  the  order  of  magnitude  of  the  bottom  friction  coefficient  was  determined 
to  be  10"2.  In  these  and  in  almost  all  analytic  and  most  numerical  models, 
the  time-averaged  bottom  friction  stress  is  linearized  by  assuming  that  the 
magnitude  of  the  longshore  current  is  much  smaller  than  that  of  the  wave- 
induced  horizontal  water  particle  orbital  velocity.  Other  approximations  used 
in  the  original  models  of  Bowen  (1969)  and  Longuet-Higgins  (1970a,  b)  include 
restriction  of  the  waves  to  small  angles,  omission  of  wave  refraction  in  the 
surf  zone,  and  simplified  formulations  of  the  eddy  viscosity  coefficient. 

Also,  the  driving  force  of  wind  and  the  wave-current  interaction  were  omitted. 
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43.  In  the  Bowen  (1969)  model,  the  bottom  friction  stress  was  linear¬ 
ized  in  the  longshore  current  velocity,  and  both  the  wave  orbital  velocity  in 
the  bottom  friction  stress  and  the  eddy  viscosity  in  the  lateral  mixing  term 
were  constant.  The  resultant  solution  of  the  second-order  differential 
equation  for  the  longshore  current  as  a  function  of  the  distance  across  the 
shore  pertains  to  a  uniformly  (plane)  sloping  beach.  The  solution  is  ex¬ 
pressed  in  terms  of  modified  Bessel  functions  and  is  not  easily  evaluated 
without  a  computer. 

44.  Longuet-Higgins  (1970a)  solved  for  the  current  profile  inside  the 
surf  zone  by  assuming  no  lateral  mixing,  obtaining  a  triangular  distribution 
for  the  current  with  a  peak  at  the  breaker  line.  In  a  companion  paper 
(Longuet-Higgins  1970b) ,  lateral  mixing  was  included  with  an  eddy  viscosity 
coefficient  proportional  to  the  distance  offshore  and  the  wave  orbital  speed. 

A  linearized  bottom  friction  stress  was  used  with  the  wave  orbital  velocity 
expressed  by  shallow-water  wave  theory.  Solution  of  the  resultant  second- 
order  ordinary  differential  equation  for  the  current  resulted  in  a  power 
series  of  only  two  terms  inside  the  surf  zone  and  one  term  for  the  region 
seaward  of  the  breaker  line.  An  important  dimensionless  parameter  that 
emerged  controlled  the  shape  of  the  profile  and  can  be  interpreted  as  express¬ 
ing  the  relative  strength  of  lateral  mixing  and  bottom  friction. 

45.  A  surprising  result  in  longshore  current  modeling  has  been  that  a 
linearized  friction  term  provides  a  good  quantitative  solution  of  the  current 
profile  even  though  the  assumption  is  violated  in  the  verification  data. 
Thornton  and  Guza  (1986)  and  the  present  work  as  discussed  below  have  shown 
that  the  success  of  the  linear  model  in  describing  virtually  all  field  and 
laboratory  data  sets  owes  to  the  fact  that  the  bottom  friction  and  wave 
driving  force  dominate,  leading  to  an  approximate  correspondence  in  linear  and 
nonlinear  solutions  by  a  reduction  in  the  value  of  the  friction  coefficient  in 
the  nonlinear  friction  model  as  compared  with  the  linear  model.  Presumably, 
the  value  of  the  bottom  friction  coefficient  pertaining  to  the  nonlinear  model 
is  more  physically  correct. 

46.  Liu  and  Dalrymple  (1978)  examined  the  functioning  and  mathematical 
range  of  validity  of  the  linear  or  weak- current  and  nonlinear  or  strong- 
current  bottom  friction  formulations  through  closed- form  solutions  for  the 
current  profile  inside  the  breaker  line,  including  wave-induced  setup  but 
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omitting  lateral  mixing.  Weak-current  and  strong- current  approximations  were 
treated  independently  to  obtain  two  different  solutions.  Although  the  formal 
validity  of  the  linear  and  nonlinear  bottom  friction  approximations  was 
demonstrated,  their  results  did  not  explain  the  robustness  of  the  linear  model 
achieved  by  simple  scaling  of  the  value  of  the  friction  coefficient.  Also, 
the  weak- current  case  of  Liu  and  Dalrymple  contains  an  error  producing  values 
of  the  current  at  the  breaker  line  that  are  too  small  by  a  factor  of  about  0.5 
(Kraus  and  Sasaki  1979,  Baum  and  Basco  1986,  Miller  1987). 

47.  Kraus  and  Sasaki  (1979)  extended  the  Longuet-Higgins  (1970a,  b) 
solution  to  include  large  wave  angles,  refraction  in  the  surf  zone,  and  the 
improved  lateral  mixing  description  of  Madsen,  Ostendorf,  and  Reyman  (1978). 
The  model  was  verified  with  four  laboratory  cases  of  Mizuguchi ,  Oshima,  and 
Horikawa  (1978)  and  their  own  field  data.  By  matching  the  position  of  the 
maximum  current  velocity,  the  value  of  the  bottom  friction  coefficient  found 
in  the  field  verification  was  0.0061,  smaller  than  the  values  ranging  from 
0.011  to  0.024  found  for  the  laboratory  cases. 

48.  Dalrymple  (1980)  investigated  the  wave-current  interaction  in  the 
surf  zone  associated  with  wave  refraction  and  larger  incident  wave  angles  in 
an  analytic  model  that  omitted  lateral  mixing.  The  wave -current  interaction 
was  included  in  Snell's  Law  governing  refraction  (wave  direction),  but  was 
omitted  in  calculating  wave  height  decay  because  a  simple  linear  decrease  in 
wave  height  in  constant  proportion  to  the  water  depth  was  assumed.  Dalrymple 
found  that  the  interaction  lowered  the  magnitude  of  the  current  and  noted  that 
conservation  of  the  force  on  the  bottom  by  the  waves  implies  that  the  width  of 
the  surf  zone  must  correspondingly  increase.  Dalrymple  concluded  that  the 
wave-current  interaction  by  refraction  was  significant  (order  of  10  percent) 
only  for  large  incident  wave  angles.  Similarly,  Thornton  and  Guza  (1986) 
found  that  the  refractive  effect  of  the  mean  current  was  less  than  1  percent 
for  their  calculated  cases. 

49.  Thornton  and  Guza  (1986)  derived  an  analytic  solution  under  the 
conditions  of  a  plane  beach,  small  wave  angle,  no  mixing,  and  linear  bottom 
friction,  but  using  an  analytic  random  wave  height  model.  Their  paper  intro¬ 
duces  randomness  in  the  wave  height  in  the  longshore  current  calculation  in  an 
efficient  way  and  compares  model  predictions  to  a  high-quality  field  data  set. 
The  model  successfully  reproduced  the  measurements  of  one  case  of  their  field 
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measurements  (2  February  1980)  by  adjustment  of  the  value  of  the  bottom 
friction  coefficient.  Thornton  and  Guza  also  present  numerical  models  with 
linear  and  nonlinear  bottom  friction  representations,  discussed  in  the  next 
section. 

50.  Smith  and  Kraus  (1987)  extended  the  Kraus  and  Sasaki  (1979)  model 
to  include  a  power  law  form  of  the  wave  height  decay  (concave  up,  concave 
down,  or  linear  with  water  depth)  depending  on  an  empirical  coefficient 
related  to  the  breaker  index  (ratio  of  wave  height  and  water  depth  at  incipi¬ 
ent  breaking)  and  beach  slope.  The  power  law  provides  a  more  realistic 
representation  of  the  wave  height  distribution  across  the  shore  than  a  simple 
linear  dependence  on  water  depth  and  was  found  to  control  the  shape  of  the 
current  profile  in  the  surf  zone  similar  to  lateral  mixing.  This  result 
indicates  that  more  accurate  representation  of  wave  height  decay  will  change 
the  values  of  the  bottom  friction  coefficient  and  lateral  mixing  strength. 
Miller  (1987)  developed  a  similar  capability  in  a  more  complex  analytic 
solution  using  an  exponential  wave  height  decay.  Smith  and  Kraus  (1988) 
investigated  wind  driving  in  an  analytic  wave-  and  wind- generated  longshore 
current  model,  but  found  it  infeasible  because  of  a  discontinuity  in  the 
solution  at  the  shoreline. 

51.  Leont'ev  (1988)  developed  analytic  solutions  for  random  wave  height 
decay  in  the  surf  zone  and  the  resultant  longshore  current  for  both  plane - 
sloping  beaches  and  beaches  having  a  bar  and  trough.  Previous  analytic  models 
of  wave  breaking  and  longshore  currents  were  also  reviewed.  Leont'ev  verified 
his  random  height  model  against  the  bar-trough  laboratory  data  of  Battjes  and 
Janssen  (1978) . 

52.  McDougal  and  Hudspeth  (1989)  developed  an  analytic  solution  for  a 
planar  slope  on  the  foreshore  joining  to  a  concave  bottom  profile.  A  concave- 
sloped  beach,  which  is  more  realistic  representation  of  the  typical  beach 
profile,  was  found  to  have  the  peak  in  the  longshore  current  located  more 
shoreward  than  on  a  planar  slope.  Their  solution  is  relatively  complicated 
and  is  expressed  in  terms  of  special  functions  involving  matching  conditions 
both  at  the  slope  intersection  and  the  breaker  line. 
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Numerical  Solutions 


53.  A  large  number  of  numerical  models  of  the  wave-induced  nearshore 
circulation  have  been  developed  since  the  first  comprehensive  model  of  this 
type  introduced  by  Noda  (1974) .  These  models  solve  for  the  horizontal 
circulation  and  water  surface  elevation  on  a  2-D  grid  by  using  either  a 
finite-difference  or  finite -element  solution  approach.  Two-dimensional  models 
are  extremely  expensive  to  run  as  compared  with  one -dimensional  (1-D)  models 
such  as  presented  here,  but  are  necessary  to  calculate  wave  and  current 
conditions  on  irregular  sea  bottoms.  Two-dimensional  models  are  computation 
expensive  for  two  reasons.  First,  calculations  must  be  performed  for  Nz  grid 
points  as  opposed  to  N  points  for  a  1-D  model,  where  N  is  the  representative 
number  of  grid  points  in  the  x-  and  y- directions .  Second,  2-D  models  typical¬ 
ly  reach  the  steady-state  solution  by  stepping  through  time,  requiring  on  the 
order  of  200  to  500  iterations  over  the  grid  for  a  general,  new  condition. 
("Hot  starts"  can  be  made  if  the  solution  for  a  similar  condition  is  avail¬ 
able,  which  reduces  the  number  of  iterations  somewhat,  although  instabilities 
can  sometimes  occur  with  this  procedure.)  In  contrast,  1-D  solutions  can 
operate  directly  with  the  steady-state  equations  of  motion,  requiring  as 
little  as  one  pass  or  iteration  on  N  grid  points  with  the  particular  solution 
scheme  for  a  linearized  model  or,  for  example,  typically  4  to  8  iterations  for 
the  present  model  in  nonlinear  form.  In  the  following,  discussion  will  focus 
on  representative  models  listed  in  Table  2. 

54.  Thornton  (1970)  developed  a  numerical  model  that  included  lateral 
mixing  and  the  bottom  friction  coefficient  of  Jonsson  (1966).  Although  the 
model  in  principle  could  be  applied  to  an  arbitrary  bottom  topography,  in 
practice  it  was  limited  by  a  simple  wave  height  decay  model  that  depends  on  a 
constant  wave -height  to  water -depth  ratio  in  the  surf  zone.  There  was  only 
one  empirical  parameter,  the  bottom  roughness,  and  model  predictions  agreed 
with  laboratory  measurements  and  approximate  field  measurements.  Details  of 
the  numerical  solution  method  were  not  given. 

55.  Jonsson,  Skovgaard,  and  Jacobsen  (1974)  developed  a  numerical  model 
with  a  sophisticated  numerical  solution  scheme  that  is  probably  not  warranted 
for  engineering  application  considering  the  empiricism  of  the  formulation  and 
the  accuracy  of  input  data.  An  interesting  feature  of  the  development  was  use 
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Table  2 


Capabilities  of  Selected  Numerical  Longshore  Current  Models 


Model  by  Author* 


Capabilitv 

T 

JSJ 

V 

WTG 

BB 

H 

S 

Present 

Model 

Barred  beach 

-- 

-- 

Yes 

Yes 

? 

Yes 

Large  wave  angles 

Yes 

Yes 

Yes 

Yes 

Yes 

Yes 

Yes 

Yes 

Wave  refraction 

Yes 

Yes 

Yes 

Yes 

Yes 

Yes 

Yes 

Yes 

Nonuniform  wave 
height  decay 

-- 

-- 

-- 

Yes 

-- 

? 

Yes 

Nonlinear  bottom 
friction 

Yes 

Yes 

-- 

Yes 

Yes 

Yes 

Lateral  mixing 

Yes 

Yes 

Yes 

Yes 

Yes 

Yes 

Wave -current 
interaction 

Yes 

-- 

Yes 

Yes 

Wind  driving 

-- 

Yes 

Yes 

Random  wave  height 

Yes 

Yes 

Yes 

Estimate  of 

? 

Med 

Med 

Low 

Med 

Low 

Med 

High 

calc,  efficiency 


*  Key  to  authors  and  models:  T  =  Thornton  (1970);  JSJ  =  Jonsson, 
Skovgaard,  and  Jacobsen  (1974);  V  =  Visser  (1984a);  Wu,  Thornton, 
and  Guza  (1985);  BB  =  Baum  and  Basco  (1986);  H  -  Hubertz  (1987); 

S  =  Southgate  (1989) . 


of  a  combined  wave  and  current  friction  factor. 

56.  Visser  (1984a,  b)  presented  a  numerical  model  "aimed  to  be  more 
reliable  than  existing  ones"  through  calibration  with  his  careful  laboratory 
measurements  (Visser  1982)  involving  seven  cases.  Four  of  the  cases  are  used 
in  model  testing  in  the  present  study.  The  Visser  model  includes  a  nonlinear 
bottom  friction  term  formulated  for  combined  waves  and  current,  lateral 
mixing,  and  refraction  in  the  surf  zone.  It  has  limitations  in  using 
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shallow-water  linear-wave  theory  and  constant  wave-height  to  water-depth  ratio 
in  the  surf  zone.  Rather  than  allow  wave  energy  dissipation  to  begin  at  the 
break  point  (point  of  highest  wave  height) ,  the  dissipation  is  started  at  the 
measured  plunge  point  based  on  his  observation  that  measured  setdown  was 
greatest  at  the  plunge  point,  not  at  the  break  point.  The  numerical  model 
gave  an  overall  excellent  fit  to  the  data  with  the  same  value  of  the  empirical 
mixing  coefficient,  individually  determined  bottom  friction  coefficients 
(roughness  lengths),  and  the  measured  plunge  point.  Details  of  the  numerical 
solution  scheme  were  not  given  except  that  iteration  was  used. 

57.  Wu,  Thornton,  and  Guza  (1985)  applied  their  2-D  finite -element 
model  to  simulate  the  alongshore  current  system  measured  at  Leadbetter  Beach, 
California  (Thornton  and  Guza  1986,  1989).  These  measurements  are  also  used 
here.  Their  model  included  linearized  bottom  friction  terms,  lateral  mixing 
after  Longuet-Higgins  (1970b),  and  nonlinear  convective  terms.  The  wave  model 
in  the  surf  zone  used  a  constant  wave-height  to  water-depth  ratio,  and  the 
measured  root-mean- square  (rms)  wave  height  was  input  as  a  monochromatic 
regular  wave.  Fitting  parameters  were  the  bottom  friction  coefficient,  mixing 
coefficient,  and  a  wave-height  to  water-depth  ratio  of  0.44  determined  in  one 
case,  which  was  held  constant  through  the  other  four  cases.  Use  of  the  small 
value  of  0.44  is  inconsistent  with  the  treatment  of  the  rms  wave  as  a  mono¬ 
chromatic  wave,  and  this  subject  is  pursued  further  in  Part  IV.  The  mean 
value  of  the  (linear)  bottom  friction  coefficient  was  0.01.  Running  of  the 
model  with  and  without  the  nonlinear  convective  terms  gave  about  the  same 
agreement  with  the  measurements  by  modification  of  the  calibration  parameters. 

58.  Baum  and  Basco  (1986)  presented  a  numerical  model  similar  to  the 
model  developed  here  in  emphasizing  calculation  of  the  longshore  current 
profile  on  a  bar  and  trough  bottom.  The  wave  decay  model  of  Dally  (1980)  and 
Dally,  Dean,  and  Dalrymple  (1985),  which  allows  multiple  wave  breaking  and 
reformation,  provided  the  wave  driving  forces.  Shallow-water  linear-wave 
approximations  were  used  everywhere  in  the  calculation  domain.  Other  features 
were  large-wave  angle,  linearized  bottom  friction  stress,  and  lateral  mixing. 
Model  predictions  agreed  well  with  the  laboratory  data  of  Mizuguchi,  Oshima, 
and  Horikawa  (1978)  for  a  plane-sloping  beach,  and  reasonable  agreement  was 
found  with  the  field  data  of  Kraus  and  Sasaki  (1979)  for  an  irregular  profile. 
Numerical  efficiency  for  this  linear  model  is  not  optimal  because  iteration  is 


28 


used  to  determine  the  mean  water  surface  elevation  and  an  iterative  "shooting" 
method  to  calculate  the  current . 

59.  Recognizing  the  importance  of  wave  randomness  and  the  wind  in 
determining  the  longshore  current  profile  in  the  field,  Hubertz  (1987) 
developed  a  model  incorporating  these  effects  to  provide  a  more  realistic  and 
accurate  alternative  to  handbook- type  calculations.  Wave  randomness  was 
calculated  with  a  spectral  wave  propagation  model  that  transforms  the  spectrum 
accounting  for  changes  in  direction,  height,  and  peak  frequency.  The  drag 
coefficient  of  Amorocho  and  DeVries  (1980)  was  used  in  the  wind  driving  term. 
The  model  contains  a  quadratic  bottom  friction  law,  but  lateral  mixing  was 
omitted  and  a  constant  wave-height  to  water-depth  ratio  applied  in  the  surf 
zone.  Model  calculations  were  verified  with  field  data  involving  both  waves 
and  wind,  including  the  current  seaward  of  the  surf  zone  arising  mainly  from 
the  wind.  The  model  was  developed  as  a  special  case  of  a  2-D  model  with  a 
time-stepping  solution  to  the  steady  state;  therefore,  several  hundred 
iterations  are  required  to  obtain  the  solution. 

60.  The  most  recent  numerical  longshore  current  model  found  in  the 
literature  review  is  part  of  a  larger  1-D  modeling  system  for  calculating 
waves,  wave-induced  longshore  current,  and  tidal  current  (Southgate  1989). 

The  model  was  developed  with  the  same  pragmatic  philosophy  as  the  present 
study,  that  engineering  situations  may  involve  thousands  of  runs  and  large 
quantities  of  data,  necessitating  use  of  1-D  models.  Features  included  in  the 
model  are  random  wave  breaking  through  the  truncated  Rayleigh  distribution 
method  of  Battjes  and  Janssen  (1978),  nonlinear  bottom  friction  law,  wave- 
current  interaction,  and  inclusion  of  the  tidal  current  as  calculated  from 
another  model  or  data.  Lateral  mixing  was  evidently  not  taken  into  account  to 
increase  calculation  efficiency,  and  wind  driving  was  not  included.  It  is 
believed  that  the  model  as  presented  is  not  applicable  to  barred  beaches.  In 
the  monochromatic  (single-wave)  mode,  agreement  of  the  current  prediction  for 
Visser  (1982)  Case  1  was  not  good  because  of  the  omission  of  lateral  mixing. 

It  was  argued  that  for  field  application  use  of  random  wave  heights  decreased 
the  need  for  lateral  mixing.  Details  of  the  numerical  solution  scheme  were 
not  given  but,  because  of  the  necessity  of  performing  time  integration  of  the 
bottom  friction  stress  and  an  apparent  absence  of  a  skillful  first 
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approximation  for  the  iteration,  efficiency  of  the  Southgate  model  is  judged 
to  be  medium. 

61.  Finally,  this  review  concludes  with  mention  of  two  longshore 
current  models  that  employed  nonlinear -wave  theory  under  simple  conditions  of 
a  plane-sloping  beach,  those  of  James  (1974)  and  Pechon  (1987) .  James 
employed  a  hybrid  hyperbolic  wave  theory  to  calculate  waves  and  current  in  the 
surf  zone,  but  such  an  approach  no  longer  seems  reasonable  with  accurate  and 
efficient  calculation  schemes  for  cnoidal  waves.  Pechon  (1987)  used  the 
shallow-water  wave  theory  of  Serre  (1953;  according  to  Pechon  1987)  and  gave 
examples  of  differences  in  wave  height  transformation,  radiation  stresses,  and 
longshore  current  calculated  with  linear-  and  nonlinear -wave  theory.  Because 
only  one  comparison  case  with  the  Visser  (1982)  data  was  shown,  it  is  diffi¬ 
cult  to  judge  superiority  of  results.  A  problem  in  use  of  nonlinear -wave 
theories  is  matching  at  the  boundary  between  shallower  and  deeper  water,  where 
discontinuities  in  some  quantities  will  occur  (Hardy  and  Kraus  1988).  Also, 
it  has  not  been  proved  that  any  nonlinear -wave  theory  is  superior  to  linear- 
wave  theory  in  the  surf  zone.  Once  special  conditions  are  imposed,  such  as 
use  of  nonlinear -wave  theory  to  calculate  shoaling,  but  not  refraction, 
modification  of  the  group  speed  inside  the  surf  zone,  but  not  outside,  and 
other  such  inconsistent  simplifications,  generality  of  the  model  is  probably 
degraded. 
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PART  III:  WAVE  MODEL 


62.  This  chapter  describes  the  wave  model  used  to  drive  the  longshore 
current  model.  Wave  model  selection  and  development  are  first  discussed, 
followed  by  introduction  of  the  governing  equations  and  a  description  of  the 
numerical  implementation.  Randomization  of  the  wave  model  output  is  accom¬ 
plished  using  a  Monte-Carlo  input  simulation  technique,  making  it  possible  to 
calculate  properties  of  the  wave  height  probability  density  function  at  any 
location  across  the  shore  on  a  multiple  bar  and  trough  profile.  The  wave 
model  is  verified  with  laboratory  and  field  measurements  in  Part  IV.  Addi¬ 
tional  wave  model  verification  is  presented  in  conjunction  with  the  longshore 
current  simulations  given  in  Part  VI,  as  interpretation  of  the  current 
calculation  results  requires  consideration  of  the  wave  input. 

Overview  of  Wave  Model 

63.  Calculation  of  the  longshore  current  as  a  function  of  distance 
across  the  shore  requires  knowledge  of  the  driving  wave  characteristics  at 
intervals  inside  and  outside  the  surf  zone.  As  a  predictive  tool,  the  wave 
model  must  calculate  the  change  in  height,  direction,  and  length  of  the  waves 
as  they  propagate  onshore  starting  from  wave  conditions  specified  at  the 
seaward  end  of  the  grid.  Such  a  calculation  involves  the  major  wave  transfor¬ 
mation  processes  of  shoaling,  refraction,  breaking  with  energy  dissipation, 
and  wave  reformation.  Wave  diffraction  does  not  enter  in  the  applications 
considered  in  this  report  since  it  produces  wave  height  and  direction  gradi¬ 
ents  perpendicular  to  the  direction  of  wave  propagation,  violating  the  basic 
assumption  of  longshore  uniformity  used  in  the  present  model.  Wave  dissipa¬ 
tion  seaward  of  the  surf  zone  is  neglected  for  the  sake  of  computational 
speed,  but  could  be  easily  incorporated  if  warranted,  such  as  for  calculation 
of  waves  and  currents  over  a  muddy  bottom. 

64.  In  the  development  of  the  wave  model,  certain  requirements  were 
imposed  to  define  a  numerically  stable,  cost-effective  engineering  model  with 
minimal  computer  memory  demands : 

a.  The  model  should  accept  arbitrary  wave  input  (time  series  of 
representative  wave  height,  period,  and  direction)  and  complex 
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bottom  topographies  (which  exhibit  uniformity  alongshore  as,  for 
example,  beaches  with  multiple  bars  and  troughs) . 

b.  The  calculation  accuracy  of  the  wave  transformation  should  be 
compatible  with  engineering  data  available  for  the  wave  input 
and  calculation  accuracy  of  the  current  distribution. 

c.  The  model  should  contain  a  minimum  number  of  empirical  parame¬ 
ters,  preferably  with  values  that  have  been  well  established 
with  data  for  general  wave  and  bottom  topography  conditions . 

d.  The  execution  time  of  the  model  should  be  sufficiently  short  and 
memory  requirements  minimal  to  allow  it  to  run  efficiently  on  a 
personal  computer. 

65.  One  representation  of  the  wave  transformation  that  satisfies  the 
criteria  is  linear-wave  theory  for  calculating  shoaling  and  refraction  seaward 
of  the  break  point,  connecting  to  a  modified  version  of  the  breaker  decay 
model  developed  by  Dally  (1980)  (see  also  Dally,  Dean,  and  Dalrymple  1985)  to 
calculate  wave  properties  in  the  surf  zone.  Linear-wave  theory  has  been 
successful  in  many  applications  for  predicting  wave  propagation  even  under 
conditions  where  it  is  formally  violated,  such  as  in  the  vicinity  of  the  break 
point  and  in  the  surf  zone,  and  it  requires  minimal  computer  memory  and 
execution  time  as  compared  with  nonlinear-wave  theories.  At  the  present  state 
of  knowledge  on  surf  zone  dynamics,  it  seems  difficult  to  significantly 
improve  the  predictive  capabilities  and  reliability  of  engineering  wave  models 
using  more  complicated  wave  theories,  although  efficient  calculation  routines 
are  available  for  the  region  seaward  of  the  break  point  (e.g.,  Kraus,  Cialone, 
and  Hardy  1987;  Hardy  and  Kraus  1988). 

66.  The  Dally  (1980)  breaker  decay  model  is  one  of  several  models 
available  for  calculating  wave  transformation  in  the  surf  zone.  Similar 
models  have  been  developed  by  Goda  (1975),  Battjes  and  Janssen  (1978), 
Mizuguchi  (1980) ,  Thornton  and  Guza  (1983) ,  and  Svendsen  (1984) .  The  main 
difference  among  these  models  is  the  formulation  of  the  energy  dissipation 
term,  where  one  common  approach  uses  the  analogy  between  wave  breaking  and  a 
moving  hydraulic  bore  for  estimating  the  rate  of  energy  dissipation 

(Le  Mehaute  1962).  Some  models  employ  a  probabilistic  approach  (Goda  1975, 
Battjes  and  Janssen  1978,  Thornton  and  Guza  1983)  for  describing  the  transfor¬ 
mation  of  a  representative  statistical  wave  height,  usually  the  rms  wave 
height.  A  discussion  of  the  differences  between  probabilistic  models  and 
deterministic  models  is  given  later  in  this  chapter. 
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67 .  The  Dally  model  was  selected  for  use  in  the  present  study  because  of 
the  extensive  verification  that  has  been  made  for  a  variety  of  wave  conditions 
and  bottom  topographies.  Both  small-scale  (Dally,  Dean,  and  Dalrymple  1985) 
and  large-scale  laboratory  (Larson  and  Kraus  1989)  measurements  have  been 
employed  in  the  verification  process,  as  well  as  field  measurements  (Ebersole 
1987,  Dally  1990).  An  appealing  feature  of  the  model  is  the  demonstrated 
reliability  of  the  values  of  its  two  empirical  parameters,  although  a  small 
and,  as  yet,  unquantified  dependence  on  beach  slope  has  been  noted  (Dally, 
Dean,  and  Dalrymple  1985;  Larson  and  Kraus  1989).  Thus,  a  satisfactory 
description  of  wave  height  transformation  in  the  surf  zone  using  average 
parameter  values  is  expected  without  additional  calibration  effort.  The  Dally 
model  is  formulated  in  terms  of  wave  energy  flux  and  is,  therefore,  general. 

In  applications  to  date,  including  the  present  report,  the  flux  is  expressed 
by  means  of  linear-wave  theory. 

68.  The  question  arises  whether  other  wave-related  quantities  are  accu¬ 
rately  described  by  the  model,  in  particular,  the  radiation  stresses,  which 
are  of  importance  in  calculating  the  longshore  current,  even  if  the  wave 
height  is  well  predicted.  Roelvink  and  Stive  (1989)  distinguished  between 
turbulent  energy  production  and  dissipation,  the  latter  governing  changes  in 
radiation  stresses  and  thus  quantities  such  as  the  setup  and  the  longshore 
current.  This  topic  is  discussed  in  Part  VI  where  longshore  current  model 
predictions  are  compared  with  laboratory  measurements  made  by  Visser  (1982) . 
Indirect  evaluation  of  the  model  can  be  made  by  comparing  measured  and 
calculated  quantities,  such  as  the  mean  water  surface  elevation  and  longshore 
current  profile,  which  depend  on  different  components  of  the  radiation  stress. 

69.  In  the  model,  wave  transformation  seaward  of  the  break  point  for  an 
individual  wave  is  calculated  using  the  wave  energy  flux  equation  without 
bottom  or  surface  dissipation  (although  dissipation  could  be  included  if 
warranted) ,  with  quantities  determined  from  linear-wave  theory.  Once  wave 
breaking  takes  place,  a  dissipation  term  representing  breaking  is  activated  in 
the  energy  flux  equation.  Wave  refraction  throughout  the  calculation  domain 
is  given  by  Snell's  Law  under  the  condition  of  straight  and  parallel  bottom 
contours.  An  empirical  relationship  is  used  to  predict  the  critical  ratio 
between  wave  height  and  water  depth  at  incipient  breaking.  The  mean  displace¬ 
ment  of  the  water  surface  (setup  or  setdown)  is  calculated  using  the 
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cross-shore  momentum  equation.  Importantly,  the  model  has  the  desirable 
feature  of  allowing  wave  reformation;  wave  breaking  can  cease,  and  a  wave  can 
again  begin  to  shoal  if  the  wave -height  to  water- depth  ratio  falls  below  a 
critical  value. 

70.  All  equations  in  the  wave  model  involve  variables  regarded  as  time 
averages;  that  is,  changes  in  the  offshore  wave  conditions  must  occur  suffi¬ 
ciently  slowly  to  allow  a  quasi- stationary  approach.  In  principle,  the 
restriction  of  slowly  varying  conditions  should  apply  to  the  calculation  of 
setup  and  setdown  as  well  as  the  current,  because  establishment  of  a  mean 
water  surface  elevation  requires  a  certain  number  of  identical  waves  before 
the  time -averaged  momentum  equations  are  valid.  The  wave  energy  flux  equation 
has  been  shown  to  be  applicable  for  individual  waves  propagating  across  the 
surf  zone  in  the  field  (Ebersole  1987)  ,  exhibiting  a  surprising  degree  of 
robustness  that  works  in  favor  of  the  approach  taken  here. 

71.  In  a  laboratory  environment  with  monochromatic  (constant  period)  and 
regular  (constant  height)  waves,  or  in  the  field  when  clean  swell- type  wave 
conditions  occur,  the  choice  of  a  representative  wave  height  is  unambiguous, 
and  the  break  point  will  be  relatively  stable  on  a  fixed  bed.  However,  if 
large  variability  in  the  wave  conditions  exists  (which  is  typically  the  case 
in  the  field  with  waves  of  different  generating  sources,  characteristics,  and 
breaking  locations),  it  is  more  difficult  to  determine  a  representative  wave 
height,  period,  and  direction.  A  probabilistic  approach  is  thus  reasonable 
for  describing  wave  transformation  in  the  surf  zone  under  these  conditions. 

72.  There  are  two  generic  approaches  for  probabilistically  describing 
surf  zone  wave  properties:  integration  of  the  wave  energy  flux  equation  with 
respect  to  an  existing  probability  density  function  (characteristics  of  the 
probability  density  function  must  be  known  at  each  point  across -shore,  such  as 
the  percentage  of  broken  waves),  and  simulation  of  the  transformation  of  each 
wave  component  from  a  known  probability  density  function  prescribed  in  deeper 
water.  The  latter  approach  is  computationally  more  intensive  and  rests  on  the 
validity  of  linearly  superimposing  the  response  from  each  wave  component,  but 
will  yield  information  about  the  local  wave  height  probability  density 
function  without  a  priori  assumptions  about  the  shape  of  the  function.  The 
individual -component  approach  is  also  expected  to  more  accurately  describe  the 
wave  height  distribution  on  an  irregular  bottom  where  wave  reformation  can 
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occur  and  both  breaking  and  reformed  waves  coexist.  Dally  (1990)  has  recently 
developed  closed- form  solutions  describing  random  wave  breaking  on  a  planar 
beach  using  the  individual -component  approach,  and  these  solutions  reproduce 
trends  reported  in  laboratory  and  field  measurements.  The  analytic  work  of 
Dally  for  planar  beaches  both  supports  and  complements  the  computer- intensive 
approach  taken  here  for  barred  beach  profiles. 

73.  A  strict  limitation  of  the  wave  and  longshore  current  models 
presented  in  this  report  arises  from  the  requirement  of  alongshore  uniformity 
in  forcing  conditions.  This  limitation  allows  a  1-D  approach  to  be  taken  in 
which  the  wave  (current)  distribution  along  a  single  cross -shore  transect  is 
sufficient  to  describe  the  local  wave  (current)  conditions.  In  situations 
violating  longshore  uniformity,  such  as  in  the  vicinity  of  structures  where 
significant  three-dimensionality  in  waves  occurs  or  in  areas  having  highly 
irregular  bathymetry,  the  present  model  may  not  be  suitable. 

Governing  Equations  for  the  Wave  Model 

74.  For  the  1-D  case  with  straight  and  parallel  bottom  contours  along¬ 
shore  (but  allowing  the  depth  to  be  nonmonotonic  with  distance  offshore,  such 
as  in  the  case  of  multiple  bars  and  troughs) ,  wave  characteristics  across  the 
profile  are  determined  by  four  fundamental  equations: 

a.  The  wave  energy  flux  equation. 

b.  The  cross -shore  momentum  equation. 

c.  The  wave  number  equation  (Snell's  Law). 

d.  The  wave  dispersion  relation. 

75.  Figure  1  is  a  definition  sketch  introducing  the  axes  convention  and 
notation  to  be  used  in  equations  governing  wave  transformation  and  the 
longshore  current.  Angle  definitions  giving  positive  orientations  are  also 
shown.  In  particular,  the  symbol  8  is  used  in  all  governing  equations  to 
denote  the  incident  wave  angle.  However,  the  sign  convention,  although 
common,  is  not  intuitive  because  a  positive  angle  will  generate  a  longshore 


For  convenience,  symbols  and  abbreviations  are  listed  in  the  Notation 
(Appendix  C) . 
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current  directed  along  the  negative  y-axis.  Therefore,  as  input  to  the 
numerical  model,  the  mirror-image  angle  &'  is  used,  which  is  defined  as 


0'  =  -6 


(1) 


cjD  ANGLE  OF  WIND 
0  ANGLE  OF  WAVE  RAY 


Figure  1.  Definition  sketch  for  wave  and  wind  direction 
Energy  flux 

76.  By  generalizing  the  breaker  decay  model  of  Dally  (1980)  originally 
developed  for  normal  wave  incidence  (see  Ebersole,  Cialone,  and  Prater  1986 
for  a  similar  treatment),  in  the  absence  of  currents,  the  2-D  equation  for 
conservation  of  energy  flux  incorporating  energy  dissipation  produced  by  wave 
breaking  is , 


(F  cos0)  +  -A  (F  sin0)  =  4  (F  -  Fs) 
ox  dy  d  s 


(2) 
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where 


x  -  cross -shore  coordinate  positive  in  the  seaward  direction  (m) 
F  =  wave  energy  flux  (Nm/m/sec) 

6  =  wave  angle  with  respect  to  shore  normal 
y  =  longshore  coordinate  (m) 

k  =  wave  decay  coefficient  ( nondimens ional  empirical  coefficient) 
d  =  total  water  depth  (m) 

Fs  =  stable  wave  energy  flux  (Nm/m/sec) 

The  wave  energy  flux  is  given  by, 


F  =  ECg 


(3) 


where 

E  =  wave  energy  density  (Nm/m2) 

Cg  =  wave  group  speed  (m/s) 

The  wave  energy  density  is  written  using  linear-wave  theory  as, 

E-±pcrH2  (4) 


where 

p  =  density  of  water  (kg/m3) 

g  =  acceleration  produced  by  gravity  (m/sec2) 

H  =  wave  height  (m) 

The  total  water  depth  d  is  the  sum  of  the  still -water  depth  h  and  the  change 
in  mean  water  level  rj  produced  by  waves  and  wind: 

d  =  h  +  T)  (5) 


77.  Assuming  the  waves  to  be  uniform  alongshore  and  the  bottom  contours 
to  be  straight  and  parallel,  Equation  2  reduces  to: 
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-0-  < F  cos0)  =  ~  Fa) 

dx  d  B 


(6) 


Equation  6  is  the  main  equation  used  in  this  report  to  describe  the  transfor¬ 
mation  of  wave  height  in  the  absence  of  a  current. 

78.  If  a  mean  current  is  present,  there  is  an  interaction  between  the 
waves  and  the  current,  and  an  additional  term  appears  in  the  wave  energy  flux 
equation  (Equation  6)  to  express  the  work  done  by  the  radiation  stress  on  the 
current.  Inclusion  of  the  wave-current  interaction  gives  (Mei  1983), 


-4-  (F  COS0) 
dx 


s*y%--d'F 


Fs) 


(7) 


where 

S^y  =  shoreward- directed  radiation  stress  alongshore  (N/m) 

V  =  longshore  current  speed  (m/sec) 

79,  Because  of  the  inclusion  of  the  dV/dx  term  in  Equation  7,  the  wave 
and  current  distributions  must  be  solved  for  simultaneously  (iteratively)  if 
waves  contribute  to  current  generation.  The  importance  of  the  radiation 
stress  or  wave-current  interaction  term  for  the  wave  height  and  the  across- 
shore  distribution  of  the  longshore  current  will  be  discussed  in  Part  VI . 

80.  The  radiation  stress  component  is  given  by  linear-wave  theory 

as , 

sxy  =  -~^p9HZsin2Q 


where 


n 


(9) 
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81.  In  Equation  7,  the  wave  decay  coefficient  k  controls  the  rate  of 
energy  dissipation,  whereas  the  stable  energy  flux  Fs  controls  the  amount  of 
energy  dissipation  necessary  for  stable  waves  to  exist  once  breaking  is 
initiated.  The  "stable  wave  condition"  refers  to  a  state  where  energy 
dissipation  produced  by  wave  breaking  ceases,  permitting  waves  to  reform. 

The  stable  energy  flux  is  expressed  as, 

F„  =  EeCg  (10) 


where  Es  is  the  wave  energy  density  corresponding  to  a  wave  that  has  decayed 
to  a  stable  form.  The  stable  wave  energy  flux  corresponds  to  a  stable  wave 
height  Hs,  which  is  a  function  of  the  water  depth  (Horikawa  and  Kuo  1966, 
Dally  1980), 

Hs  =  Yd  (11) 


where  T  is  a  nondimens ional  empirical  coefficient. 

82.  Two  empirical  coefficients,  therefore  enter  in  the  breaker  decay 
model,  k  and  T.  Values  of  these  coefficients  have  been  found  to  vary  little 
over  a  wide  range  of  wave  and  bottom  conditions  in  the  laboratory  and  on 
natural  beaches.  Dally,  Dean,  and  Dalrymple  (1985)  recommend  the  values 

k  =  0.15  and  T  =  0.40  for  a  varying  bottom  slope  on  the  basis  of  favorable 
comparison  of  model  predictions  to  laboratory  and  prototype -scale  tank  data. 
These  values  are  used  as  defaults  in  modeling  results  presented  here. 

Ebersole  (1987)  evaluated  the  performance  of  the  breaker  decay  model  using 
field  measurements  and  found  good  results  for  engineering  applications.  He 
also  tested  the  breaker  decay  model  developed  by  Svendsen  (1984)  and  found  it 
to  be  comparable  but  slightly  inferior  in  minimizing  the  rms  error  between 
measurements  and  calculation. 

Dispersion  relation 

83.  The  wave  group  speed  C&  is  related  to  the  phase  speed  of  an  individ¬ 
ual  wave  C  through  the  factor  n,  which  is  a  function  of  the  water  depth  and 
the  wavelength  L  (or  wave  period  T ) : 
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(12) 


Cg  =  nC 

The  phase  or  form  speed  of  a  wave  is  determined  through  the  dispersion 
relation, 

C  =  C0  tanh  (13) 


where 


(14) 


The  quantity  C0  is  the  wave  speed  in  deep  water.  Equation  13  is  one  form  of 
the  dispersion  relation  for  linear  waves;  it  is  an  implicit  equation  that  is 
solved  numerically  since  C  depends  on  the  wavelength,  C  —  L/T. 

Cross -shore  momentum  equation 

84.  As  waves  propagate  over  a  sloping  bottom,  a  variation  in  the  flux  of 
momentum  (variation  in  radiation  stress)  arises  because  of  shoaling,  refrac¬ 
tion,  and  breaking  and  alters  the  mean  water  elevation.  Wave  shoaling  prior 
to  breaking  produces  an  increase  in  wave  height  and  a  corresponding  increase 
in  momentum  flux.  This  flux  increase  is  balanced  by  lowering  of  the  mean 
water  elevation,  called  setdown  (Longuet-Higgins  and  Stewart  1962,  1963). 
Inside  the  surf  zone,  as  waves  break  and  decrease  in  height,  the  momentum  flux 
decreases  and  an  increase  in  mean  water  elevation  occurs,  known  as  setup. 

85.  The  displacement  of  the  mean  water  surface  (setup  or  setdown) 
produced  by  waves  (neglecting  wind  for  present  discussion)  is  determined  from 
the  momentum  equation, 


dS^ 

dx 


-pgd 


dr] 

dx 


(15) 


where  S ^  =  radiation  stress  in  the  direction  of  the  waves  (N/m) .  The 
radiation  stress  component  Sxx  is  given  by  linear-wave  theory  for  an  arbitrary 
wave  angle  of  incidence  as: 
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=  ~pgH 2 


n(c  os20  +  1) 


1 

2 


(16) 


Snell's  Law 

86.  The  energy  flux  conservation  equation  (Equation  6,  or  Equation  7  if 
the  wave-current  interaction  is  included) ,  the  wave  dispersion  relation 
(Equation  13) ,  and  the  cross-shore  momentum  equation  (Equation  15)  are  not 
sufficient  to  specify  the  wave  height  if  the  waves  arrive  at  an  oblique 
(nonzero)  angle  with  respect  to  the  bottom  contours.  In  this  case  an  addi¬ 
tional  equation  expressing  the  conservation  of  wave  number  is  needed  to  obtain 
the  wave  angle.  The  wave  number  conservation  equation  is  written  in  general 
form  as : 

d  I  sin8\  _ 
cbc\  L  /  dy 


87.  For  straight  and  parallel  bottom  contours,  as  assumed  throughout 
this  report,  there  is  no  longshore  variation  in  wave  properties,  and  Equa¬ 
tion  17  simplifies  to: 


_d_ 

dx 


=  0 


(18) 


If  Equation  18  is  integrated,  Snell's  Law  is  obtained,  which  states  that  the 
quantity  sind/L  is  constant.  Equations  6  (or  7),  13,  15,  and  18  determine  the 
wave  height  distribution  after  appropriate  boundary  conditions  are  estab¬ 
lished.  If  the  wave-current  interaction  is  included,  Equation  18  is  modified 
to  describe  additional  refraction  produced  by  the  mean  current.  In  this  case 
Snell's  Law  becomes: 


d  I  C 
dx  \  sinG 


+  V\ 


=  0 


(19) 
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Numerical  Solution  Scheme 


88.  An  explicit  finite-difference  scheme  is  used  to  numerically  solve 
the  governing  equations  determining  the  wave  height  as  a  function  of  distance 
across  the  shore.  The  following  presents  the  solution  scheme  and  the  differ¬ 
ence  equations  employed.  Figure  2  gives  a  definition  sketch  of  the  grid  and 
defines  quantities  entering  the  difference  equations . 


Figure  2.  Definition  sketch  for  the  numerical  grid 

89.  The  numerical  calculation  starts  at  the  most  seaward  point  on  the 
grid  and  proceeds  onshore,  where  quantities  known  at  one  grid  point  are  used 
to  determine  corresponding  quantities  at  the  next  grid  point  closer  to  shore. 
In  the  difference  equations ,  the  index  i  denotes  the  number  of  a  specific  grid 
point.  A  staggered  grid  is  used  with  quantities  defined  at  either  points  in 
the  middle  of  a  calculation  cell  or  at  the  boundaries  between  cells  (Fig¬ 
ure  2).  The  fundamental  quantity  in  the  middle  of  a  cell  is  the  water  depth, 
and  these  grid  points  will  be  referred  to  as  "h-points"  (depth  points).  At 
the  boundaries  of  calculation  cells,  the  main  quantity  is  the  wave  height,  and 
grid  points  located  here  are  known  as  "//-points"  (wave  points) .  All  wave- 
related  quantities  are  defined  at  //-points.  The  numbering  convention  is  such 
that  the  wave  height  Hi  is  defined  on  the  left  boundary  of  calculation  cell 
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number  i.  Finally,  it  is  noted  that  if  the  grid  contains  N  h-points,  it  will 
have  N+l  //-points. 

90.  The  incident  wave  height,  period,  and  direction  must  be  available  at 
the  most  seaward  grid  point  (directly  specified  or  calculated  from  some 
reference  depth  offshore)  prior  to  starting  the  calculation.  From  knowledge 
of  the  wave  properties  at  the  seaward-most  grid  point  (an  //-point) ,  wave 
setdown  at  the  seaward-most  point  is  determined  from  the  analytic  solution  to 
Equation  15  first  presented  by  Longuet-Higgins  and  Stewart  (1962): 


*1 


-71 H2 


AL  sinh 


(20) 


The  mean  water  surface  elevation  is  calculated  at  //-points,  and  depths  at  cell 
boundaries  are  obtained  by  linear  interpolation,  except  for  the  first  and  last 
//-point  were  the  depth  at  the  nearest  cell  is  used.  Similarly,  the  energy 
flux  and  the  radiation  stress  in  the  direction  of  the  waves  are  determined 
from  depth  and  wave  quantities  at  the  first  //-point,  proceeding  shoreward. 

91.  In  going  from  one  //-point  to  the  next,  wave  refraction  is  first 
determined  if  the  incident  wave  approaches  at  a  nonzero  angle  to  the  bottom 
contours.  From  Equation  18,  the  angle  dL  between  wave  crests  and  bottom 
contours  at  the  next  shoreward  //-point  is  given  by: 


0;  =  arc  sin 


inf-^-sin0^. 

V  ^1*1 


(21) 


Note  that  the  grid  point  numbering  decreases  in  the  shoreward  direction  since 
the  x-axis  points  offshore.  The  water  depth  at  grid  point  I  in  Equation  21  is 
corrected  with  the  setdown  or  setup  value  calculated  at  //-point  i+1  (half  a 
cell  seaward)  in  order  to  eliminate  the  need  for  iteration  (Dally  1980). 
Wavelengths  are  calculated  using  the  dispersion  relation  (Equation  13) , 
numerically  solved  by  a  Pade'  approximation  (Hunt  1979).  Wave  angles  in 
Equation  21  are  defined  at  H-points . 
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92.  The  next  step  in  the  calculation  is  to  determine  the  change  in 
energy  flux  and  thus  the  wave  height.  Equation  6  is  written  in  difference 
form  as , 


F  =  Eui  (cos6i+1  -  O.S Aci)  +  AciFsi 
1  cosGj  +  0 . 5 Acl 


(22) 


where 


kAx 

hi  +  *1  i*i 


(23) 


in  which  Ax  is  the  length  step  or  grid  cell  length  (meters) .  The  stable  wave 
energy  flux  is  determined  from: 


Fsi  =  -§  P^[r(h_;  +  tli+1)]2 


( Cgi  +  j) 
2 


(24) 


An  average  value  of  the  group  speed  over  a  cell  is  used  since  this  quantity  is 
defined  at  tf-points,  whereas  the  stable  wave  energy  flux  is  evaluated  at  h- 
points.  Seaward  of  the  break  point,  k  should  be  set  to  zero,  implying  that 
Aci  is  also  zero,  because  energy  dissipation  by  breaking  does  not  occur. 

93.  After  the  energy  flux  has  been  calculated  at  a  specific  point,  the 
corresponding  wave  height  is  determined  from  Equations  3  and  4: 


Hi 


F;  \ 


a  P  3Cgi 


(25) 


Using  this  wave  height,  the  radiation  stress  is  calculated  from  Equation  16, 
and  the  mean  water  surface  displacement  is  given  from  Equation  15  expressed  in 
difference  form: 
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(26) 


Tli  =  Ai+i  + 


P  9(ht  +  T1i+1) 


As  before,  r?  is  evaluated  at  an  ff-point  half  a  cell  seaward  of  the  h-point  in 
calculating  the  depth  to  avoid  use  of  an  implicit  scheme. 

94.  The  alongshore  radiation  stress  component  S ^  is  calculated  at  H- 
points  from  Equation  8  for  use  in  the  longshore  current  model.  Thus,  if  the 
wave-current  interaction  is  neglected  all  wave  quantities  needed  to  calculate 
the  wave  height  can  be  computed  by  a  single  pass  with  an  explicit  solution 
scheme  as  described  above.  However,  if  the  influence  of  a  steady  current 
cannot  be  neglected,  the  wave  and  current  distributions  have  to  be  solved  for 
simultaneously  or  iteratively.  To  include  the  wave-current  interaction,  the 
numerical  representations  of  the  wave  energy  flux  equation  and  Snell's  Law 
must  be  modified.  Equation  21  describing  wave  refraction  is  replaced  by  the 
following  equation  according  to  Equation  19: 


Qi  =  arc  sin 


r  Li 

-  f>i-i  -  V^) 

sin0i+1 

(27) 


Longshore  current  speeds  are  evaluated  at  h-points  (see  Part  V)  and  are 
interpolated  to  obtain  values  at  //-points  for  the  refraction  calculation.  To 
modify  the  difference  form  of  the  energy  flux  equation,  an  extra  term  appears 
in  the  numerator  of  Equation  22,  according  to  Equation  7: 


dV 


dx 


*  [(^>M 


(-V; 


( 

4 


V- 


4  Ax 


(28) 


(The  Ax  factor  in  the  denominator  of  Equation  28  will  be  canceled  by  the 
corresponding  factor  from  the  right  side  of  Equation  7.) 

95.  Because  the  longshore  current  velocity,  initially  unknown  if  it  has 
a  wave -induced  contribution,  occurs  both  in  the  refraction  calculation  and  the 
energy  flux  equation,  independent  solution  of  the  wave  height  distribution  is 
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not  possible.  An  iterative  scheme  must  be  employed.  The  calculation  is 
started  by  assuming  no  current;  thereafter,  the  current  distribution  is 
determined  and  the  wave  computation  updated,  proceeding  in  this  manner  until 
the  wave  and  current  calculations  converge  with  some  specified  accuracy 
(typically,  1  percent)  at  every  grid  point.  The  influence  of  the  wave-current 
interaction  is  small  for  all  but  extremely  high  current  speeds,  as  will  be 
discussed  in  Part  VI,  so  convergence  is  rapid. 

96.  At  every  calculation  step,  a  check  is  made  to  determine  if  depth- 
limited  wave  breaking  has  occurred  according  to  a  criterion  that  specifies  the 
maximum  permissible  ratio  between  wave  height  and  water  depth  yh  (called  the 
breaker  index)  for  an  unbroken  wave.  This  index  typically  varies  between  0.7 
and  1.2,  with  the  values  0.78  and  1.0  most  commonly  used  in  modeling.  A 
dependence  of  yb  on  wave  steepness  and  beach  slope  has  been  reported  by 
several  authors  (e.g.,  Galvin  1969,  Weggel  1972,  Singamsetti  and  Wind  1980, 
Sunamura  1980,  Larson  and  Kraus  1989,  Smith  and  Kraus  1990).  Beach  slope  and 
wave  steepness  may  be  combined  to  yield  the  surf  similarity  parameter  £ 
(Battjes  1974)  used  to  form  a  predictive  equation  for  yb,  as  reported  by 
Sunamura  (1980)  and  Larson  and  Kraus  (1989), 


H, 


Yb=  =1-145 


0.21 


(29) 


where  the  surf -  similarity  parameter  is  defined  as: 

{ -  “■»(£)'* 


(30) 


In  Equations  29  and  30,  the  subscripts  b  and  o  denote  breaking  and  deepwater 
conditions,  respectively,  and  tan/3  is  the  average  slope  of  the  sea  bottom 
along  approximately  one -third  of  the  local  wavelength  seaward  of  the  break 
point . 

97.  Once  breaking  is  initiated,  the  wave  decay  coefficient,  k  is 
assigned  the  value  0.15,  and  energy  dissipation  takes  place.  Also,  if 
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breaking  has  occurred  and  if  the  wave  energy  flux  falls  below  the  stable 
energy  flux,  a  check  is  made  for  wave  reformation.  The  stable  wave  height 
coefficient  is  set  to  T  =  0.40  and,  if  there  is  wave  reformation,  k  is  set  to 
zero.  An  arbitrary  number  of  break  points  with  intermediate  zones  of  wave 
reformation  can  be  represented  by  the  numerical  model . 
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PART  IV:  VERIFICATION  OF  WAVE  MODEL 


98.  The  Dally  (1980)  breaker  decay  model  has  been  extensively  verified 
with  laboratory  and  field  data  and  demonstrated  to  provide  a  good  description 
of  wave  height  decay  in  the  surf  zone  (Dally,  Dean,  and  Dalrymple  1985, 
Ebersole  1987).  Another  positive  feature  of  the  model  is  that  variability  in 
the  two  empirical  coefficients  for  different  wave  and  topographic  conditions 
is  small.  In  this  section,  model  performance  for  complex  bathymetries 
involving  bar  and  trough  formations  is  evaluated  using  data  from  a  large  wave 
tank  experiment  (Kajima  et  al.  1983a,  b) .  Also,  the  capability  of  the  model 
to  reproduce  changes  in  the  statistical  characteristics  of  a  wave  height 
probability  density  function  is  demonstrated.  Additional  wave  model  verifica¬ 
tion  for  laboratory  and  field  data  sets  for  cases  involving  longshore  current 
measurements  is  presented  together  with  verification  of  the  longshore  current 
model  in  Part  VI . 


Large  Wave  Tank  Data  Simulations 

99.  The  accurate  and  controlled  measurements  reported  by  Kajima  et  al . 
(1983a,  b)  are  well  suited  for  verifying  simulation  models  of  waves  transform¬ 
ing  over  complex  bottom  bathymetries.  The  measurements  were  made  in  the  large 
wave  tank  of  the  Central  Research  Institute  for  Electric  Power  Industry 
(CRIEPI)  in  Japan.  The  tank  is  205  m  long,  3.4  m  wide,  and  6  m  deep,  and  wave 
measurements  were  made  using  resistance  wave  gages.  The  primary  objective  of 
the  CRIEPI  program  was  to  investigate  beach  profile  change  under  wave  action. 
Most  tests  were  started  from  a  plane  slope  (1/50  to  1/10) ,  and  profile 
evolution  caused  by  the  regular  breaking  waves  was  measured  through  time. 
Twenty- four  tests  were  made  with  wave  periods  ranging  from  3  to  12  sec  and 
wave  heights  from  0.3  to  1.8  m  in  the  horizontal  portion  of  the  flume.  Tests 
were  performed  separately  with  two  grain  sizes,  0.47  and  0.27  mm  (quartz 
sand) . 

100.  Because  the  initial  beach  profile  was  not  in  equilibrium  with  the 
imposed  waves,  wave  action  redistributed  the  sand  to  approach  a  stable  profile 
over  the  course  of  the  run.  Most  of  the  tests  exhibited  offshore  sand 
movement  and  corresponding  bar  formation  with  the  waves  breaking  on  the  bar . 
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These  tests  thus  constitute  an  excellent  source  for  evaluating  the  performance 
of  the  wave  model  for  barred  profiles,  one  of  the  main  requirements  entering 
model  development  in  the  present  study. 

101.  Larson  and  Kraus  (1989)  used  a  portion  of  the  CRIEPI  data  set  to 
investigate  the  sensitivity  of  the  wave  model  to  the  two  empirical  parameters 
k  and  T  in  Equations  6  and  11.  A  slight  dependence  on  beach  slope  was  noted 
(Dally,  Dean,  and  Dalrymple  1985),  but  a  predictive  relationship  could  not  be 
developed  with  any  level  of  significance.  In  the  present  report,  the  values 
«  =  0.15  and  F  =  0.40,  as  recommended  by  Dally,  Dean,  and  Dalrymple  (1985), 
were  employed  to  test  the  model.  The  reasoning  was  that  if  these  values 
determined  for  plane  beaches  hold  well  for  complex  bottom  profiles  and 
prototype-scale  waves,  then  the  recommended  parameter  values  could  be  confi¬ 
dently  used  as  default  values. 

102.  Wave  height  measurements  made  at  the  end  of  the  runs  were  used  in 
the  present  study  since  the  profile  was  close  to  equilibrium  (changing 
little) ,  and  the  bars  and  troughs  were  most  prominent.  A  large  number  of 
cases  were  simulated  to  test  the  wave  model  in  the  present  study;  here, 
typical  cases  showing  the  predictive  capability  of  the  model  are  described. 

The  numbering  of  the  cases  follows  the  notation  of  Kajima  et  al .  (1983a).  The 
elapsed  times  bracketing  the  wave  height  measurements  and  the  time  of  the 
bottom  profile  measurement  input  to  the  wave  model  are  given  for  each  case. 

The  grid  spacing  was  0.5  m. 

Case  3-1 

103.  This  case  exhibited  bar  development  with  material  composing  the  bar 
removed  from  the  foreshore  and  from  the  area  seaward  of  the  bar.  The  wave 
period  was  relatively  long  (T  =  9.1  sec,  H  =  1.07  m) ,  promoting  development  of 
a  broad,  flat  bar.  Figure  3  shows  the  comparison  between  measured  (67.1  to 
67.6  hr)  and  calculated  wave  height  together  with  the  bottom  profile  measured 
at  71  hr.  The  model  reproduces  the  bipeaked  wave  height  distribution  across 
the  shore.  Equation  29,  developed  with  all  the  CRIEPI  breaking  wave  data 
(Larson  and  Kraus  1989) ,  was  used  to  determine  the  breaker  index  and  gave  a 
good  prediction  of  the  location  of  the  break  point  for  this  case. 
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Figure  3.  Comparison  for  CRIEPI  Case  3-1 


Case  3-4 

104.  Extremely  steep  waves  were  applied  in  this  case  with  a  generated 
wave  height  of  1.54  m  and  a  period  of  3.1  sec.  A  very  peaked  bar  developed  in 
the  offshore  with  much  steeper  flanks  than  are  formed  in  the  field  under 
varying  waves  and  water  level,  and  a  second  small  bar  was  created  near  the 
shoreline.  This  bottom  profile  provides  a  severe  test  of  the  wave  model  for 
which  steep  shoaling  is  required  at  the  outer  bar  followed  by  wave  reformation 
in  the  trough  and  second  breaking  on  the  inner  bar. 

105.  Figure  4  compares  the  measured  (69.6  to  70.1  hr)  and  calculated 
wave  height  distributions  along  the  tank,  displayed  with  the  profile  measured 
after  76.1  hr.  The  measured  wave  height  distribution  is  typical  of  the 
extremely  peaked  bottom  profile  cases.  The  model  reproduces  the  location  of 
the  break  point  and  the  initial,  steep  portion  of  the  wave  height  decay  after 
breaking.  However,  the  calculated  waves  reformed  earlier  in  their  travel 
across  the  trough  landward  of  the  sharp  bar,  the  model  underpredicting  the 
energy  dissipation  and  wave  decay.  An  improved  simulation  can  be  obtained 
without  modifying  model  parameter  values  by  calculating  with  a  smoothed 
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profile  (Larson  and  Kraus  1989)  or  by  requiring  wave  dissipation  to  occur  over 
a  minimum  distance  before  reformation  can  take  place  (introducing  an  addition¬ 
al  empirical  parameter) .  Both  of  these  modifications  can  be  physically 
justified. 
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Figure  4.  Comparison  for  CRIEPI  Case  3-4 


Case  4-4 

106.  The  initial  profile  for  this  case  was  the  final  profile  from  a 
previous  case  involving  steep  waves,  in  which  a  peaked  bar  developed  in  the 
seaward  part  of  the  profile.  Milder  (accretionary)  wave  conditions 
(H  =  1.0  m,  T  =  9.0  sec)  were  employed  in  Case  4-4,  which  flattened  the 
offshore  bar  considerably,  and  the  break  point  was  located  at  a  smaller 
double-bar  feature  that  formed  inshore.  The  result  of  the  model  simulation  is 
illustrated  in  Figure  5,  which  plots  the  measured  wave  height  (8.6  to  9.0  hr), 
the  profile  measured  after  10  hr  of  wave  action,  and  the  calculated  wave 
height  distribution  across  the  tank.  The  overall  model  prediction  is  in  good 
agreement  with  the  measurements  at  and  shoreward  of  the  break  point.  Dif¬ 
ferences  occur,  however,  at  the  two  peaks  in  the  measurements  appearing  at 
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Figure  5.  Comparison  for  CRIEPI  Case  4-4 

about  65  and  95  m.  These  local  maxima  are  believed  to  be  a  result  of  wave 
reflection  on  the  relatively  steep  bottom  slope;  reflection  is  not  described 
in  the  wave  model . 

Case  5-2 

107.  This  case  is  an  example  of  a  beach  with  a  mild  slope  (1/50)  exposed 
to  moderately  steep  waves  ( H  =  0.74  m,  T  =  3.1  sec),  causing  a  small  bar  to 
develop  under  the  breaking  waves.  Figure  6  displays  the  results  of  the  model 
simulation  by  comparing  the  measured  waves  (97.7  to  98.4  hr),  bottom  profile 
after  103  hr,  and  the  calculated  wave  height  across  the  tank.  Agreement  is 
good  overall,  with  small  differences  between  calculation  and  measurements 
appearing  in  the  intermediate  region  where  wave  reformation  took  place.  These 
differences  are  caused  by  the  peaked  bar,  which  induces  rapid  wave  reformation 
in  the  model  as  was  described  for  Case  3-4. 

Summary 

108.  The  wave  model  satisfactorily  reproduced  measured  wave  heights  over 
the  barred  bathymetry  generated  in  a  large  tank.  For  some  cases  with  very 
steep  waves,  reformation  was  predicted  too  early  (too  far  seaward)  in  the  wave 
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Figure  6.  Comparison  for  CRIEPI  Case  5-2 


model,  resulting  in  overprediction  of  the  wave  height  in  zones  of  reformation. 
Premature  reformation  could  be  eliminated  by  slight  smoothing  of  the  profile, 
which  is  physically  justified  because  a  real  wave  probably  responds  over  a 
certain  length  of  the  bottom  and  not  to  abrupt,  highly  localized  changes  in 
bathymetry. 


Random  Wave  Simulations 

109.  Two  different  approaches  may  be  taken  to  simulate  the  transform¬ 
ation  of  random  waves.  Both  approaches  employ  a  deterministic  description  of 
the  wave  transformation  process  itself,  but  the  input  to  the  calculation 
procedure  or  model  is  a  wave  height  probability  density  function  with  known 
statistical  properties  specified  by  different  procedures.  It  is  assumed  here 
that  the  waves  are  narrow  banded  in  frequency  and  direction  so  that  a  single 
frequency  and  direction  are  sufficient  to  characterize  the  sea  state.  Thus, 
randomness  in  the  waves  only  enters  through  variability  in  wave  height. 
Because  the  wave  height  variability  is  described  by  a  probability  density 
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function,  the  waves  are  assumed  to  be  independent  and  have  no  correlation  with 
precedent  or  antecedent  waves.  This  procedure  can  be  extended  to  an  arbitrary 
number  of  wave  direction  and  period  bins  if  the  directional  frequency  spectrum 
is  known. 

Modified  probability  density  approach 

110.  This  method  was  introduced  by  Goda  (1975)  and  further  developed  by 
Battjes  and  Janssen  (1978),  Thornton  and  Guza  (1983),  and  Leont'ev  (1988). 

The  wave  energy  flux  equation  is  integrated  with  respect  to  the  probability 
density  function  for  the  wave  height,  in  which  the  rms  wave  height  is 

used  as  the  single  statistical  descriptor.  If  the  wave  energy  dissipation 
term  is  integrated  to  yield  the  average  dissipation  rate,  an  independent 
procedure  or  criterion  must  be  applied  to  specify  the  percentage  of  broken 
waves  at  each  calculation  point  across  the  shore.  Battjes  and  Janssen  (1978) 
truncated  a  Rayleigh  distribution  above  the  height  of  the  breaking  wave  and 
equated  the  percentage  of  broken  waves  to  the  truncated  probability  mass. 
Thornton  and  Guza  (1983)  assumed  that  the  waves  break  in  proportion  to  the 
distribution  for  all  waves  and  introduced  an  empirical  weighting  function 
which  depends  on  the  ratio  of  wave  height  to  water  depth.  Earlier,  Goda 
(1975)  introduced  a  model  in  which  breaking  occurred  over  a  small  range  of 
water  depths  with  redistribution  of  waves  in  a  truncated  probability  density 
function. 

Monte -Carlo  approach 

111.  The  other  approach  is  to  randomly  select  offshore  waves  according 
to  height  from  a  known  wave  height  distribution  and  transform  each  as  an 
individual  wave  with  a  propagation  and  breaker  decay  model.  A  wave  height 
probability  density  function  is  obtained  at  each  calculation  point  across  the 
shore  from  transformation  of  a  large  number  of  waves  randomly  selected  from 
the  distribution.  The  percentage  of  broken  waves  is  automatically  given  as 
the  number  of  broken  waves  at  the  specific  location  in  relation  to  the  total 
number  of  waves.  Because  each  wave  component  is  calculated  separately,  a 
breaker  criterion  such  as  Equation  29  may  be  used  without  appeal  to  further 
assumptions  as  required  by  the  truncated  probability  density  approach.  The 
disadvantage  of  the  Monte-Carlo  method  is  the  computational  intensiveness, 
because  a  large  number  of  waves  (order  of  100)  must  be  simulated  to  obtain  an 
accurate  representation  of  the  probability  density  function.  A  more  efficient 
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but  less  accurate  way  of  constructing  the  probability  density  function  at 
points  across  the  shore  is  to  discretize  it  with  respect  to  typical  wave 
heights  and  transform  each  such  bin  with  appropriate  weighting  (Dally  1987, 
1990;  Nairn  1988) . 

112.  The  Monte-Carlo  approach  will  be  used  here  instead  of  the  modified 
probability  density  approach  because  limited  information  exists  on  weighting 
functions  suitable  for  describing  the  proportion  of  broken  waves  and  how  they 
vary  across  the  surf  zone,  particularly  for  an  irregular  or  barred  profile. 
Thornton  and  Guza  (1983)  demonstrated  using  field  data  that  the  Rayleigh 
distribution  holds  as  a  good  approximation  primarily  in  deep  water,  and  in 
some  cases  also  surprisingly  well  in  shallow  water.  Ebersole  and  Hughes 
(1987)  and  Hughes  and  Borgman  (1987)  examined  deviations  from  the  Rayleigh 
distribution  of  the  wave  height  in  and  just  outside  the  surf  zone  using  a  data 
set  obtained  with  an  accurate  photogrammetric  technique  to  measure  the  water 
surface  elevation.  Portions  of  the  data  sets  of  Thornton  and  Guza  (1983)  (as 
given  in  Thornton  and  Guza  1989)  and  Ebersole  and  Hughes  (1987)  will  be  used 
here  to  evaluate  the  accuracy  of  the  Monte -Carlo  approach  using  the  Dally 
(1980)  breaker  decay  model  to  predict  characteristic  changes  in  the  statis¬ 
tical  properties  and  distribution  of  wave  height  across  the  profile. 
Verification  of  Monte-Carlo  approach 

113.  It  is  assumed  that  the  Rayleigh  distribution  is  an  adequate 
description  of  the  input  wave  height  distribution  at  the  most  seaward  grid 
point,  presumed  to  be  in  relatively  deep  water  (depth-limited  breaking  does 
not  occur) .  If  p  is  a  uniformly  distributed  random  number  having  a  value 
between  zero  and  one,  the  Rayleigh-distributed  wave  height  corresponding  to 
this  level  of  probability  is  given  by, 


in  which  Hp  =  wave  height  threshold  associated  with  the  exceedance  probability 
of  1  -  p,  and  In  =  natural  logarithm. 

114.  By  numerical  experimentation  using  a  random-number  generator,  it 
was  established  that  approximately  100  wave  selections  are  necessary  to 
reasonably  reproduce  the  Rayleigh  distribution,  and  500  waves  give  good 
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agreement.  Typical  Monte-Carlo  simulations  of  the  cumulative  probability 
distributions  for  100,  200,  and  500  waves  selected  randomly  are  shown  in 
Figure  7.  In  the  figure,  the  simulated  distribution  was  plotted  using  the 
Weibull  plotting  position  formula.  A  sufficiently  large  selection  of  random 
numbers  will  yield  an  accurate  representation  of  a  target  Rayleigh  cumulative 
distribution  once  the  rms  wave  height  is  specified. 

115.  Field  data  on  the  wave  height  transformation  across  the  surf  zone 
collected  by  Ebersole  and  Hughes  (1987)  at  CERC's  FRF  in  Duck,  North  Carolina, 
were  used  to  test  the  Monte-Carlo  simulations.  Ebersole  and  Hughes  filmed, 
with  16 -mm  synchronized  movie  cameras,  the  water  surface  elevation  at  14  poles 
installed  at  nominal  6-m  intervals  across  the  surf  zone  of  a  near-planar 
beach.  This  method,  called  the  "photopole"  method  by  CERC  researchers, 
provides  a  direct  measurement  of  the  wave  height  and  period  and  does  not  rely 
on  a  transfer  function  as  is  necessary  if  the  sea  surface  elevation  is 
inferred  from  another  measured  variable  such  as  pressure  or  wave  orbital 
velocity.  Photopole  data  collection  runs  lasted  approximately  12.5  min  with 
film  exposures  made  every  0.2  sec. 

116.  Ebersole  (1987)  verified  the  applicability  of  the  Dally  breaker 
decay  model  with  field  data.  His  study  concerned  tracking  of  individual  waves 
through  the  surf  zone,  and  he  did  not  apply  the  breaker  decay  model  in  a 
probabilistic  context  to  investigate  the  transformation  of  the  probability 
density  function  or  statistical  wave  parameters  such  as  Hrms.  Hughes  and 
Borgman  (1987)  noted  that  measured  wave  heights  in  the  surf  zone  tend  to  be 
skewed  toward  higher  values  than  predicted  by  the  Rayleigh  distribution.  They 
proposed  a  "Beta-Rayleigh"  distribution  to  describe  the  deviation  from  the 
pure  Rayleigh  distribution  in  shallow  water  in  terms  of  three  empirical 
parameters.  Such  a  distribution  could  easily  be  incorporated  in  the  present 
model  if  further  research  indicates  the  correction  is  warranted. 

117.  To  evaluate  the  Monte -Carlo  approach,  photopole  run  859041510 

(4  September  1985  at  15:10)  from  Ebersole  and  Hughes  (1987)  was  chosen  as  a 
typical  case  for  simulation.  Also,  the  data  provided  the  opportunity  to 
determine  the  degree  to  which  the  probability  density  function  could  be 
predicted  across  the  shore.  As  mentioned  previously,  use  of  the  Monte -Carlo 
technique  requires  no  a  priori  assumptions  about  the  shape  of  the  probability 
density  function  at  locations  along  the  profile.  However,  the  probability 


56 


b .  200  waves 


Figure  7.  Monte-Carlo  simulation  of  Rayleigh 
wave  height  distribution 
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density  function  must  be  specified  at  some  seaward  location  for  input  to  the 
model . 

118.  The  beach  profile  had  a  step- type  shape  with  a  steep  foreshore. 
Quantitative  information  on  the  local  incident  wave  angle  is  not  available, 
but  visual  observation  indicated  that  the  waves  arrived  at  a  small  angle  to 
the  shore  (taken  to  be  zero  in  the  present  work)  as  long-crested  swell.  The 
energy-based  significant  wave  height  measured  at  FRF  Gage  630  located  at  18 -m 
depth  was  0.60  m,  and  the  peak  spectral  wave  period  was  11.3  sec. 

119.  Two  different  approaches  were  taken  in  the  simulations  of  the  field 
data.  The  first  method  used  the  waves  measured  at  the  most  seaward  photopole 
during  the  run  and  transformed  them  individually  across  the  shore.  Sixty- 
three  waves  were  identified  at  this  pole  during  the  12-min  run.  The  other 
approach,  the  true  Monte -Carlo  technique  as  would  be  used  in  a  typical 
engineering  application,  used  the  measured  rms  wave  height  at  the  seaward-most 
photopole  to  define  a  Rayleigh  distribution  from  which  a  large  number  of  waves 
were  randomly  selected  and  transformed  across  the  shore.  The  main  difference 
between  the  two  approaches  is  that  a  number  of  very  small  waves  are  expected 
to  appear  in  the  selection  process  from  the  Rayleigh  distribution;  these  waves 
will  produce  lower  values  of  statistically  determined  wave  heights  at  some 
locations  along  the  profile. 

120.  Parameters  in  the  breaker  decay  model  were  set  to  the  default 
values  of  k  =  0.15  and  T  =  0.40,  and  the  breaker  index  was  set  to  1.0  for  all 
waves  because  the  random  wave  height  simulation  can  produce  unreasonably  high 
wave  steepness  values  in  Equation  29  (breaker  criterion) .  Wave  heights 
measured  at  the  most  seaward  photopole  (Pole  13)  located  in  1.44-m  depth  were 
transformed  seaward  to  3-m  depth  by  linear-wave  theory  for  starting  the 
calculation.  It  is  noted  that  a  few  of  the  waves  at  Pole  13  had  already 
broken,  introducing  a  small  error  in  the  calculations.  The  depths  available 
from  the  photopole  record,  which  are  total  depths,  were  entered  directly  in 
the  governing  calculations,  so  that  no  prediction  of  setup  or  setdown  was 
made . 

121.  Figure  8  shows  a  comparison  between  the  calculated  and  measured  rms 
wave  height  across  the  shore  at  the  photopoles.  The  difference  in  results  of 
the  two  simulation  approaches  is  minor,  with  the  actual  measured  waves  giving 
a  larger  Hims  than  the  random  selection  procedure  seaward  of  the  area  where 
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most  waves  break.  The  maximum  in  the  measured  Hrms  is  located  shoreward  of 
the  calculated  maxima,  but  the  general  trend  in  the  wave  height  decrease  is 
well  reproduced.  Use  of  a  larger  breaker  index  would  give  improved  agreement 
in  the  location  of  the  maxima,  but  this  procedure  is  considered  artificial 
since  the  discrepancy  is  probably  related  to  factors  not  included  in  the 
model,  in  particular,  nonlinear-wave  shoaling. 
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Figure  8.  Comparison  between  measured  and  calculated  Hims 

122.  Output  from  the  wave  model  can  be  used  to  calculate  statistical 
quantities  at  any  point  on  the  calculation  grid,  such  as  the  percentage  of 
broken  waves  and  the  ratio  of  the  rms  wave  height  to  the  average  water  depth. 
The  empirical  probability  density  function  or  normalized  wave  height  histogram 
is  also  readily  constructed  from  the  individual  wave  heights  at  each  point  for 
comparison  to  the  measured  density  function  at  any  photopole. 

123.  Figures  9a  and  b  respectively  give  comparisons  of  measured  and 
calculated  probability  density  functions  at  photopoles  9  and  3;  Pole  9  was 
located  just  shoreward  of  the  area  where  the  higher  waves  began  breaking,  and 
Pole  3  was  located  close  to  shore  where  almost  all  waves  had  broken.  Because 
most  waves  were  calculated  to  break  farther  from  shore,  the  calculated  wave 
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Figure  9.  Calculated  and  measured  wave  height 
probability  density  functions 
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heights  at  Pole  9  are  smaller  than  the  measured  heights ,  and  the  maxima  in  the 
calculated  probability  density  functions  are  displaced  toward  the  lower  waves. 
At  Pole  3,  located  close  to  the  shoreline,  the  maxima  in  the  measured  and 
calculated  probability  functions  coincide,  but  the  shape  of  the  calculated 
function  is  more  peaked.  Figures  9a  and  b  show  that  the  calculations  well 
reproduce  the  overall  shape  and  range  of  the  measured  probability  density 
function  at  the  two  very  different  regions  of  breaking  waves  in  the  surf  zone. 
Summary 

124.  The  Monte-Carlo  technique  of  transforming  individual  wave  com¬ 
ponents  was  successful  in  reproducing  the  rms  wave  height  measured  in  the 
field  on  an  irregular  bottom  topography.  The  main  deviation  between  calcula¬ 
tion  and  measurements  appeared  at  the  location  of  the  maximum  height,  for 
which  the  calculated  waves  tended  to  break  seaward  of  the  actual  region  of 
significant  wave  breaking.  The  two  calculation  approaches  of  propagating  in¬ 
dividual  waves  from  the  3-m  depth  and  randomly  selecting  waves  from  a  Rayleigh 
distribution  defined  by  the  rms  wave  height  gave  almost  the  same  results .  It 
is  concluded  that  the  Raleigh  wave  height  distribution  is  an  adequate  descrip¬ 
tor  of  waves  entering  the  surf  zone,  and  specification  of  the  rms  wave  height 
or  similar  statistical  wave  height  at  a  location  seaward  of  the  first  breakers 
is  sufficient  for  accurately  estimating  random  wave  transformation  over  a  surf 
zone  with  irregular  bottom  topography. 
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PART  V:  LONGSHORE  CURRENT  MODEL 


125.  This  chapter  presents  the  numerical  model  of  the  longshore  current 
developed  in  this  study.  The  model  calculates  the  longshore  current  across  a 
barred  profile  under  forcing  by  arbitrary  combinations  of  waves  arriving  at  an 
oblique  angle  to  shore  and  a  steady  wind  blowing  over  the  sea  surface  at  an 
arbitrary  angle.  The  longshore  components  of  the  wave  and  wind  forces 
generate  a  longshore  current,  and  the  shore-normal  components  change  the  mean 
water  surface  elevation,  resulting  in  setup  and  setdown.  The  main  restriction 
of  the  model  is  that  the  beach  contours  must  be  uniform  alongshore. 

126.  First,  the  theoretical  background  is  given,  and  the  governing  equa¬ 
tions  of  the  model  are  derived  starting  from  the  time-averaged  and  depth- 
integrated  equations  of  motion.  The  formulation  of  the  bottom  friction  term 
in  the  alongshore  momentum  equation  is  examined  with  various  degrees  of 
approximation  to  determine  the  accuracy  of  simplifications  made  to  increase 
computational  efficiency.  Two  mathematical  descriptions  of  the  bottom 
friction  are  presented,  resulting  in  a  linear  and  a  nonlinear  longshore 
current  model.  Representation  of  lateral  mixing  is  also  discussed,  and  an 
eddy  viscosity  coefficient  is  proposed  for  calculating  currents  generated  over 
a  barred  profile.  Finally,  the  discretization  procedure  for  the  governing 
equations  and  the  numerical  solution  scheme  are  described,  leading  to  the 
computationally  efficient  double -sweep  implicit  technique  for  obtaining  the 
velocity.  The  current  model  is  verified  in  Part  VI. 


Governing  Equations  for  Nearshore  Currents 

127.  The  vertically  integrated,  time-averaged  momentum  equations  may  be 
written  (Mei  1983) , 


VfX  4  Vfy  -  ~<>fX  4  ^ 


(32) 


+  v^V  =  _  aq  + 


U  dx  +  V  dy 


3  0y  +  ^hy  +  Ly  +  Rby  +  RSy 


(33) 
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where 


U  =  time-averaged  (mean)  cross-shore  current  (m/sec) 

V  =  time-averaged  (mean)  longshore  current  (m/sec) 

Fbx  =  cross -shore  component  of  bottom  friction  term  (m/sec2) 

Fby  —  longshore  component  of  bottom  friction  term  (m/sec2) 

Lx  =  cross -shore  component  of  lateral  mixing  term  (m/sec2) 

Ly  =  longshore  component  of  lateral  mixing  term  (m/sec2) 

Rbx  =  cross -shore  component  of  the  wave  driving  term  (m/sec2) 

Rby  =  longshore  component  of  the  wave  driving  term  (m/sec2) 

Rsx  =  cross -shore  component  of  the  wind  (surface)  driving  term  (m/sec2) 
Rsy  =  longshore  component  of  the  wind  (surface)  driving  term  (m/sec2) 


Equation  32  is  a  more  general  form  of  the  momentum  equation  in  the  x- direction 
than  Equation  15  used  in  the  wave  model  in  Part  III  to  compute  setup  and 
setdown.  After  a  number  of  simplifying  assumptions,  Equation  15  is  obtained 
from  Equation  32  as  described  below.  In  the  wave-related  terms,  the  subscript 
"5"  is  used  to  denote  the  force  produced  by  shoaling  and  breaking  waves,  since 
the  more  natural  notation  "w"  might  cause  confusion  between  waves  and  wind. 
Wave -driving  terms 

128.  The  wave-driving  terms  are  a  function  of  the  change  in  radiation 
stresses  that  takes  place  as  waves  propagate  toward  shore  and  transform  by 
shoaling,  refraction,  and  breaking.  They  are  expressed  as, 


R 


hx 


=  -1  l  dS*x 

pd^  dx 


dSxy 

dy 


R 


by 


-1  (  dS*y 
p d\  dx 


dy 


(34) 


(35) 


where  the  radiation  stress  components  Sxx  and  Sx y  are  given  by  Equations  16  and 
8,  respectively.  By  linear-wave  theory,  the  radiation  stress  component 
alongshore  Syy  is : 


S^=  -±p  g# 


n(sinz0  +  l) 


(36) 
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Bottom  friction  terms 


129.  Bottom  friction  is  expressed  as  a  quadratic  stress  law  in  the  total 
local  fluid  velocity  and  is  produced  by  the  combination  of  steady  motion  from 
the  mean  wave-  and  wind- induced  current  and  unsteady  motion  from  oscillatory 
waves.  The  instantaneous  shear  stress  components  at  the  bottom  are  written, 

Zbx  ~  cfp^u2  +  Vs  u  (37) 

xby  =  cfP'lu2  +  V^V  (38) 


where  Ct  is  an  empirical  bottom  friction  coefficient,  and  the  x-component  u 
and  y-component  v  of  the  total  current  velocity  (evaluated  at  the  bottom)  are 
given  by: 


u  =  U  +  ub  (39) 

v  =  V  +  vb  (40) 


The  horizontal  velocity  components  ub  and  vb  of  oscillatory  wave  motion  at  the 
bottom  are  given  from  linear-wave  theory  as, 

ub  =  ^cospt^]cos6  (41> 

vb  =  ^cos(^p)  sin0  (42> 


where  the  amplitude  of  the  horizontal  component  of  the  wave  orbital  velocity 
is 


u.  = 


gHT 


2Lcosh 


2nd' 

L 


(43) 
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130.  From  Equations  37  and  38,  the  bottom  friction  terms  in  the  x-  and 
y-  momentum  equations  (Equations  32  and  33)  may  be  obtained,  noting  that  Fbx 
and  Fby  are  not  identical  to  the  shear  stress  but  are  by  definition  the  shear 
stress  components  divided  by  water  density  and  total  depth: 


Fbx  =  ~~Uu2  +  v2  u  > 


(44) 


Fby  = 


(45) 


In  these  equations  and  elsewhere  in  this  report,  triangular  brackets  denote  a 
time -averaging  operation  over  the  interval  of  a  wave  period. 

131.  As  seen  from  the  equations  for  the  oscillatory  velocity  components 
(Equations  41  and  42),  u  and  v  are  time  dependent,  and  a  time  average  will 
produce  a  nonzero  value  because  of  the  existence  of  the  steady  current.  Thus, 
in  solving  the  governing  equations  to  obtain  the  current  velocity,  intensive 
numerical  computations  are  required  not  only  because  of  the  nonlinear  (quad¬ 
ratic)  stress  law,  but  also  to  perform  the  time -averaging  operation. 

Lateral  mixing  terms 

132.  The  lateral  mixing  terms  describe  the  diffusion  of  momentum  and  are 
related  to  the  turbulent  Reynolds  stresses.  Lateral  mixing,  which  is  not  well 
understood  in  the  nearshore,  is  modeled  using  an  eddy  viscosity  approach.  The 
mixing  terms  may  be  expressed  in  general  form  as 


1 

\*l 
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(46) 


(47) 


where  the  £irj  (i,j  =  x,y)  are  the  components  of  an  eddy  viscosity  tensor. 

133.  Knowledge  of  momentum  diffusion  in  the  surf  zone  by  turbulence  is 
limited,  but  observations  of  the  growth  of  a  point  source  of  dye  injected  in 
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the  surf  zone  show  that  its  dimensions  increase  much  more  rapidly  across -shore 
as  compared  with  alongshore,  indicating  that  the  magnitudes  of  the  viscosity 
components  involving  the  cross -shore  component  are  much  greater  than  that  of 
Syy  (Kraus,  Mimura,  and  Horikawa  1980).  In  the  present  case  involving 
homogeneity  alongshore  and  no  cross -shore  current  in  the  steady- state  case 
(U  =  0)  ,  the  only  eddy  viscosity  component  of  interest  is  eyx,  abbreviated  as 
e  hereafter. 

134.  Many  forms  of  the  eddy  viscosity  coefficient  e  have  been  proposed, 
as  summarized  by  Kraus,  Mimura,  and  Horikawa  (1980);  McDougal  and  Hudspeth 
(1986);  and  O'Conner  and  Yoo  (1987).  Traditionally,  a  Prandtl  approach  is 
taken  in  which  the  eddy  viscosity  coefficient  is  assumed  proportional  to 
appropriate  velocity  and  length  scales.  Longuet-Higgins  (1970b)  set  e  ~ 
x(gh)1/2,  which  is  an  intuitively  reasonable  representation  in  the  surf  zone 
and  on  a  uniformly  sloping  beach,  since  the  eddy  size  should  increase  with 
distance  offshore  (or,  equivalently,  with  the  wave  height,  which  is  depth- 
limited  and  proportional  to  distance  offshore) ,  and  the  velocity  scale  should 
be  related  to  um,  which  is  proportional  to  ( gh)1,z  in  shallow  water.  Madsen, 
Ostendorf,  and  Reyman  (1978)  introduced  the  form  e  ~  xum  to  give  the  eddy 
viscosity  coefficient  a  decreasing  value  seaward  of  breaking  waves;  otherwise, 
for  the  surf  zone  their  form  of  e  is  similar  to  that  of  Longuet-Higgins. 

135.  Kraus,  Mimura,  and  Horikawa  (1980)  and  McDougal  and  Hudspeth  (1986) 
developed  generalized  analytic  solutions  of  the  longshore  current  profile  as  a 
function  of  the  form  of  the  eddy  viscosity  coefficient.  Both  groups  concluded 
that  the  current  profile  on  a  plane  beach  was  insensitive  to  the  form  of  the 
coefficient.  On  a  barred  profile,  parameterization  of  e  with  the  cross -shore 
distance  x  is  not  logical  because  wave  breaking  in  the  vicinity  of  bars 
located  at  different  distances  from  the  shore  will  create  local  areas  of 
agitated  white  water  (strong  turbulence) ,  with  possible  quiescent  green  water 
(weak  turbulence)  in  between  bars  where  wave  reformation  can  occur.  Based  on 
their  measurements  of  the  longshore  current  profile  on  a  barred  beach  in  the 
Great  Lakes,  Greenwood  and  Sherman  (1986)  concluded  that  lateral  mixing  was 
enhanced  by  current  velocity  gradients  in  the  vicinity  of  bars. 

136.  In  this  report,  the  local  bottom  orbital  velocity  is  taken  as  the 
characteristic  velocity  and  the  local  wave  height  as  the  typical  length  scale 
to  give , 
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(48) 


e  =  A  umH 


where  A  is  an  empirical  coefficient.  In  this  formulation  the  eddy  viscosity 
has  two  rational  properties : 

a.  It  decays  with  distance  seaward  of  the  surf  zone  as  the  mean 
depth  increases. 

b.  It  is  sensitive  to  local  wave  conditions  as  exist  in  breaking 
over  bars  and  reformation  over  troughs,  in  that  larger  waves 
will  produce  greater  mixing. 

Wind- driving  terms 

137.  The  local  wind  is  another  important  driving  force  in  the  nearshore. 
Birkemeier  and  Dalrymple  (1975)  were  among  the  first  to  consider  the  wind  in 
nearshore  circulation  numerical  modeling.  If  a  wind  blows  over  a  water 
surface,  a  current  will  be  generated,  and  the  surface  will  tilt  in  adjustment 
of  the  water  body  to  the  transfer  of  momentum  at  the  air- sea  interface.  Wind¬ 
generated  currents  have  been  detected  at  depths  on  the  order  of  100  m  (Long 
and  Hubertz  1988),  and  strong  wind  can  produce  significant  setup  and  setdown 
at  the  shoreline  during  storms,  depending  on  the  direction  of  its  cross -shore 
component  (e.g.,  Bretschneider  1966). 

138.  The  wind  stress  appears  as  a  forcing  term  in  the  momentum  equations 
and  is  expressed  by  a  quadratic  drag  law.  For  a  wind  blowing  over  a  water 
surface,  the  shear  stress  rs  at  the  air-water  interface  is  written, 

^  =  CDpaW2  (49) 


where  the  subscript  s  refers  to  forcing  by  wind  at  the  water  surface.  The 
shear  rs  has  the  units  of  N/m2,  and 
CD  =  drag  coefficient 
pa  =  density  of  air  (kg/m3) 

W  =  wind  speed  (m/sec) 

If  the  wind  blows  at  an  angle  to  the  shoreline,  the  wind- induced  shear  stress 
in  Equation  49  is  expressed  as  a  cross -shore  component  rsx  and  a  longshore 
component  rsy, 
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=  CDpa\w\Wcos<p 

(50) 

=  CDpa\w\W3inip 

(51) 

where  <j>  is  the  incident  angle  of  the  wind.  The  angle  definition  employed  in 
Equations  50  and  51  is  shown  in  Figure  1. 

139.  In  the  general  case  of  a  wind  blowing  at  a  nonzero  angle  to  the 
shoreline,  the  driving  terms  describing  the  wind  in  the  momentum  equations, 
Equations  32  and  33,  are: 
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Cp  Pa 
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I  Vf  \  Wcos  <p 
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(53) 


In  a  strict  sense,  in  the  wind-driving  terms  the  velocity  difference  between 
the  wind  and  the  water  surface  W  -  V  should  be  used  instead  of  W  alone.  The 
speed  of  a  coastal  current  induced  by  a  persistent  wind  is  about  3  percent  of 
the  mean  wind  speed  (e.g.,  Hsu  (1988)  and  the  discussion  in  the  next  chapter). 
Because  the  condition  W  »  V  is  expected  to  hold  in  all  applications,  V  can  be 
neglected  in  comparison  to  IV  .  Inclusion  of  V  in  the  wind-driving  term  would 
add  another  nonlinearity  and  complicate  the  numerical  solution,  although  the 
iterative  solution  scheme  used  in  the  nonlinear  model  could  accommodate  such 
an  expression,  but  at  greater  computational  overhead.  From  discussion  in 
Part  III,  the  uncertainty  in  the  value  of  the  drag  coefficient  appears  to 
obviate  further  any  practical  need  to  account  for  the  relative  speed  of  the 
wind  and  current  in  the  wind- driving  terms. 

140.  If  a  steady  wind  blows  over  the  open-ocean  surface,  it  is  well 
known  (Pond  and  Pickard  1977,  Bishop  1984)  that  the  resultant  equilibrium 
current  will  not  be  in  the  direction  of  the  wind,  but  an  angle  to  the  right 
(left)  in  the  Northern  (Southern)  Hemisphere  because  of  the  Coriolis  force. 

The  Coriolis  force  is  not  included  in  the  equations  of  motion  used  in  this 
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report  because  of  the  relatively  narrow  seaward  extent  of  the  longshore 
current,  but  the  possible  necessity  of  including  the  Coriolis  force  in  the 
offshore  boundary  condition  is  noted. 

141.  At  some  sites,  the  coastal  wind  boundary  layer  (Hsu  1988)  should  be 
considered  to  account  for  changes  in  the  local  wind  speed  and  direction. 
Further  research  on  this  subject  is  planned. 

Governing  Equations  for  Longshore  Current  and  Setup 

142.  The  full  governing  equations,  Equations  32  and  33,  are  reduced 
based  on  the  following  assumptions: 

a.  Linear-wave  theory  is  applicable  everywhere,  both  inside  and 
outside  the  surf  zone. 

b.  The  time-averaged  cross-shore  current  U  is  zero. 

c.  The  bottom  contours  are  straight  and  parallel,  indicating 
uniformity  in  the  y-direction  ( d/dy  =0). 

d.  Bottom  friction  in  the  cross-shore  direction  is  small  in 
comparison  with  dS^/dx .  (Retaining  the  cross -shore  friction 
term  couples  the  momentum  equations.) 

These  assumptions  imply  that  all  convective  acceleration  terms  in  Equations  32 
and  33  are  zero.  The  friction  and  lateral  mixing  terms  vanish  from  the 
momentum  equation  in  the  x-direction  (Equation  32),  leaving: 


dS 

XX 

dx 


-  Wcos<p 


(54) 


Neglecting  wind,  this  equation  reduces  to  Equation  15  used  in  the  wave  model 
to  compute  setup  and  setdown.  Because  the  variables  in  Equation  54  depend 
only  on  the  x-coordinate ,  the  partial  derivatives  in  Equation  32  were  changed 
to  full  derivatives. 

After  simplifications  based  on  the  above -listed  assumptions,  the 
remaining  terms  in  the  y-  or  longshore  component  momentum  equation  are  wave 
driving,  wind  driving,  bottom  friction,  and  lateral  mixing.  Equation  33  may 
thus  be  written,  using  Equation  35  for  the  wave-driving  term  and  Equation  47 
for  the  lateral  mixing  term, 
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(55) 
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dx 


-  f, 
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1  dSxy 
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'D-^IMWsin<P 


where  fhy  is  defined  as: 


fby  =  ~  Fbyd  (56) 

Equation  55  is  a  second- order  nonlinear  ordinary  differential  equation  in  the 
current  velocity.  The  assumptions  decouple  the  momentum  equations  (if  the 
wave-current  interaction  is  neglected),  eliminating  the  need  for  simultaneous 
solution  of  the  waves  and  longshore  current,  significantly  reducing  computa¬ 
tional  time  and  improving  stability.  The  nonlinearity  enters  through  the 
bottom  friction  term  containing  a  dependence  on  current  velocity  squared.  The 
next  section  discusses  different  simplifications  and  procedures  for  evaluating 
the  bottom  friction  terms,  and  their  validity. 


Evaluation  of  Bottom  Friction 


144.  The  computational  difficulty  in  evaluating  the  bottom  friction  term 
lies  not  only  in  its  dependence  on  the  current  velocity  squared,  introducing  a 
strong  nonlinearity,  but  also  in  the  time-average  that  has  to  be  calculated  at 
each  model  grid  point  across  the  shore  because  of  the  oscillatory  motion  of 
the  waves  (Equation  45).  From  Equations  41  and  42  for  the  orbital  velocities, 
the  y-component  of  the  friction  term  may  be  expressed  as: 
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145.  Evaluation  of  fby  at  each  model  grid  point  through  the  surf  zone  is 
computation  intensive;  a  16-point  Simpson  integration  or  an  8-  or  16-point 
Gauss  integration  is  typically  used  to  perform  the  time  average.  In  analytic 
treatments  and  in  many  numerical  circulation  models,  simplifying  approaches 
have  been  taken  to  linearize  Equation  57  by  putting  restrictions  on  the  ratio 
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of  the  magnitude  of  the  current  velocity  to  the  wave  orbital  velocity,  as  well 
as  restricting  the  incident  wave  angle. 

146.  As  discussed  in  the  literature  review  in  Part  II,  reliable  theore¬ 
tical  formulas  for  the  friction  factor  cf  are  lacking  for  the  surf  zone, 
although  the  value  in  the  field  seems  to  be  fairly  well  established  in  the 
range  of  approximately  0.005  to  0.01.  Numerical  experimentation  with  NMLONG 
did  not  produce  improved  agreement  between  measurements  and  calculations  using 
the  wave - dependent  friction  factor  of  Jonsson  (1966),  although  the  Jonsson 
factor  did  provide  values  of  the  correct  order  of  magnitude.  To  reduce  model 
calculation  time,  a  constant  value  of  cf  is  used  that  can  be  a  default  value 
or  a  value  determined  in  calibration. 

Weak-current  (linear)  approximation 

147.  The  most  common  approximation  to  Equation  57  (Bowen  1969,  Longuet- 
Higgins  1970a)  results  from  the  assumptions  that  the  current  is  weak  with 
respect  to  the  magnitude  of  the  wave  orbital  velocity  and  that  the  incident 
wave  angle  is  small.  With  these  restrictions  Equation  57  can  be  expanded  in  a 
Taylor  series  and  time-averaged  term  by  term;  after  higher  order  terms  are 
dropped,  a  linear  form  of  fby  in  V  is  obtained: 


fby  =  ^CfUmV 


(58) 


A  higher  order  approximation  was  presented  independently  by  Liu  and  Dalrymple 
(1978)  and  Kraus  and  Sasaki  (1979),  still  assuming  a  weak  current,  but 
allowing  large  incident  wave  angles.  The  time-averaged  friction  term  becomes: 

fby  =  ^cfumV{  1  +  sin28)  (59) 


Strong-current  approximation 

148.  The  case  of  a  strong  current  and  small  waves  is  simple  and  reduces 
Equation  57  to  a  quadratic  velocity  dependence  on  the  current  according  to: 
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fby  =  Cf\V\V 


(60) 


Square -wave  approximation 

149.  A  more  applicable  approach  for  numerical  modeling  was  introduced  by 
Nishimura  (1982,  1988),  who  approximated  the  sinusoidal  time  variation  in  wave 
orbital  velocity  by  a  square  wave  with  the  same  area.  This  construct  allows 
an  explicit  solution  of  Equation  57  to  be  obtained  with  respect  to  the  time¬ 
averaging  procedure,  the  details  of  which  are  given  in  Appendix  A.  For 
negligible  cross-shore  current  (17=0),  Nishimura's  expression  for  the  time- 
averaged  friction  term  is, 

fby  =  cf[z  +  -^sin2ejy  (61) 


where  the  auxiliary  quantities  Z  and  w  are  defined  as: 


Z  =  +  w2  +  2 wVsinB  +  yjv2  +  w2  -  2 wVsinB) 


(62) 


lu- 


(63) 


150.  According  to  Nishimura  (1988)  the  maximum  relative  error  in  this 
approximation  is  10  percent  as  compared  with  high-order  numerical  integration 
of  Equation  57.  (The  accuracy  of  Equation  61  is  investigated  in  Appendix  A.) 
Equation  61  significantly  reduces  the  computation  time  for  the  longshore 
current  because  the  time  averaging  was  performed  analytically;  also,  the 
necessity  of  invoking  the  assumptions  of  a  weak  longshore  current  and  small 
wave  angles  is  eliminated.  A  quadratic  dependence  on  the  velocity  in  Equa¬ 
tion  55  still  remains,  requiring  an  iterative  approach  to  solve  for  the 
current . 

151.  In  the  development  of  NMLONG ,  both  a  linear  (Equation  59)  and  a 
nonlinear  (Equation  61)  bottom  friction  law  were  considered,  yielding  two 
different  numerical  solution  schemes  or  models.  To  investigate  the  validity 
of  the  different  approximations  of  the  bottom  friction  term,  Equation  57  was 
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integrated  numerically  using  20-point  Gauss  quadrature,  yielding  a  quantity 
(-^by^on  accurate  to  about  the  fourth  decimal  place.  The  bottom  friction  term 
given  by  Equation  57  was  normalized  with  the  linearized  bottom  friction 
(-^byhin  according  to  Equation  58.  The  ratio  (fby)non/(^by)iin  depends  on  two 
nondimens ional  variables  as  shown  in  Figure  10,  the  ratio  of  magnitudes  of  the 
current  and  the  wave  orbital  velocity  V/um,  and  the  incident  wave  angle  6. 


fnon  /  f  lin 


Figure  10.  Comparison  of  linear  and  quadratic  bottom  friction 


152.  Figure  10  displays  the  deviation  between  the  numerically  evaluated 
(complete)  friction  term  and  the  linearized  term  as  a  function  of  V/um  and  6. 
The  linear  friction  formulation  gives  errors  as  great  as  50  percent  for  strong 
currents  and  large  wave  angles.  The  influence  of  the  angle  is  more  effective 
at  small  values  of  V/um,  because  the  term  V 2  dominates  for  strong  currents. 

If  the  magnitude  of  the  current  velocity  equals  that  of  the  wave  orbital 
velocity,  the  bottom  friction  stress  for  the  linearized  formulation  is  only 
half  the  correct  value  if  the  same  value  of  cf  is  applied  in  both  cases.  This 
discrepancy  can  be  avoided  in  the  pragmatic  sense  by  calibrating  the  linear 
current  model  with  a  larger  value  of  cf.  Thus  linearized  models  have  been 
applied  successfully  at  F/um-ratios  well  beyond  the  strict  limit  of  mathemati¬ 
cal  validity  of  the  weak- current  approximation. 
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Approximation  to  the  Nonlinear  Solution 


153.  If  lateral  mixing  is  omitted  it  is  possible  to  derive  a  correction 
factor  to  be  applied  to  the  solution  of  the  linearized  model  to  obtain  the 
current  velocity  from  the  nonlinear  model.  This  factor,  called  a,  is  a 
function  of  V/um  and  8,  and  must  for  the  general  case  be  represented  in 
tabular  or  graphical  form.  Omitting  mixing  in  Equation  55,  the  driving  force 
is  balanced  by  bottom  friction.  Therefore,  for  a  certain  driving  force  (waves 
and/or  wind) ,  the  bottom  friction  force  will  be  identical  whether  a  linear  or 
a  quadratic  friction  law  is  used.  By  equating  the  friction  formulations 
specified  by  Equations  57  and  58,  a  relation  is  derived  between  the  velocity 
predicted  by  the  nonlinear  model,  numerically  evaluated  as  Vnon ,  and  that  pre¬ 
dicted  by  the  linear  model  Wlin.  An  equation  for  a  is  obtained, 


sjcoanz  (cositz  +  2at|rO  +  a24r2 


+ 


£  CQSTt  Z' 

onir  , 


(64) 


where  z  =  2 t/T  is  a  dummy  integration  variable  for  the  transformation  made  to 
obtain  convenient  limits  for  Gauss  integration,  and 

a  =  (65) 

vlin 


*  = 


' lin 


(66) 


C  =  sin8 


(67) 


154.  A  Newton -Raphs on  technique  was  employed  to  find  the  solution  of 
Equation  64.  In  Figure  11  the  factor  a  is  displayed  as  a  function  of  Rlin/um 
and  8.  As  also  seen  in  Figure  10,  the  incident  wave  angle  is  only  important 
for  small  values  of  Flin/um,  because  the  solution  approaches  the  strong-current 
case  for  large  values  on  Vlin/um. 

155.  As  previously  stated,  a  may  be  used  to  correct  the  solution 
obtained  with  a  linearized  friction  term  in  the  absence  of  lateral  mixing. 
Without  mixing,  the  current  velocity  at  any  point  across  the  shore  is  indepen¬ 
dent  of  the  velocity  at  neighboring  points  and  may  be  solved  for  independently 
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Figure  11.  Correction  factor  for  nonlinear  friction  term 

once  local  wave  and  wind  properties  are  known.  The  factor  a  is  determined  by- 
local  wave  properties  that  enter  through  the  wave  orbital  velocity  and  wave 
angle . 

156.  The  factor  a  may  be  interpreted  in  another  way,  namely,  as  a 
function  of  the  bottom  friction  coefficient  cf.  Identical  current  velocities 
may  be  obtained  from  the  nonlinear  and  linear  models  if  different  values  of 
the  friction  coefficient  (cf)non  and  (cf)lin  are  employed  in  the  respective 
models  (still  neglecting  lateral  mixing)  such  that  (cf)non  =  a  (cf)Un,  with  a 
defined  as  before.  It  is  noted  that  even  if  (cf)lin  is  taken  as  a  constant, 
the  apparent  (Cf)non  will  vary  in  accordance  with  the  local  wave  properties 
because  a  is  a  function  of  those  properties. 

157.  Limiting  cases  of  Equation  64  provide  insights  to  previous  work. 

If  V  approaches  zero  (weak-current  assumption),  the  limiting  value  is: 
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(68) 


If  dimensional  quantities  are  reintroduced,  the  solution  for  the  current  will 
be  (dropping  the  index  "non"  for  the  improved  solution  of  the  current) : 


V  = 


Rby  +  Rsy 


cfum( 1  +  sin20) 


(69) 


Equation  69  is  identical  to  the  solution  obtained  with  the  large -angle  and 
weak-current  friction  formulation  used  by  Liu  and  Dalrymple  (1978)  and  Kraus 
and  Sasaki  (1979) ,  omitting  the  wind-driving  term.  This  expression  for  the 
bottom  friction  involves  no  restrictions  in  terms  of  the  wave  angle,  but  only 
depends  on  the  assumption  of  a  weak  current  with  respect  to  the  magnitude  of 
the  bottom  wave  orbital  velocity  and  is,  in  that  sense,  exact. 

158.  The  other  limiting  case  of  Equation  64  to  consider  is  ->  «>(  the 
so-called  strong- current  approximation.  In  this  case,  Equation  64  may  be 
solved  to  yield: 

“  ■  (^) 


Expressed  in  dimensional  quantities,  Equation  70  simplifies  to: 

Rby  +  jggyjj 


(71) 


This  solution  corresponds  to  the  formulation  given  in  Equation  60. 

159.  Lateral  mixing  is  expected  to  enter  in  applications,  which  elimin¬ 
ates  the  possibility  of  using  a  simple  approach  with  a  correction  factor  to 
obtain  a  rapid  solution  of  the  nonlinear  current  equation.  However,  compu¬ 
tation  time  is  reduced  by  applying  the  correction  factor  to  improve  the 
initial  estimate  in  an  iterative  scheme  involving  the  quadratic  bottom  stress 
term. 
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Evaluation  of  the  Wind  Drag  Coefficient 


160.  Two  problems  are  associated  with  applying  Equation  49  for  the  wind 
stress,  determining  the  value  of  the  drag  coefficient  CD  and  choosing  a  proper 
wind  speed .  These  problems  are  related  because  Cq  is  expected  to  vary  with 
the  reference  height  of  the  wind  according  to  the  variation  of  the  wind  speed 
with  elevation.  In  most  applications  the  drag  coefficient  is  referenced  to 
the  wind  speed  given  at  a  level  of  10  m.  The  drag  coefficient  is  a  bulk 
description  of  the  flow  resistance  offered  by  the  water  surface  and  is 
composed  of  two  parts,  frictional  drag  and  form  resistance.  Frictional  drag 
is  caused  by  the  shear  along  the  water  surface,  and  form  drag  is  produced  by 
pressure  differences  associated  with  the  shape  of  the  surface.  Thus  CD  is 
expected  to  be  a  function  of  several  parameters  including  a  representative 
roughness  length  scale  of  the  surface  and  the  general  air  flow  conditions. 

A  difficulty  in  determining  the  value  of  the  drag  coefficient  at  an  air- sea 
interface  is  the  mobility  of  the  water  surface  and  the  ease  in  which  it 
deforms.  Strong  winds  create  waves  and  thus  increase  the  surface  roughness  as 
compared  with  more  moderate  winds. 

161.  A  large  number  of  studies  have  been  performed  with  the  objective  of 
determining  the  drag  coefficient  over  a  water  surface  (or  roughness  length 
scale  of  the  surface) .  Garratt  (1977)  and  Long  and  Hubertz  (1988)  review 
such  studies  and  summarize  measured  values  of  CD.  Typical  values  are  in  the 
range  of  1  10  3  to  3  10  3  with  the  associated  wind  speed  in  the  interval  5  to 
20  m/sec,  although  large  scatter  in  the  data  is  common.  Long  and  Hubertz 
(1988)  reported  a  factor  of  32  between  the  lowest  (0.2  10~3)  and  highest 

(6.4  10  3)  value  of  CD  found  in  the  literature.  A  dependence  of  CD  on  the  wind 
speed  is  also  found  (Garratt  1977,  Large  and  Pond  1981,  Amorocho  and  DeVries 
1980,  The  WAMDI  Group  1988),  at  least  for  higher  wind  speeds,  and  some  recent 
studies  also  attempt  to  quantify  the  influence  of  waves  (Huang  et  al .  1986). 

162.  The  classical  method  for  determining  the  drag  coefficient  is  based 
on  a  logarithmic  velocity  profile  above  the  water  surface, 
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(72) 


where 

u*  —  shear  velocity  (=  (rs/p)1/2  (m/sec) 
k  =  von  Karman's  constant  (about  0.41) 
z  =  elevation  above  water  surface  (m) 

zD  =  roughness  length  scale  (or  location  where  W  =  0)  (m) 

From  measured  wind  velocity  profiles,  it  is  possible  to  empirically  determine 
u*  and  zQ,  and,  thus,  CD.  The  drag  coefficient  is  written  using  Equation  49 
and  the  definition  of  shear  velocity, 


( io 


(73) 


where  the  subscript  "10"  denotes  values  evaluated  at  the  10-m  elevation. 

163.  The  logarithmic  wind  velocity  profile  represented  by  Equation  72 
is  valid  for  neutral  (adiabatic)  atmospheric  stability.  Under  non-neutral 
conditions  of  either  stable  and  unstable  thermal  stratification,  the  wind 
profile  departs  from  logarithmic  form.  However,  the  effect  of  thermal 
stratification  diminishes  with  increasing  wind  speed,  more  so  for  stable  than 
unstable  conditions.  Hsu  (1972)  measured  the  vertical  wind  profile  at  three 
locations  in  the  surf  zone  over  a  2 -week  period  and  obtained  nearly  1,500 
15-min  average  profiles.  Analysis  of  the  profile  data  under  adiabatic  and 
onshore  wind  conditions  showed  that  more  than  90  percent  of  the  wind  measure¬ 
ments  were  logarithmic.  This  result  and  others  referenced  by  Hsu  (1972) 
indicates  that  stress  laws  developed  from  logarithmic  wind  profiles  are 
adequate  for  routine  applications  of  engineering  models . 

164.  Most  investigations  of  £7D  transform  the  calculated  drag  coefficient 
to  neutral  or  stable  atmospheric  conditions  using  environmental  data  at  the 
measurement  site.  In  the  simulations  presented  here,  it  is  assumed  that  the 
atmospheric  conditions  are  stable  and  that  values  of  UD  from  the  literature 
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may  be  applied  without  corrections.  Also,  all  wind  velocities  and  drag  coef¬ 
ficients  refer  to  the  10-m  elevation  if  not  otherwise  stated. 

165.  To  establish  the  order  of  magnitude  of  the  drag  coefficient  within 
the  context  of  longshore  current  modeling  and  to  qualitatively  investigate  the 
relationship  between  the  transfer  of  momentum  at  the  air- sea  interface  and  the 
resulting  current,  the  simple  situation  without  lateral  mixing  and  wave 
driving  was  examined.  From  the  momentum  equation,  Equation  55,  it  is  seen 
that  the  wind  stress  is  balanced  by  bottom  friction.  Since  the  momentum 
equation  is  vertically  integrated,  the  wind  effectively  acts  on  the  entire 
water  body  and  generates  an  average  current  velocity  V.  In  the  nearshore 
region  with  comparably  shallow  water,  this  conceptual  model  should  provide  a 
good  first  approximation.  Thus,  equating  the  wind  stress  to  the  bottom 
friction  stress  yields  the  following  relationship  between  wind  speed  and  wind¬ 
generated  current  velocity  (which  is  a  general  velocity  and  not  just  the 
longshore  current  velocity)  in  the  situation  of  equilibrium  momentum  transfer: 


The  ratio  between  the  density  of  air  (pa  =  1.2  kg/m3)  and  density  of  seawater 
(p  =  1.03  103  kg/m3)  is  approximately  1/850.  With  a  ratio  of  C”D/cf  =  0.2,  the 
magnitude  of  the  wind- generated  current  velocity  is  approximately  1.5  percent 
that  of  the  wind  speed,  whereas  for  a  ratio  of  1.0  the  corresponding  value  is 
3.5  percent.  In  the  presence  of  waves,  the  generated  current  for  a  specific 
ratio  of  CD/cf  will  be  smaller  because  of  the  additional  friction  force 
contributed  by  the  wave  bottom  orbital  velocity. 

166.  To  obtain  quantitative  information  on  the  ratio  CD/cf,  field  data 
from  the  FRF  as  described  in  Hubertz  et  al .  (1987)  and  Birkemeier  et  al . 

(1987)  were  used.  These  studies  describe  results  of  a  comprehensive  study  of 
nearshore  waves  and  currents  conducted  during  September  1985  in  which  the 
current  was  measured  at  11  locations  in  the  nearshore  together  with  a  large 
number  of  other  oceanographic  and  meteorological  parameters.  The  longshore 
current  data  set  was  not  used  in  the  present  modeling  on  longshore  currents 
since  most  of  the  current  meters  were  outside  the  surf  zone  during  the 
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measurement  time  (although  a  storm  with  high  waves  passed  the  area  on  12  to 
14  September,  providing  some  potentially  useful  data).  Also,  it  was  difficult 
to  obtain  reliable  information  on  the  incident  wave  angle  defining  conditions 
favorable  to  current  modeling;  the  times  currents  were  generated  only  by  local 
winds  were  not  known.  However,  the  data  do  provide  the  capability  of  evaluat¬ 
ing  wind- generated  currents  in  an  approximate  fashion. 

167.  The  objective  of  the  present  analysis  is  to  obtain  a  rough  estimate 
of  the  ratio  CD/ct  from  wind  and  current  measurements  through  regression 
analysis.  A  data  record  from  7  to  14  September  1985  containing  hourly 
averages  of  wind  speed  and  direction,  longshore  current  velocity,  wave  height, 
and  wave  period  was  analyzed.  On  or  after  14  September,  a  rip  current 
appeared  at  the  location  of  the  main  array  of  the  current  meters,  making  the 
measurements  after  that  date  inappropriate  for  the  present  purpose.  The 
current  was  taken  from  the  most  seaward  current  meter  (Gage  SD25)  in  the  array 
(in  an  average  water  depth  of  6.1  m)  to  minimize  the  influence  of  currents 
generated  by  breaking  waves .  At  this  location  the  current  was  measured  at 
three  vertical  stations,  and  averages  of  those  stations  were  used.  The  wind 
speed  and  direction  were  measured  at  the  end  of  the  FRF  pier  at  elevations  of 
14.02  m  (Gage  682)  and  12.97  m  (Gage  683),  respectively.  No  effort  was  made 
to  transform  the  wind  measurements  to  the  10-m  level  or  correct  for  instabili¬ 
ty  in  the  atmosphere . 

168.  Although  the  current  was  measured  outside  the  breaker  zone  with 
presumably  little  contamination  by  the  wave -generated  longshore  current,  the 
wave  orbital  velocity  was  still  an  important  factor  determining  the  bottom 
friction  because  the  current  meter  was  situated  in  relatively  shallow  water. 
Wave  information  from  Gage  640  in  8  m  of  water  was  used  to  calculate  the 
bottom  orbital  velocity.  Because  reliable  information  on  the  incident  wave 
angle  was  lacking,  a  small-angle  assumption  was  made,  meaning  the  wave  angle 
was  neglected  in  determining  the  bottom  friction.  The  mean  wave  height  was 
chosen  as  a  representative  measure  for  calculating  the  bottom  orbital  velo¬ 
city,  and  the  bottom  friction  was  computed  from  the  Nishimura  (1982)  square- 
wave  approximation  (Equation  61) . 

169.  The  empirical  relationship  between  V[VZ  +  (2/n)umz]1/2  and  |IVjt/y  is 
determined  from  the  trend  in  the  153  available  data  points  obtained  under  the 
aforementioned  assumptions  and  plotted  in  Figure  12,  for  which  Wy  =  |W| sin<j>  is 
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the  longshore  component  of  the  wind.  The  least-squares  fitted  straight  line 
through  the  measurement  points  is  also  shown.  Although  considerable  scatter 
is  present  (partly  due  to  simplifications  introduced  in  the  analysis),  a  clear 
trend  indicates  that  stronger  currents  are  associated  with  stronger  winds. 

The  slope  of  the  straight  line  represents  the  ratio  pjp  CD/cfI  from  which  an 
average  value  can  be  evaluated.  Assuming  a  density  ratio  between  air  and 
seawater  of  1/850,  the  slope  of  the  line  in  Figure  12  gives  CD/cf  -  1.2. 

Thus,  according  to  this  exploratory  analysis,  the  drag  and  bottom  friction 
coefficients  were  approximately  equal  for  the  conditions  studied.  This  result 
is  within  the  range  of  variability  of  the  parameter  values  as  reported  in  the 
literature,  but  higher  than  typical  values  reported.  Most  values  of  CD  given 
in  the  literature  are  considerably  smaller  than  values  of  cf. 


V[V2+(|um)2]  (m/sec)2 


IWlWy  (m/sec)2 


Figure  12.  Wind- induced  longshore  current  versus  wind  velocity 

170.  Although  the  above  analysis  was  crude,  it  shows  that  CD  and  cf 
have  the  same  order  of  magnitude.  The  coefficient  of  determination  for  the 
regression  line  in  Figure  12  is  r2  =  0.82,  indicating  that  more  than  80  per¬ 
cent  of  the  variation  in  the  measured  data  is  explained  by  the  regression 
model.  Inclusion  of  an  intercept  in  the  regression  only  marginally  improved 
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the  fit  (r2  =  0.84),  thus  supporting  the  assumption  that  the  wind  is  the 
dominant  driving  force  for  the  currents  studied  in  this  analysis.  However, 
although  a  simple  linear  regression  model  between  wind  and  current  provides  a 
reasonable  description  according  to  the  data  analysis,  the  cause  and  effect 
process  between  wind  and  current  is  complex.  Other  factors  are  involved,  such 
as  the  contribution  to  the  current  from  obliquely  incident  breaking  waves 
which  were  in  part  or  entirely  generated  by  the  wind  at  the  site. 

171.  The  WAMDI  Group  (1988)  found  that  their  hindcast  of  waves  in  the 
North  Sea  was  sensitive  to  the  value  of  the  drag  coefficient  for  high  wind 
speeds  and  developed  the  following  expression  used  in  NMLONG: 


cDm 


'1.2875  lO'3  , 

(0.8  +  0.065|f/|)  IQ’3 


|jy|  <7.5  m/ sec 
\w\  2:  7 . 5  m/ sec 


(75) 


172.  Hsu*  compared  predictions  of  Equation  75  with  his  field  measure¬ 
ments  in  the  Gulf  of  Mexico  made  at  both  low  wind  speeds  and  at  very  high  wind 
speeds  in  a  hurricane.  Equation  75  fit  his  data  well  through  the  range  of 
wind  speeds.  This  independent  agreement  provides  validation  of  the  WAMDI  drag 
coefficient  formula  at  a  site  with  very  different  wind  and  wave  characteris¬ 
tics,  but  both  sites  are  relatively  shallow-water  basins,  a  situation  occur¬ 
ring  in  nearshore  applications  involving  longshore  currents  in  and  near  the 
surf  zone. 


Numerical  Solution  of  the  Longshore  Current  Equation 

173.  If  the  wave-current  interaction  is  neglected,  iteration  between  the 
wave  and  current  calculations  is  not  necessary,  and  all  wave  properties  can  be 
determined  prior  to  solving  for  the  current.  These  wave  properties  enter  in 
the  current  calculation  through  the  wave  driving  term  dS^/dx,  the  bottom 
friction  term  that  depends  on  um  and  6 ,  and  the  lateral  mixing  term  that 
depends  on  H  and  um.  The  x- component  momentum  equation  yields  the  mean 
displacement  of  the  water  surface  and  thus  the  time -averaged  mean  water  depths 


Personal  Communication,  November  1989,  Dr.  S.  A.  Hsu,  Professor,  Coastal 
Studies  Institute,  Louisiana  State  University,  Baton  Rouge,  LA. 
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used  in  calculation  of  the  current.  The  numerical  grid  for  the  current 
calculation  is  superposed  with  the  wave  model  grid  as  discussed  next. 
Finite-difference  equations 

174.  Figure  13  is  a  definition  sketch  for  the  discretization  and  numeri¬ 
cal  solution  scheme  used  in  the  longshore  current  models.  (A  corresponding 
sketch  for  the  wave  model  is  given  in  Figure  2.)  The  same  approach  is  used 
for  the  linear  and  nonlinear  current  models,  except  that  the  nonlinear  model 
requires  iteration  between  the  current  and  bottom  stress.  The  numerical 
solution  will  first  be  described  for  the  linear  model  and  additional  comments 
given  regarding  modifications  necessary  for  solving  the  nonlinear  model 
equations . 

175.  For  convenience,  the  quantities  A  and  B  are  defined  to  simplify 
notation: 


A  =  ed 


(76) 


B  =  —cfum  (1  +  sin20) 

i 


(77) 


Equation  55  is  written,  using  the  bottom  friction  term  from  Equation  59, 


A 


d2V 

dx2 


+ 


dAdv  -  bv=  1  dSxy 
dx  dx  p  dx 


|^|  Wsincp 


(78) 


Equation  78  is  expressed  in  difference  form  in  Equation  79  for  each  calcula¬ 
tion  cell  running  from  number  2  to  N  -  1  for  a  grid  encompassing  N  cells. 
Boundary  conditions  on  V  must  be  provided  at  cells  1  and  N.  Difference 
equation  representations  of  Equation  78  are  coupled  by  the  mixing  terms  (terms 
in  Equation  78  involving  A)  and  form  a  system  of  N  -  2  equations  with  N  -  2 
unknown  current  velocities  V. 

176.  The  system  of  equations  is  tridiagonal  and  can  be  solved  with  a 
double-sweep  (implicit  solution  scheme)  technique.  In  the  following,  the 
subscript  i  denotes  (as  before)  the  number  of  a  specific  calculation  cell,  and 
a  staggered  grid  will  be  used  with  certain  quantities  evaluated  in  the  middle 
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of  a  cell  and  others  at  the  cell  boundaries,  according  to  Figure  13.  The 
current  velocity  is  given  in  the  middle  of  a  cell  (h-point)  together  with  the 
quantity  Bi,  whereas  A±  and  (Sxy)1  are  evaluated  at  the  boundaries  of  calcula¬ 
tion  cells  (ff-points) .  Because  B  contains  the  wave  orbital  velocity  and  wave 
angle,  its  value  is  interpolated  from  neighboring  ff-points.  A  staggered  grid 
of  this  type,  which  involves  averaging  and  overlapping  of  calculated  quanti¬ 
ties,  enhances  the  stability  of  the  numerical  scheme  by  calculating  deriva¬ 
tives  at  the  most  appropriate  position. 


Figure  13.  Definition  sketch  for  longshore  current  calculation 

177.  The  canonical  form  of  the  difference  equation  for  calculation  cell 
number  i  is , 


-AH1 V1_1  +  BH1Vi  -  CH±VU1  =  RHi 


(79) 


where  the  coefficients  are  defined  as: 


84 


AH-  = 


Aj 


(80) 
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BH,  =  - 


B-  + 


(Ait  i  +  Ai) 


(Ax) : 
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CH, 


ji+i 
(Ax) : 
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1  (S*y)i*l  -  (gw)j 
p  Ax 
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Analogous  to  the  wave  model,  //-point  number  i  lies  on  the  left  boundary  of  the 
corresponding  calculation  cell. 

178.  By  introducing  two  auxiliary  coefficients,  EEi  and  FFit  a  double¬ 
sweep  solution  scheme  may  be  derived.  It  is  assumed  that  a  recursive  solution 
to  the  system  of  equations  exists  of  the  form: 


Vi  =  EEiVi+1  +  FF i 


(84) 


By  substituting  Equation  84,  written  for  index  i-1,  into  Equation  79  the 
following  expressions  are  obtained  for  the  double-sweep  coefficients, 


EE; 


CHi 

DNj 


FF ,  = 


AHiFF^i  +  RHi 
DN, 


(85) 


(86) 
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where 


DN1  =  BH1  -  AH1EE1_1 


(87) 


The  solution  procedure  is  to  determine  EE t  and  FFA  starting  from  1  =  2 
through  the  whole  grid  in  a  forward  sweep  from  knowledge  of  the  shoreward 
boundary  condition  at  grid  point  1-1.  The  current  velocities  are  then 
calculated  in  a  backward  sweep  from  point  N  -  1  by  using  Equation  84  and  the 
seaward  boundary  condition  at  1  —  N. 

Boundary  conditions 

179.  Specification  of  the  seaward  boundary  condition  is  straightforward: 
zero  wave-induced  current  velocity  at  the  end  of  the  grid,  i.e.,  only  a  wind- 
induced  current  if  a  wind  is  blowing.  This  specification  requires  that  the 
seaward-most  calculation  cell  (1  =  N )  be  located  sufficiently  distant  from  the 
break  point  to  truly  have  negligible  wave- induced  current.  If  the  seaward 
grid  boundary  is  located  too  close  to  the  break  point,  the  current  profile 
will  be  artificially  skewed  toward  the  low  side.  One  surf  zone  width  seaward 
of  the  break  point  is  considered  sufficient  distance  for  the  situation  of  no 
wind.  If  a  strong  wind  that  has  a  substantial  cross-shore  component  is 
blowing,  the  seaward-most  calculation  cell  must  be  located  much  farther 
offshore  since  setup  or  setdown  may  extend  to  deeper  water. 

180.  The  most  frequently  specified  shoreward  boundary  condition  in  the 
literature  is  zero  current  at  the  shoreline,  although  this  contradicts  field 
observations,  which  commonly  show  a  longshore  current  in  the  swash  zone.  The 
equations  used  in  the  present  and  other  longshore  current  models  have  limited 
validity  in  the  swash  zone,  where  complex  fluid  motion  and  strong  fluid- 
sediment-bed  interactions  not  described  by  the  conventional  governing  equa¬ 
tions  employed  here  take  place.  To  improve  the  current  prediction  in  shallow 
water  close  to  the  shoreline  or  swash  zone,  a  different  boundary  condition  is 
employed  in  the  present  models.  A  finite-valued  current  at  the  shoreward 
boundary  is  determined  from  the  governing  equation  neglecting  mixing,  which 
then  allows  an  explicit  solution  to  be  found.  The  current  velocity  at  the 
most  shoreward  grid  point  is  determined  for  the  linear  case  as, 
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(88) 


y  _  _  1  ^xy)  2 

1  psx  Ax 


in  which  the  subscript  "1"  refers  to  the  most  shoreward  cell.  By  using 
Equation  84,  the  shoreward  boundary  condition  can  be  expressed  in  terms  of  the 
double-sweep  coefficients.  It  is  found  that  EE1  =  0  and  FF1  = 

Nonlinear  model 

181.  If  the  nonlinear  bottom  friction  term  is  used,  a  double-sweep 
solution  is  not  directly  applicable  since  terms  arise  containing  V2  which 
prevents  deriving  a  simple  recursive  formula  such  as  Equation  84.  However,  if 
an  iterative  technique  is  employed  and  the  quadratic  friction  term  approxi¬ 
mated  using  velocity  values  from  the  preceding  iteration,  the  double-sweep 
scheme  can  still  be  employed.  Initially,  the  current  profile  calculated  with 
a  linear  friction  term  is  used  as  a  first  approximation  in  the  iteration  for 
the  nonlinear  current  profile.  The  only  difference  as  compared  with  the 
previously  discussed  solution  scheme  is  the  quantity  B  (Equation  77) ,  which  is 
replaced  by, 


B  =  cf\Z  + 


w2 


sinz0 


(89) 


where  Z  and  w  are  defined  in  Equations  62  and  63,  respectively,  and  the 
current  velocity  V  needed  in  Z  is  taken  from  the  previous  iteration.  The 
current  profile  is  recalculated  in  this  manner  until  the  maximum  difference 
between  consecutive  iterations  at  any  point  is  less  than  some  predefined 
percentage  of  the  current,  typically  1  percent. 

182.  The  shoreward  boundary  condition  requires  an  iterative  approach  in 
the  nonlinear  model.  (The  seaward  boundary  condition  also  requires  iteration 
if  a  wind- induced  current  is  present.)  Equation  88  is  still  applicable,  but 
the  factor  B1  contains  the  velocity  a  quantity  that  is  required  as  part  of 
the  solution.  A  Newton- Raphson  solution  method  is  used  to  determine  this 
velocity  at  the  shoreward  boundary  with  a  starting  value  from  the  linear 
solution. 
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183.  To  reduce  the  number  of  iterations  required  to  obtain  the  solution 
to  the  nonlinear  equation,  values  obtained  from  the  linear  solution  are 
multiplied  by  a  correction  factor  a  as  given  in  Appendix  A  to  provide  the 
initial  approximation.  Multiplication  by  a  typically  saves  one  to  three 
iterations  out  of  typically  four  to  eight  total  iterations,  a  reduction  of 
about  one  third  to  one  half  the  number  of  required  calculations.  If  both 
waves  and  wind  generate  the  longshore  current  and  if  the  wind  contribution  is 
dominant  in  some  areas  of  the  surf  zone,  the  number  of  iterations  may  increase 
to  10  or  20.  Such  a  situation  is  rare,  because  a  strong  wind  will  ordinarily 
be  accompanied  by  high  waves . 

Wave -current  interaction 

184.  If  the  wave -current  interaction  is  included,  calculations  of  wave 
properties  and  the  current  become  coupled  and  must  be  solved  iteratively. 
Because  the  effect  of  the  current  on  the  wave  height  distribution  was  found  to 
be  small  for  all  current  velocities  tested,  the  calculation  is  started  by 
assuming  no  wave-current  interaction.  After  the  current  computation  is  made, 
the  wave  height  distribution  is  updated,  and  a  new  current  profile  is  ob¬ 
tained.  The  calculation  including  wave -current  interaction  normally  converges 
after  two  or  three  iterations  of  this  outer  loop  to  satisfy  a  1  percent 
tolerance  for  the  wave  height  and  current  at  each  grid  point. 


88 


PART  VI:  VERIFICATION  OF  LONGSHORE  CURRENT  MODEL 


185.  The  present  model  is  a  generalized  version  of  other  analytic  and 
numerical  models  of  the  wave- induced  current  which  have  many  times  been 
demonstrated  to  provide  accurate  descriptions  of  the  form  and  magnitude  of  the 
longshore  current;  therefore,  verification  is  assured  for  the  situation  of 
regular  waves  shoaling  on  beaches  with  monotonically  increasing  depths  and 
with  no  wind  present.  Since  the  model  has  more  capabilities  than  previous 
ones,  the  more  germane  questions  to  be  addressed  in  verification  are: 

a.  What  is. the  difference  in  predictions  between  the  linear  and 
nonlinear  current  models? 

b.  Are  predictions  reasonable  for  the  longshore  current  (and 
waves)  generated  over  a  barred  bottom  profile? 

c.  Are  values  of  empirical  parameters  in  the  current  models 
sufficiently  stable  that  the  model  can  be  used  in  engineering 
applications  with  minimal  or  no  changes? 

d.  Are  the  wind- induced  current  and  effect  of  wind  on  the  mean 
water  level  properly  calculated,  and  does  the  representation  of 
the  wind  produce  reasonable  results? 

e.  How  important  is  the  wave -current  interaction? 

f.  Is  the  model  sufficiently  efficient  to  be  implemented  in  a 
desk- top  computer  environment? 

186.  To  verify  the  applicability  of  the  linear  and  nonlinear  longshore 
current  models  and  to  establish  parameter  values,  the  models  were  tested  to 
reproduce  measurements  from  both  laboratory  and  field.  The  advantage  of  a 
laboratory  environment  is  collection  of  data  under  controlled  conditions  with 
high  resolution  in  time  and  space.  In  the  field  the  forcing  conditions  change 
continuously,  making  it  difficult  to  quantify  the  relationship  between  cause 
and  effect,  and  processes  unaccounted  for  in  the  model  may  be  acting.  How¬ 
ever,  field  data  sets  provide  the  ultimate  test  of  numerical  models  and  are 
thus  of  greatest  importance. 

187.  The  numerical  models  were  verified  using  three  data  sets:  the 
wave- induced  longshore  current,  wave  height,  and  mean  water  level  from  the 
laboratory  experiment  of  Visser  (1982);  the  field  data  on  the  current  from 
Kraus  and  Sasaki  (1979);  and  the  field  data  on  the  waves  and  current  from 
Thornton  and  Guza  (1989) .  Example  calculations  were  also  made  for  a  case 
involving  the  hypothetical  barred  topography  used  by  Ebersole  and  Dalrymple 
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(1980)  in  their  nearshore  modeling  simulation.  For  the  field  data  simula¬ 
tions,  the  random  wave  model  was  used  to  drive  the  longshore  current  model. 

The  influence  of  wind  on  the  predicted  longshore  current  and  mean  water 
surface  elevation  was  also  investigated  in  sensitivity  tests. 

Laboratory  Data  from  Visser  (1982) 

188.  The  numerical  current  model  was  first  applied  to  simulate  four 
cases  from  the  high-quality  measurements  reported  by  Visser  (1982) .  This 
excellent  data  set  provides  a  wide  range  of  wave  and  current  conditions  for 
longshore  current  model  testing.  Visser  conducted  his  experiments  in  a  large 
rectangular  wave  basin  with  a  plane -sloping  beach  using  waves  with  fixed 
height,  period,  and  incident  wave  angle.  Extraordinary  care  was  taken  to 
minimize  boundary  effects  through  use  of  lateral  discharge  controls  to 
minimize  circulation  in  the  basin,  and  it  is  believed  that  artificial  process¬ 
es  associated  with  the  finite  size  of  the  basin  were  reduced  to  a  negligible 
level  in  the  middle  section  of  the  basin.  Seven  cases  were  performed  with 
different  values  of  one  or  more  of  the  quantities  of  wave  height,  period, 
incident  angle,  beach  slope,  and  bed  roughness.  Four  cases  were  chosen  in  the 
present  study  as  representative  of  the  variation  in  the  above-mentioned  vari¬ 
ables,  namely,  Cases  1,  3,  4,  and  7  in  Visser' s  notation. 

189.  Table  3  summarizes  the  four  cases,  and  Appendix  B  gives  a  listing 
of  the  corresponding  wave  height  and  current  data  that  were  compared  with 
model  simulation  results.  The  current  profile  data  in  Appendix  B  were 
obtained  by  averaging  Visser 's  measurements  from  transects  located  in  the 
middle  of  the  wave  basin,  located  far  from  the  basin  boundaries,  where  the 
waves  and  current  were  most  uniform  alongshore. 

190.  Since  the  primary  objective  in  the  present  verification  was  to 
assess  the  performance  of  the  current  and  setup/setdown  calculations  and  not 
the  wave  calculation,  the  breaker  index  yb  in  the  wave  model  was  set  to  the 
value  measured  by  Visser  (1982)  to  eliminate  an  additional  model  parameter.  A 
comparison  is  given  in  Table  3  between  measured  breaker  indices  and  the 
indices  predicted  by  Equation  29,  which  could  be  used  as  an  option  in  the  wave 
model.  The  magnitude  and  the  trend  in  values  with  wave  steepness  is  reason¬ 
ably  well  reproduced.  It  is  interesting  to  note  the  difference  between  the 
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measured  yb  in  Cases  4  and  7 ,  which  have  the  same  wave  conditions  and  beach 
slope,  differing  only  with  respect  to  the  beach  roughness.  Visser  did  not 
measure  the  mean  water  level  in  Case  7,  assuming  that  it  was  the  same  as  in 
Case  4. 

191.  In  the  wave  calculation,  the  two  empirical  parameters  in  the 
breaker  decay  model  were  assigned  the  recommended  values  of  k  -  0.15  and 
T  “  0.40.  The  calculation  grid  was  extended  to  the  horizontal  portion  of  the 
wave  basin  to  make  the  current  computation  independent  of  the  location  of  the 
most  seaward  grid  point. 


Table  3 

Summary  of  Experimental  Cases  from  Visser  (1982-) 


Case 

H* 

T 

e 

Beach 

No. 

ra 

sec 

deg 

Slope 

7meas 

^calc 

Comment 

1 

0.072 

2.01 

31.1 

0.101 

1.00 

1.09 

Large  wave  angle 

3 

0.089 

1.00 

15.4 

0.101 

0.85 

0.95 

Large  wave  height 

4 

0.078 

1.02 

15.4 

0.050 

0.83 

0.83 

Smooth  bottom 

7 

0.078 

1.02 

15.4 

0.050 

0.74 

0.83 

Rough  bottom 

*  Wave  conditions  refer  to  the  horizontal  portion  of  the  basin. 


192.  The  two  empirical  parameters  in  the  current  model,  the  bottom 
friction  coefficient  and  the  coefficient  A  in  the  lateral  mixing  term,  were 
optimized  to  provide  the  best  fit  to  the  data.  As  a  criterion  for  the  best 
fit,  the  sum  of  squares  of  the  difference  between  the  measured  and  predicted 
current  velocities  was  used.  Both  the  linear  and  nonlinear  bottom  friction 
terms  were  used  in  the  simulations  to  evaluate  the  friction  formulation.  As 
described  by  Visser  (1982) ,  the  assumption  of  a  weak  current  with  respect  to 
the  bottom  orbital  velocity  was  to  a  great  extent  violated  in  his  experiments. 
Thus,  in  principle,  the  linearized  friction  model  should  produce  inferior 
predictions  of  the  current  distribution.  However,  as  will  be  shown,  the 
difference  between  optimal  modeling  results  for  the  linear  and  nonlinear 
models  is  not  significant. 
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193.  The  modeling  results  for  the  four  cases  are  shown  in  Figures  14  to 
17.  Each  figure  contains  the  measured  and  calculated  wave  height  distribution 
across -shore  and  measured  and  calculated  current  profiles  calculated  using  the 
linear  and  nonlinear  models.  Table  4  summarizes  the  modeling  efforts  in  terms 
of  the  optimal  parameter  values  for  the  linear  and  nonlinear  models.  In  the 
following,  a  summary  of  results  for  each  case  is  given. 

Table  4 

Optimal  Model  Parameters  for  the  Visser  Q982)  Cases 


Case  No. 

(Cf)lin 

A]_in 

( Cf )non 

Anon 

1 

0.014 

0.60 

0.009 

0.20 

3 

0.012 

0.30 

0.007 

0.15 

4 

0.009 

0.60 

0.005 

0.30 

7 

0.019 

0.60 

0.014 

0.50 

Case  1  (Figure  14) 

194.  This  case  proved  the  most  difficult  to  reproduce  because  the 
maximum  in  the  measured  current  profile  was  located  far  shoreward  of  the  break 
point  in  the  inner  half  of  the  surf  zone,  shoreward  of  the  region  where 
maximum  decay  in  the  wave  height  occurred.  Visser  (1984a,  b)  developed  a 
numerical  model  of  the  current  which  incorporated  the  assumption  that  the 
energy  dissipation  starts  at  the  plunge  point,  not  at  the  break  point. 

However,  this  concept  introduces  at  least  one  additional  parameter  in  a  model, 
the  plunge  distance,  which  is  related  to  the  breaker  type  or  how  the  transfer 
of  kinetic  turbulent  energy  occurs. 

195.  Smith  and  Kraus  (1990)  found  that  the  plunge  distance  for  barred 
profiles  was  approximately  half  that  for  a  plane  slope  with  the  same  incident 
wave  conditions.  Breaker  type  was  also  different  between  barred  profile  and 
plane  slope.  Reliable  prediction  of  plunge  distance  is  at  present  difficult, 
but  recent  progress  holds  promise  that  advances  will  be  made  for  use  in  engi¬ 
neering  models . 
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Calculated  Height  Measured  Height  Linear  Model  Nonlinear  Model  Measured  Current 

.  *  .  .  □ 

a.  Wave  height  and  longshore  current  distributions 

Mean  Water  Surface  Elev.  (m) 

0.05 

0.04 

0.03 

0.02 

0.01 

0 

-0.01 

Figure  14.  Model  simulations  for  Visser  (1982)  Case  1 


b.  Mean  water  surface  elevation 
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196.  Model  simulations  for  Case  1  thus  have  the  maximum  in  the  current 
profile  located  seaward  of  the  measured  maximum,  although  the  wave  height 
distribution  is  well  predicted,  especially  in  the  surf  zone.  The  results 
shown  in  Figure  14a  illustrate  the  applicability  of  the  breaker  decay  model 
modified  to  include  wave  refraction.  The  prediction  of  the  nonlinear  model 
had  an  overall  inferior  least  squares  fit  as  compared  with  the  linear  model. 
However,  the  tail  of  the  current  profile  outside  the  surf  zone  is  better 
reproduced  by  the  nonlinear  model.  The  general  trend  for  the  nonlinear  model 
is  to  place  the  maximum  in  the  current  somewhat  more  seaward  than  does  the 
linear  model,  but  to  have  a  more  rapid  decay  in  current  magnitude  seaward  of 
the  break  point. 

197.  By  employing  a  smaller  value  of  the  mixing  coefficient,  a  current 
distribution  more  similar  to  the  measured  distribution  was  obtained  with  the 
models,  but  the  main  body  of  the  distribution  was  located  too  far  seaward. 

This  location  produces  a  larger  least  squares  error  compared  with  the  higher 
value  of  the  mixing  coefficient.  The  optimal  mixing  coefficient  for  the 
nonlinear  model  is  significantly  smaller  than  for  the  linear  model,  which  is 
typical  for  all  cases  studied.  The  optimum  for  the  parameter  values  is 
located  in  a  rather  flat  region  in  mixing/friction  coefficient  space,  implying 
that  slight  changes  in  the  optimum  parameter  values  cause  only  small  changes 
in  the  sum  of  squares.  As  an  example,  increasing  or  decreasing  the  mixing 
coefficient  by  50  percent  from  the  optimal  value  of  0.2  for  the  nonlinear 
model  produced  a  change  in  the  sum  of  squares  of  about  1  percent.  The 
sensitivity  to  the  friction  factor  is  somewhat  greater  since  this  parameter 
determines  the  magnitude  of  the  current  distribution.  Overall,  predictions  of 
the  longshore  current  model  were  not  particularly  sensitive  to  small  variation 
in  the  two  empirical  parameters  cf  and  A  about  the  optimum. 

198.  The  calculation  of  the  mean  water  surface  elevation  is  well  repro¬ 
duced  (Figure  14b),  particularly  in  the  inner  surf  zone.  The  point  of 
calculated  maximum  setdown  occurs  somewhat  seaward  of  the  measured  location, 
causing  the  setup  to  begin  more  seaward.  The  general  trend  in  the  setup  and 
setdown  is  well  reproduced. 
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Case  3  (Figure  IS') 

199.  The  current  measurements  are  well  reproduced  by  the  models, 
although  the  position  of  the  calculated  maximum  lies  slightly  seaward  of  the 
measured  maximum.  The  wave  height  calculation  agrees  with  the  trend  in  the 
measurements  up  to  the  break  point  and  for  some  distance  into  the  surf  zone, 
but  the  measured  wave  height  decay  is  steeper  closer  to  shore.  A  better  fit 
could  have  been  achieved  by  optimizing  parameter  values  in  the  breaker  decay 
model.  Compared  with  the  linear  model,  the  nonlinear  current  model  required  a 
smaller  valued  optimum  mixing  coefficient,  giving  a  maximum  in  the  current 
located  more  seaward  than  for  the  linear  model.  The  linear  model  calculation 
has  an  overall  better  visual  fit  than  the  nonlinear  model.  The  optimal 
friction  coefficient  for  the  nonlinear  model  is  considerably  smaller  than  the 
linear  model  (Table  4) . 

200.  The  calculated  point  of  maximum  setdown  again  lies  seaward  of  the 
measured  point,  and  the  more  seaward  start  of  setup  in  the  model  causes  the 
mean  water  surface  elevation  to  be  slightly  overestimated  in  the  surf  zone. 
Case  4  (Figure  16) 

201.  In  this  case  the  bottom  was  smoother  than  for  Cases  1  and  3,  and 
the  location  of  the  maximum  current  is  better  reproduced  by  the  current  models 
than  in  those  cases.  The  measured  wave  height  distribution  decays  more 
rapidly  than  the  calculated  distribution,  as  in  Case  3.  Overall,  the  wave  and 
current  models  well  simulate  the  respective  measurements. 

202 .  The  measured  setdown  was  remarkably  flat  in  this  case ,  whereas  the 
numerical  model  gave  the  typical  sharp  minimum.  The  measured  setup  was  well 
predicted  by  the  model,  including  the  value  at  the  mean  shoreline. 

Case  7  (Figure  l?) 

203.  The  bottom  roughness  was  increased  for  this  case  as  compared  with 
Case  4,  and  a  larger  value  of  the  optimum  friction  coefficient  (Table  4)  was 
required  in  the  models.  Best  agreement  between  measurements  and  model 
predictions  for  all  four  Visser  (1982)  data  cases  was  found  for  this  case  with 
the  simulated  maximum  in  the  current  lying  almost  on  top  of  the  measured 
maximum.  The  wave  height  distribution  showed  the  least  agreement  for  the  four 
cases,  with  the  predicted  decay  beginning  shoreward  of  the  measured  decay. 
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Current  Velocity  (m/sec)  Wave  Height  (m) 


Distance  Offshore  (m) 

Calculated  Height  Measured  Height  Linear  Model  Nonlinear  Model  Measured  Current 

.  *  .  .  □ 


a.  Wave  height  and  longshore  current  distributions 


Distance  Offshore  (m) 


b.  Mean  water  surface  elevation 
Figure  15.  Model  simulations  for  Visser  (1982)  Case  3 
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Calulated  Height  Measured  Height  Linear  Model  Nonlinear  Model  Measure  Current 

*  .  .  □ 

a.  Wave  height  and  longshore  current  distributions 


b.  Mean  water  surface  elevation 
Figure  16.  Model  simulations  for  Visser  (1982)  Case  4 
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Current  Velocity  (m/sec)  Wave  Height  (m) 


Calculated  Height  Measured  Height  Linear  Model  Nonlinear  Model  Measured  current 
- - -  .  .  □ 

Figure  17.  Model  simulations  for  Visser  (1982)  Case  7 


This  offset  in  the  wave  height  distribution  promoted  the  better  agreement  in 
the  current  calculations,  similar  to  Visser' s  (1984a,  b)  use  of  the  plunge 
point  rather  than  the  break  point  to  start  wave  dissipation.  The  difference 
between  the  optimum  mixing  coefficient  in  the  linear  and  nonlinear  models  was 
smaller  than  in  the  previous  cases.  Accordingly,  the  tail  of  the  current 
distribution  for  the  nonlinear  model  is  located  slightly  seaward  of  the  linear 
model  as  opposed  to  the  other  cases . 

Summary 

204.  In  conclusion,  the  following  summarizes  model  validation  for  the 
wave  height  distribution,  longshore  current  profile,  and  water  surface 
elevation  using  laboratory  measurements  of  Visser  (1982) : 

a.  The  longshore  current  profile  was  well  reproduced  by  the  linear 
and  nonlinear  models,  but  the  maximum  in  the  measured  current 
is  located  farther  inshore  than  predicted  by  the  numerical 
models,  probably  related  to  the  wave  breaker  type  and  location 
of  the  plunge  point,  which  are  not  parameters  in  the  models. 
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b.  The  wave  height  distribution  was  well  predicted  with  the 
breaker  decay  model  modified  to  include  refraction  and  using 
previously  established  parameter  values  for  all  simulations  (no 
calibration) . 

c.  The  overall  fit  of  the  linear  model  was  somewhat  better  than 
for  the  nonlinear  model  even  though  the  weak- current  approxima¬ 
tion  was  violated. 

d.  The  optimum  friction  coefficient  for  the  nonlinear  model  was 
significantly  smaller  for  all  cases,  on  the  average  of  about 
35  percent  less  than  that  of  the  linear  model  (Table  4) . 

Values  for  the  friction  coefficient  for  the  linear  model  were 
0.012  (Cases  1,  3,  and  4)  and  0.019  (Case  7),  whereas  respec¬ 
tive  values  for  the  nonlinear  model  were  0.0070  and  0.014. 

e.  The  optimum  mixing  coefficient  for  the  nonlinear  model  was  sig¬ 
nificantly  smaller  for  all  cases,  on  the  average  of  about 

45  percent  less  than  that  of  the  linear  model.  The  average 
optimum  mixing  coefficient  was  0.5  for  the  linear  model  and  0.3 
for  the  nonlinear  model. 

f .  The  optimum  mixing  and  friction  parameter  values  in  the  current 
models  were  located  in  relatively  flat  regions,  implying  that 
model  predictions  are  insensitive  to  slight  changes  around  the 
optimum. 

g.  The  measured  change  in  mean  water  surface  elevation  produced  by 
the  obliquely  incident  waves  was  well  reproduced  by  the  model. 

205.  The  linear  current  model  gave  a  somewhat  overall  better  fit  than 
the  nonlinear  model.  This  result  is  not  considered  to  reflect  a  physically 
better  description  of  the  phenomenon,  but  indicates  that  the  linear  model  is 
more  flexible  in  calibration.  Smaller  values  of  the  mixing  coefficient 
obtained  with  the  nonlinear  model  could  be  interpreted  as  less  necessity  for 
smoothing,  implying  a  better  physical  description. 


Field  Data  from  Kraus  and  Sasaki  (1979) 

206.  Kraus  and  Sasaki  (1979)  measured  the  longshore  current  profile 
along  seven  transects  on  a  sandy  beach  facing  the  Japan  Sea.  The  measurements 
were  made  simultaneously  on  a  transect  by  divers  positioned  at  5-m  intervals. 
After  completing  measurements  on  one  transect,  the  divers  moved  to  the  next 
transect,  located  25  m  away.  The  current  was  measured  by  timing  the  movement 
of  almost  neutrally  buoyant  floats  with  current- sensing  fins  located  at  about 
middepth;  the  floats  were  tethered  to  the  wrists  of  the  divers  by  a  2 -m- long 
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line.  An  average  current  velocity  was  formed  of  three  successive  measurements 
at  each  location.  The  incident  waves  arrived  as  clean  swell,  and  repeated 
visual  measurements  of  the  breaking  waves  using  a  graduated  rod  gave  a  signif¬ 
icant  height  of  1.0  m,  period  of  4.1  sec,  and  angle  of  9  deg.  The  water  depth 
at  each  sampling  location  was  measured  by  rod  and  transit.  For  the  numerical 
simulations  performed  in  the  present  study,  average  values  of  the  water  depth 
and  current  velocity  were  used.  Appendix  B  contains  a  summary  of  the  measure¬ 
ments.  The  bottom  topography  had  a  step-type  shape,  and  secondary  breaking 
was  observed  far  shoreward  of  the  main  break  point  located  in  about  1-m  water 
depth. 

207.  As  shown  in  Figure  18,  the  measured  current  profile  for  the  average 
bottom  bathymetry  had  a  peak  located  somewhat  inshore  of  the  main  break  point 
but  was  otherwise  uniform  throughout  the  surf  zone.  A  flat  profile  is 
difficult  to  reproduce  in  a  model  using  only  monochromatic  waves  unless  strong 
mixing  is  applied;  in  contrast,  inclusion  of  wave  reformation  and  secondary 
breaking  can  improve  modeling  results  in  such  a  situation.  A  secondary  break 
point  produces  a  second  local  maximum  in  the  current  velocity  inshore, 
producing  the  uniform- like  current  profile  inshore  between  the  two  break 
points.  On  the  other  hand,  since  the  depth  monotonically  increases  with 
distance  offshore,  wave  reformation  will  not  occur  unless  forced  by  a  lower 
value  of  T  in  Equation  24  (Dolan  and  Dean  1984,  Larson  and  Kraus  1989). 

208.  If  a  distribution  of  wave  heights  is  employed  to  drive  the  long¬ 
shore  current  model,  breaking  will  occur  more  uniformly  over  the  surf  zone  and 
produce  a  more  uniform  current  profile.  As  described  in  Part  III,  one 
procedure  for  obtaining  such  a  distribution  is  to  randomly  select  wave  heights 
from  a  Rayleigh  distribution  specified  at  some  depth  offshore  and  transform 
the  individual  waves  with  the  breaker  decay  model.  At  each  point  across  the 
shore  the  probability  density  function  of  the  wave  height  can  be  constructed 
from  all  the  individual  transformed  wave  heights. 

209.  An  intuitive  method  for  obtaining  the  mean  longshore  current 
distribution  is  to  calculate  a  current  corresponding  to  each  randomly  picked 
wave  height  and  average  these  distributions.  Although  statistically  appeal¬ 
ing,  this  approach  is  questionable  from  a  physical  point  of  view  since  a 
single  wave  cannot  develop  a  steady  current;  a  large  number  of  consecutive 
identical  waves  are  required.  However,  for  a  long  averaging  interval  or 
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Figure  18.  Simulations  for  the  Kraus  and  Sasaki  (1979)  data; 

linear  model 


measurement  period  with  little  correlation  between  consecutive  waves,  a  Monte- 
Carlo  simulation  technique  is  considered  a  valid  procedure  for  obtaining  an 
average  current  value.  Concerning  the  equation  for  calculating  the  mean  water 
surface  elevation  (Equation  15),  the  question  can  also  be  raised  about  the 
physical  meaning  of  an  average  based  on  different  individual  waves.  Neverthe¬ 
less,  as  an  engineering  expedient,  this  procedure  is  used  in  NMLONG. 

210.  Kraus  and  Sasaki  (1979)  did  not  measure  the  wave  height  across  the 
shore,  but  only  at  the  break  point.  In  the  present  work,  the  significant  wave 
was  transformed  seaward  to  about  4-m  depth  using  linear-wave  theory,  where  no 
waves  in  a  distribution  would  break  and  the  wave  calculation  could  be  started. 
The  rms  wave  height  was  determined,  and  wave  heights  were  randomly  picked  from 
a  Rayleigh  distribution  defined  by  this  wave  height.  It  is  noted  that  since 
conditions  at  the  break  point  were  used  to  derive  the  offshore  wave  height, 
uncertainty  was  introduced  in  the  wave  height  calculations.  The  incident  wave 
period  and  angle  were  kept  constant  throughout  the  simulation,  together  with 
the  breaker  index  yb.  Randomization  of  wave  direction  and  period  in  a  joint 
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probability  distribution  would  be  an  extension  of  this  procedure;  for  the 
present,  the  assumption  was  made  that  a  variation  in  wave  height  captures  the 
main  effect  of  wave  randomness.  A  constant  breaker  index  yb  of  individual 
waves  was  chosen  as  opposed  to  Equation  29  because  this  equation  was  developed 
for  monochromatic  waves  and  involves  a  moderately  strong  dependence  on  wave 
period.  Since  in  the  present  approach  wave  period  was  constant  and  wave 
height  varied,  low  wave  steepnesses  will  occur  that  produce  unrealistically 
high  values  of  yb  by  Equation  29.  Thus,  a  constant  value  (yb  =  1.0)  was 
chosen  since  this  was  the  measured  breaker  index  for  the  visually  significant 
wave  at  the  break  point. 

Linear  model 

211.  Both  the  linear  and  nonlinear  current  models  were  tested  against 
the  data,  and  optimum  parameter  values  were  sought  with  respect  to  minimizing 
the  sum  of  squares  as  described  for  the  laboratory  tests  above.  The  current 
distribution  was  taken  as  an  average  for  500  waves,  giving  a  stable  estimate 
of  the  mean  current  at  points  across  the  shore.  The  simulations  started  by 
seeking  the  optimum  value  of  the  friction  coefficient  for  a  value  of  the 
mixing  coefficient  corresponding  to  the  average  value  (A  =  0.5)  determined 
from  the  laboratory  data  of  Visser  (1982).  For  the  linear  model,  a  friction 
coefficient  of  0.0042  gave  the  best  agreement  with  the  field  measurements. 

This  value  is  significantly  lower  than  what  was  obtained  for  the  laboratory 
measurements  of  Visser  (1982),  but  close  to  the  value  cf  =  0.0061  obtained  by 
Kraus  and  Sasaki  (1979)  using  an  analytic  model  of  the  longshore  current  and 
fitting  the  location  of  the  maximum  in  the  profile.  Figure  18  shows  the 
measured  and  calculated  current  distributions  for  the  linear  model  together 
with  the  calculated  rms  wave  height  across  the  shore. 

212.  In  a  true  optimization  with  the  friction  coefficient  and  mixing 
parameter  allowed  to  range  freely,  a  significantly  higher  value  of  A  was  ob¬ 
tained,  namely,  A  =  1.2,  whereas  cf  remained  the  same.  However,  these 
calculations  with  higher  valued  mixing  coefficients  showed  that  the  "improve¬ 
ment"  in  an  rms  sense  of  defining  error  was  an  artifact  of  extreme  flattening 
of  the  entire  distribution  and  prolongation  of  the  tail  in  deeper  water  where 
no  measurements  were  available.  In  contrast,  extrapolation  of  the  measure¬ 
ments  to  deeper  water  shows  a  sharp  decrease  in  the  current.  Thus,  additional 
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measurement  points  in  this  area  would  degrade  the  solution  with  greater  mixing 
and  translate  the  optimum  mixing  coefficient  toward  a  lower  value. 

213.  As  noted  by  Battjes  (1972)  and  Thornton  and  Guza  (1986),  the  need 
for  lateral  mixing  to  smooth  the  current  distribution  near  the  break  point  is 
reduced  for  random  waves  as  compared  with  monochromatic  waves.  There  is  a 
difference,  however,  between  obtaining  a  smooth  velocity  profile  and  obtaining 
the  distribution  that  shows  best  agreement  with  measurements.  Since  Thornton 
and  Guza  (1986)  found  little  difference  in  their  random  longshore  current 
model  whether  or  not  lateral  mixing  was  incorporated  in  simulations  of  their 
field  data,  calculations  were  carried  out  with  the  present  model  with  a  small 
value  of  A  to  investigate  this  effect.  Using  the  linear  current  model,  the 
sum  of  squares  increased  significantly  for  the  weak-mixing  case  (A  =  0.05), 
and  the  value  of  the  optimal  friction  coefficient  increased  to  cf  —  0.0050. 

The  conclusion  is  that  lateral  mixing  should  be  included  in  the  present  model 
to  obtain  good  agreement  with  field  measurements . 

Nonlinear  model 

214.  Predictions  with  the  nonlinear  model  were  also  compared  with  the 
current  measurements  (Figure  19),  starting  the  series  of  simulations  using  the 
value  of  the  mixing  coefficient  given  by  the  laboratory  data,  that  is, 

A  =  0.3.  A  cf-value  of  0.0022  gave  the  best  least  squares  agreement  and  had  a 
sum  of  squares  about  40  percent  greater  than  the  linear  model.  Similar  to 
the  linear  model  case,  the  optimum  fit  occurred  for  a  larger  mixing  coeffi¬ 
cient  value,  but  the  difference  in  the  sum  of  squares  was  only  slightly  less, 
because  the  cf  and  A  surface  was  relatively  flat. 

Summary 

215.  The  linear  and  nonlinear  models  predicted  the  general  two-peaked 
form  of  the  longshore  current  field  measurements  of  Kraus  and  Sasaki  (1979) , 
providing  a  much  superior  fit  than  possible  with  analytic  models.  The 
calculated  current  profile  was  somewhat  low  in  the  shoreward  part  of  the  surf 
zone  and  too  high  seaward  of  the  maximum  current.  The  fit  could  be  improved 
by  increasing  the  value  of  the  breaker  index  yb  to  1.1,  but  it  was  desired  to 
minimize  the  number  of  free  parameters.  Using  a  value  for  the  mixing  coeffi¬ 
cient  determined  by  the  laboratory  data  in  finding  the  optimal  friction 
coefficient  gave  reasonable  results.  In  the  laboratory  experiments,  mono¬ 
chromatic  waves  were  used,  thus  facilitating  evaluation  of  the  mixing 
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Figure  19.  Simulations  for  the  Kraus  and  Sasaki  (1979)  data; 

nonlinear  model 

coefficient  under  ideal  conditions.  For  the  field  comparison,  where  the 
method  of  averaging  many  current  distributions  was  employed  and  each  distribu¬ 
tion  generated  by  different  wave  conditions,  inferences  about  appropriate 
values  of  A  through  calibration  are  probably  less  reliable  than  for  the 
Visser  (1982)  laboratory  data. 

216 .  The  technique  of  averaging  a  large  number  of  current  distributions 
obtained  from  randomly  generated  wave  heights  was  successful  in  reproducing  a 
uniform  current  distribution  across  the  beach.  In  general,  it  was  possible  to 
achieve  a  better  fit  to  the  measurements  with  the  linear  model  as  compared 
with  the  nonlinear  model  for  the  same  reasons  as  discussed  for  the  laboratory 
data . 

Field  Data  from  Thornton  and  Guza  (1986) 


Previous  work 

217.  To  further  verify  the  applicability  of  the  longshore  current  models 
to  describe  field  conditions  and  to  establish  model  parameter  values  relevant 


104 


to  the  field,  the  wave  and  longshore  current  data  of  Thornton  and  Guza  (1986, 
1989)  were  used.  The  random  wave  model  is  tested  by  predicting  the  transfor¬ 
mation  of  the  rms  wave  height  measured  at  fixed  points  across  the  profile  in 
the  experiment.  The  data  were  collected  during  a  field  experiment  conducted 
at  Leadbetter  Beach,  Santa  Barbara,  California,  as  a  part  of  the  Nearshore 
Sediment  Transport  Study.  At  this  site  the  bottom  contours  are  relatively 
straight  and  parallel  (and  with  no  appreciable  bars  and  troughs) ,  thus  making 
the  data  suitable  for  applying  a  current  model  that  assumes  uniform  conditions 
alongshore.  The  data  used  for  the  simulations  in  this  report  derive  primarily 
from  Thornton  and  Guza  (1989) . 

218.  Four  cases  were  used  in  the  present  study  originating  from  data 
collected  on  3-6  February  1980.  These  cases  will  be  referred  to  in  this 
report  by  the  date  of  data  collection.  Table  5  provides  a  summary  of  the 
general  wave  characteristics  for  each  case,  and  a  complete  listing  of  the  data 
is  given  in  Appendix  B.  The  data  from  2  February  presented  in  Thornton  and 
Guza  (1986)  were  excluded  from  the  present  model  tests  with  reference  to  "the 
lack  of  observations  at  the  nearshore"  stated  by  Thornton  and  Guza  (1986). 

Wave  information  presented  in  Thornton  and  Guza  (1986)  was  used  and  provided 
wave  characteristics  in  about  the  4-m  water  depth  as  opposed  to  information 
given  in  Wu,  Thornton,  and  Guza  (1985)  which  refers  to  a  water  depth  of  about 

9  m.  The  bathymetry,  measured  rms  wave  height,  and  measured  current  velocity 
across-shore  are  listed  in  Thornton  and  Guza  (1989). 

219.  Several  model  simulation  studies  have  used  the  data  set  of  Wu, 
Thornton,  and  Guza  (1985)  and  Thornton  and  Guza  (1986,  1989).  For  example, 

Wu,  Thornton,  and  Guza  used  the  data  to  test  a  2-D  model  based  on  the  steady- 
state  momentum  equations  (Equations  32  and  33)  with  a  linearized  bottom 
friction  term  and  without  wind  forcing.  The  rms  wave  height  was  taken  as  the 
representative  measure  and  transformed  as  a  monochromatic  wave.  In  the  surf 
zone,  the  wave  height  was  simply  taken  as  being  proportional  to  the  water 
depth  =  7bh) .  In  another  use  of  the  data,  Thornton  and  Guza  (1986) 

verified  a  1-D  random  wave  and  current  model  based  on  ensemble -averaged 
governing  equations.  The  percentage  of  broken  waves  was  assumed  to  be 
proportional  to  the  ratio  between  wave  height  and  water  depth.  Neither  of 
these  modeling  studies  described  wave  reformation  in  the  surf  zone. 
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Table  5 

Wave  Conditions  from  the  Thornton  and  Guza  (1986)  Data 


^rms 

T 

6 

h 

Case 

m 

sec 

deg 

m 

3  Feb 

0.55 

14.3 

7.8 

4.0 

4  Feb 

0.56 

14.2 

9.0 

3.8 

5  Feb 

0.45 

12.8 

8.4 

3.6 

6  Feb 

0.26 

11.6 

8.3 

3.5 

Present  modeling 

220.  In  the  following,  a  short  description  is  given  of  the  modeling 
procedure,  and  Figures  20  to  23  present  the  simulation  results  for  the  wave 
and  longshore  current  models  together  with  the  measurements.  For  all  cases, 
the  parameter  values  k  =  0.15  and  T  =  0.40  were  used  in  the  breaker  decay 
model,  and  a  breaker  index  of  yb  =  0.78.  Only  values  of  cf  and  A  were  changed 
in  the  longshore  current  models  to  find  an  optimal  solution.  The  cell  length 
was  1  m,  and  the  numerical  grid  extended  to  a  point  where  no  influence  on  the 
current  distribution  was  expected  (6-m  water  depth) .  To  start  the  simulation, 
offshore  wave  measurements  were  transformed  to  the  seaward  end  of  the  grid 
using  linear-wave  theory  and  neglecting  setdown. 

221.  The  general  approach  in  the  simulations  was  to  assume  a  Rayleigh 
distribution  for  the  wave  height  at  the  seaward  end  of  the  grid  and  randomly 
select  500  waves  for  transformation  across  the  beach.  The  rms  wave  height  at 
each  point  across -shore  was  then  obtained  together  with  the  current  velocity 
for  each  wave,  and  the  current  velocities  for  the  500  waves  were  averaged. 

The  optimum  parameter  values  in  the  current  model  were  first  estimated  using 
the  mixing  value  from  the  laboratory  simulations  and  determining  the  friction 
coefficient  yielding  the  smallest  sum  of  squares.  Thereafter,  the  true 
optimum  was  sought,  and  the  sensitivity  of  the  model  to  changes  in  the  parame¬ 
ter  values  established.  The  response  of  the  current  model  to  very  small 
mixing  values  was  also  investigated.  Table  6  summarizes  the  simulation 
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results  for  optimum  parameter  values  using  mixing  values  from  the  laboratory- 
simulations  . 


Table  6 

Optimum  Parameter  Values  for  the  Thornton 


and 

Guza  ('1986') 

Data 

Linear 

model 

Nonlinear 

Model 

Case 

C£ 

A 

cf 

A 

3  Feb 

0.0060 

0.50 

0.0050 

0.30 

4  Feb 

0.0050 

0.50 

0.0035 

0.30 

5  Feb 

0.0040 

0.50 

0.0030 

0.30 

6  Feb 

0.0040 

0.50 

0.0035 

0.30 

222.  3  February  case.  This  was  the  most  difficult  case  to  reproduce  and 
had  the  largest  deviation  between  calculation  and  measurements.  The  case 
clearly  displays  differences  between  model  predictions  and  measurements  that 
occurred  to  some  extent  in  all  simulations  with  this  data  set.  Figure  20a 
shows  the  simulated  and  measured  rms  wave  height  across-shore ,  and  Figure  20b 
gives  the  calculated  and  measured  current  distribution.  The  wave  height 
across  the  profile  is  rather  well  predicted. 

223.  The  predicted  current  profile  for  this  case  has  an  overall  shape  in 
agreement  with  the  measurements,  but  is  translated  shoreward.  This  tendency 
was  true  for  the  linear  and  nonlinear  models,  with  both  models  giving  about 
the  same  sum  of  squares  error  estimate.  The  trend  is  opposite  to  what  was 
noted  for  the  laboratory  cases  of  Visser  (1982),  for  which  the  predicted 
maximum  in  the  current  was  always  located  seaward  of  the  measured  maximum. 
Thus,  in  this  case  introduction  of  a  plunge  point  concept  as  used  by  Visser 
(1984a,  b)  would  degrade  the  fit.  Optimum  values  of  the  friction  coefficient 
for  the  linear  and  nonlinear  models  were  0.0060  and  0.0050,  respectively,  for 
the  mixing  parameter  set  to  values  determined  in  the  laboratory  tests.  The 
optimum  for  freely  varying  combinations  of  cf  and  A  was  located  at  a  higher 
value  of  the  mixing  coefficient  and  was  in  a  very  flat  region. 
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Wave  Height  (m) 


a.  Wave  height  distribution  across -shore 


Current  Velocity  (m/sec) 


b.  Longshore  current  profile 

Figure  20.  Simulations  for  Thornton  and  Guza  (1986)  data 

(3  February) 
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224.  4  February  case.  The  optimum  friction  coefficient  was  smaller  than 
for  3  February  case  with  cf  =  0.0050  and  0.0035  for  the  linear  and  nonlinear 
models  using  optimum  mixing  values  from  the  respective  laboratory  simulations. 
The  comparison  between  linear-  and  nonlinear-model  calculations  and  measured 
wave  height  distributions  and  current  profiles  is  shown  in  Figures  21a  and  b. 
The  wave  height  distribution  is  well  reproduced,  although  it  is  overestimated 
in  the  shoreward  half  of  the  surf  zone .  The  maximum  in  the  calculated  current 
profile  was  located  shoreward  of  that  in  the  measurements.  A  small  value  of 
the  mixing  coefficient  gave  an  inferior  fit,  but  the  increase  in  the  sum  of 
squares  was  not  as  marked  as  for  the  Kraus  and  Sasaki  (1979)  data  set  simula¬ 
tions.  The  linear  and  nonlinear  models  gave  an  almost  identical  fit  according 
to  the  sum  of  squares . 

225.  5.  February  case.  Model  predictions  and  measurements  are  compared 
Figures  22a  and  b.  The  tails  of  the  velocity  profiles  show  good  agreement, 

but  the  calculated  peak  currents  lie  closer  to  shore  than  the  measurements , 
and  the  peak  magnitudes  are  overpredicted.  Optimum  values  of  the  friction 
coefficient  were  0.0040  and  0.0030  for  the  linear  and  nonlinear  models, 
respectively,  with  mixing  coefficients  as  determined  in  the  laboratory  tests. 
Thus,  smaller  values  for  the  friction  coefficient  were  obtained  than  for  the 
previous  cases  in  the  data  set. 

226.  6_ February  case.  For  this  case,  the  wave  height  was  quite  small 
(tfrms  =  0.26  m) ,  and  the  associated  current  was  weakest  of  the  four  cases 
studied.  Figures  23a  and  b  show  the  simulation  results  from  the  wave  and 
current  models,  respectively.  The  value  of  the  optimum  friction  coefficient 
was  identical  to  the  5  February  case  for  the  linear  model,  but  somewhat  higher 
for  the  nonlinear  model,  namely,  cf  —  0.0035,  with  the  mixing  coefficients  as 
determined  from  the  laboratory  tests.  The  true  optimum  with  respect  to  both 
cf  and  A  was  located  at  a  higher  A  value,  but  gave  little  improvement  in  the 
fit.  Smaller  values  of  the  mixing  parameter  gave  a  larger  sum  of  squares, 
with  the  optimum  occurring  for  a  larger  value  of  the  friction  coefficient. 
Discussion 

227.  For  the  Thornton  and  Guza  (1986)  data  set,  the  most  notable  dif¬ 
ference  between  model  predictions  and  measurements  was  not  the  shape  of  the 
current  profile,  but  the  strength  and  location  of  the  peak  current,  where  the 
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Wave  Height  (m) 


a.  Wave  height  distribution  across -shore 


Current  Velocity  (m/sec) 


b.  Longshore  current  profile 

Figure  21.  Simulations  for  Thornton  and  Guza  (1986)  data 

(4  February) 
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Wave  Height  (m) 


a.  Wave  height  distribution  across -shore 


Current  Velocity  (m/sec) 


b.  Longshore  current  profile 


Figure  23.  Simulations  for  Thornton  and  Guza  (1986)  data 

(6  February) 
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maximum  in  the  predicted  current  profile  for  all  four  cases  was  located 
shoreward  of  the  peak  in  the  measurements.  The  calculated  wave  height  across 
the  profile  was  in  reasonable  agreement  with  the  measurements.  The  agreement 
between  wave  model  and  measurements  could  be  improved  by  reducing  the  breaker 
index  to  about  0.4  to  0.5;  however,  this  reduction  is  physically  unrealistic 
because  yb  represents  the  breaker  index  for  a  monochromatic  wave  in  the 
present  context.  In  contrast  to  the  definition  of  yb  used  here,  Thornton  and 
Guza  (1983,  1986)  defined  a  quantity  that  might  be  denoted  as  7  (without  the 
subscript  b)  to  express  the  average  ratio  between  wave  height  and  water  depth 
of  all  waves,  breaking,  broken,  and  unbroken,  in  the  surf  zone.  Thus,  this 
ratio  is  a  kind  of  integrated  description  of  wave  conditions  after  the  broken 
wave  has  stabilized  in  the  form  of  a  bore,  in  some  sense  comparable  to  the 
stable  wave  height  concept  included  in  the  Dally  (1980)  breaker  decay  model, 
and  is  not  equivalent  to  the  breaker  index  of  an  individual  wave. 

228.  It  should  be  kept  in  mind  that  the  field  environment  may  have 
contributions  that  influence  the  current  measurements  in  ways  not  accounted 
for  in  the  present  models.  For  example,  the  tidal  current,  atmospheric 
pressure  differences,  small  changes  in  offshore  and  nearshore  bathymetry  and 
bottom  roughness,  and  limitations  in  the  measurements  are  all  unquantified. 

In  this  respect,  very  close  agreement  between  numerical  calculations  and  field 
measurements  should  not  be  sought  by  unrealistic  manipulation  of  model  parame¬ 
ters  . 

229.  One  way  of  improving  the  predictive  capability  of  the  longshore 
current  models  is  to  manipulate  parameters  in  the  breaker  decay  model.  A 
larger  wave  energy  dissipation  rate  would  produce  a  steeper  decay  in  wave 
height  and  translate  the  peak  of  the  wave  height  distribution  seaward,  and  the 
body  of  the  current  profile  would  accordingly  move  seaward.  The  effect  of 
moderate  changes  in  wave  model  parameters  was  not  considered  significant 
enough  to  warrant  calibration. 

230.  In  the  random  wave  selection  process  from  a  Rayleigh  distribution, 

a  certain  number  of  small  waves  that  break  close  to  shore  are  produced.  These 
waves  create  narrow  surf  zones  only  a  few,  perhaps  even  one,  calculation  cells 
wide.  In  the  longshore  current  model,  wave  breaking  very  close  to  shore  can 
create  a  sharp  peak  in  the  current  profile  and  unduly  influence  the  averaging 
of  all  current  distributions.  Trial  simulations  were  carried  out  where  wave 
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conditions  were  neglected  if  the  number  of  surf  zone  cells  was  less  than  a 
predetermined  value.  A  criterion  up  to  five  cells  was  used  (in  other  words,  a 
surf  zone  of  at  least  5-m  width  was  required) ,  which  marginally  improved  the 
fit  of  the  calculation  to  the  measurements.  Because  of  the  relatively  small 
effect  of  the  cutoff,  it  was  not  used  in  the  simulations,  and  the  calculations 
described  above  were  carried  out  with  all  waves  obtained  from  the  Rayleigh 
distribution.  However,  the  phenomenon  of  very  small  wave  heights  appearing  in 
the  Monte-Carlo  process  should  be  kept  in  mind  in  application  of  the  random 
wave  procedure . 

231.  In  summary  the  following  observations  were  made  concerning  model 
simulations  for  the  Thornton  and  Guza  (1989)  data  set: 

a.  The  rms  wave  height  variation  across -shore  was  in  general  well 
predicted,  with  the  greatest  discrepancy  in  the  region  of  steep 
decay. 

b.  The  overall  longshore  current  profile  was  correctly  reproduced 
by  the  models ,  but  the  predicted  location  of  the  maximum 
current  was  more  shoreward  than  in  the  measurements . 

c.  The  linear  and  nonlinear  models  gave  about  the  same  agreement 
in  a  least  squares  sense  in  calibration  against  the 
measurements . 

d.  The  average  value  of  the  optimum  friction  coefficient  for  the 
linear  model  was  cf  =  0.0050  and  cf  =  0.0040  for  the  nonlinear 
model  using  mixing  values  given  from  the  laboratory  simula¬ 
tions  . 

e.  The  true  optimum  with  for  free  variation  of  both  cf  and  A  was 
located  at  higher  A,  but  the  decrease  in  the  sum  of  squares  was 
small  as  compared  with  using  mixing  values  from  the  laboratory 
cases . 

f .  A  mixing  value  of  A  =  0  gave  a  smooth  velocity  profile  but 
increased  the  sum  of  squares.  The  optimum  value  of  the  fric¬ 
tion  coefficient  increased  if  A  =  0  was  employed,  giving  cf  = 
0.0065  for  the  linear  model  and  cf  =  0.0041  for  the  nonlinear 
model . 

232.  Compared  with  previous  independent  modeling  attempts,  the  present 
model  gives  agreement  similar  to  that  in  the  Thornton  and  Guza  (1986)  study, 
but  not  as  good  as  that  in  the  Wu,  Thornton,  and  Guza  (1985)  study.  The 
reason  for  the  good  agreement  obtained  by  Wu,  Thornton,  and  Guza  is  believed 
to  be  the  small  value  of  the  breaker  index  employed  (yb  —  0.44),  causing  the 
waves  to  break  farther  offshore.  In  the  present  work,  tests  indicated  that 
equally  good  fits  resulted  with  this  value  of  the  breaker  index,  but  such  a 
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small  value  of  yb  in  a  monochromatic  wave  model  is  not  justified.  In  Thornton 
and  Guza,  the  monochromatic  wave  approach  used  in  Wu,  Thornton,  and  Guza  was 
abandoned  and  a  randomized  model  presented.  The  random  wave  model  as  a 
driving  model  resulted  in  maxima  in  the  current  profiles  located  shoreward  of 
those  of  the  measured  profiles,  similar  to  the  present  results. 

233.  Values  of  friction  coefficients  obtained  with  the  present  models 
differ  from  those  of  Thornton  and  Guza  (1986),  being  on  the  average  about 
50  percent  smaller.  One  possible  explanation  may  lie  in  the  method  of 
ensemble  averaging  of  the  governing  momentum  equation.  Neglecting  lateral 
mixing  (which  would  tend  to  require  increased  values  of  the  bottom  friction 
coefficient) ,  Thornton  and  Guza  (1986)  averaged  the  multiplying  factor  in  the 
linearized  friction  term  (Equation  58),  which  is  a  function  of  wave  height, 
separately  from  the  driving  term.  If  instead,  the  driving  term  is  divided  by 
the  factor  that  multiplies  the  current  velocity  in  the  bottom  friction  term 
and  the  averaging  procedure  then  performed,  the  calculated  current  would  be  a 
factor  of  1.5  smaller.  The  latter  approach  is  considered  by  the  authors  of 
this  report  to  be  the  statistically  appropriate  method  and  results  if  the 
Monte-Carlo  simulation  is  employed.  Thus,  for  the  two  averaging  methods  to 
yield  the  same  current  magnitude,  the  friction  coefficient  in  the  present 
model  should  be  taken  as  approximately  two -thirds  that  obtained  by  Thornton 
and  Guza  (1986)  (for  the  simple  case  of  no  mixing  and  use  of  a  linearized 
friction  term) .  Inclusion  of  lateral  mixing  or  use  of  a  nonlinear  friction 
term  precludes  specification  of  a  simple  factor  connecting  the  two  methods, 
but  does  not  change  the  qualitative  relation  between  the  methods. 


jothetical  Case  from  Ebersole  and  Dalr\ 


234.  Because  the  available  data  sets  for  which  the  current  models  were 
®v&luated  encompassed  only  beach  profiles  with  monotonically  increasing  depth 
with  distance  offshore  or  a  step- type  beach,  a  hypothetical  case  was  chosen  to 
qualitatively  investigate  the  performance  of  the  model  on  more  complex  topo¬ 
graphies.  A  simple  bar-type  profile  was  selected  as  presented  by  Ebersole  and 
Dalrymple  (1980) ,  and  the  same  deepwater  wave  conditions  were  employed:  wave 
height  of  2.0  m,  period  of  8.0  sec,  and  incident  angle  of  30  deg.  The  breaker 
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index  was  set  to  yb  =  0.78,  and  the  calculation  grid  extended  out  to  7-m  water 
depth.  The  depth  values  for  the  bottom  topography  are  given  in  Appendix  B. 

235.  The  incident  monochromatic  waves  broke  on  the  bar,  reformed  in  the 
trough  shoreward  of  the  bar,  and  broke  a  second  time  farther  inshore. 

Figure  24  illustrates  changes  in  the  longshore  current  distribution  over  the 


Current  Velocity  (m/sec) 


Figure  24.  Influence  of  lateral  mixing  for  a  barred  profile 

barred  profile  for  different  values  of  the  lateral  mixing  coefficient  in  the 
linear  model.  A  friction  coefficient  of  cf  =  0.01  was  specified  to  agree  with 
Ebersole  and  Dalrymple  (1980) .  For  the  typical  wave  height  distribution  over 
a  barred  profile,  reference  is  made  to  Part  III  where  the  wave  model  was 
investigated.  An  increase  in  the  lateral  mixing  coefficient  lowers  the  peak 
in  the  current  profile  and  elongates  the  tail.  For  moderate  changes  in  A,  the 
difference  is  small,  not  being  notable  until  very  low  values  of  the  mixing 
coefficient  are  reached.  For  the  case  of  no  mixing  (A  —  0) ,  two  distinct 
peaks  result  with  the  outer  one  displaying  a  larger  maximum  current  velocity. 
The  simulation  with  the  limiting  case  A  =  0  with  monochromatic  waves,  which 
yields  an  almost  discontinuous  current  profile,  also  demonstrates  the  numeri¬ 
cal  stability  of  the  model. 
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236.  Figure  25  shows  the  difference  in  predictions  between  the  linear 
and  nonlinear  models  for  the  same  value  of  the  friction  and  lateral  mixing 
coefficients  (cf  =  0.01,  A  =  0.5).  As  expected,  the  quadratic  friction  term 
lowers  the  velocity  profile.  Apart  from  the  difference  in  magnitude  between 
the  two  distributions,  the  general  shapes  are  similar.  However,  the  effect  of 
the  nonlinear  stress  term  on  the  velocity  profile  is  greater  for  the  larger 
current  velocities  because  of  quadratic  dependence  on  current  magnitude.  In 
deeper  water,  as  the  velocity  goes  to  zero,  the  linear  and  nonlinear  model 
calculations  gradually  converge. 

Current  Velocity  (m/sec) 


Figure  25.  Model  simulations  for  a  barred  profile 


Longshore  Current  Driven  by  Wind 


Longshore  current 

237 .  To  examine  the  direct  influence  of  wind  on  the  longshore  current 
profile  and  mean  water  surface  elevation  in  the  model,  several  hypothetical 
cases  were  simulated.  Two  different  basic  profile  shapes  were  used,  a  plane- 
sloping  beach  and  a  barred  beach.  The  barred  profile  and  incident  wave 
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conditions  (deepwater  wave  height  Ha  =  2.0  m,  T  =  8.0  sec,  and  deepwater  wave 
angle  8a  =  30  deg)  were  taken  from  Ebersole  and  Dalrymple  (1980),  and  the 
plane  beach  was  given  a  slope  tanS  =  0.024.  The  wave  model  was  applied  with 
the  same  parameter  values  used  in  previous  simulations,  and  the  linear  and 
nonlinear  current  models  were  applied  with  (c£  =  0.0050,  A  =  0.50),  and 
(cf  =  0.0040,  A  =  0.30),  respectively.  The  smaller  value  of  the  friction 
coefficient  in  the  nonlinear  model  produced  current  velocities  comparable  to 
the  linear  model,  facilitating  comparison  of  model  predictions. 

238.  A  wind  speed  of  15  m/sec  was  used,  and  the  incident  angle  of  the 
wind  4>  was  taken  as  +30  and  -30  deg  to  obtain  conditions  that  reinforced  and 
opposed  the  wave- induced  longshore  current  (negative  incident  wave  angle), 
respectively.  Different  ratios  of  CD/e£  were  investigated,  and  monochromatic 
wave  conditions  were  employed.  Figures  26a  and  b  show  the  influence  of  wind 
forcing  on  the  calculated  longshore  current  distribution  for  the  plane -sloping 
beach  as  predicted  by  the  linear  and  nonlinear  current  models .  The  distribu¬ 
tion  occurring  without  wind  is  also  shown  for  comparison.  In  the  figures,  the 
current  distributions  for  CD/c£  ratios  of  1.0  and  2.0  are  shown  for  selected 
incident  wind  angles. 

239.  A  wind  with  an  alongshore  component  in  the  same  direction  as  the 
current  induced  by  breaking  waves  simply  produces  an  overall  increase  in 
current  speed.  The  linear  model  predicts  a  change  in  current  similar  in 
magnitude  across  the  shore,  whereas  the  nonlinear  model  shows  relatively 
smaller  changes  in  the  vicinity  of  maximum  in  the  distribution,  a  result  of 
the  quadratic  dependence  of  the  velocity  in  the  friction  term.  Similar 
behavior  is  seen  for  the  opposing  wind  and  current  situation,  with  the 
nonlinear  model  showing  less  reduction  in  current  velocity  around  the  peak  in 
the  distribution.  If  the  opposing  wind  is  sufficiently  strong,  it  may 
dominate  the  wave -induced  contribution  and  reverse  the  current.  A  change  in 
current  direction  in  the  opposing  wind  and  wave  situation  is  particularly 
likely  at  the  shoreward  and  seaward  ends  of  the  current  where  the  wave -induced 
current  is  weak,  as  shown  in  Figures  26a  and  b  for  some  of  the  simulations. 

240.  As  is  discussed  in  the  following  paragraphs,  the  modification  of 
the  boundary  conditions  to  include  wind  results  in  a  slightly  increasing 
current  velocity  with  distance  offshore  for  the  linear  model,  whereas  the 
nonlinear  model  produces  a  constant  current  offshore.  This  trend  is  seen  in 
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a.  Linear  model 


-50  0  50  100  150  200  250  300 

Distance  Offshore  (m) 


b.  Nonlinear  model 

Figure  26.  Effect  of  wind  on  the  longshore  current  (plane  beach) 
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Figures  26a  and  b  in  the  seaward  tail  of  the  profiles  where  the  linear  model 
result  displays  a  larger  gradient  than  the  nonlinear  model.  The  current 
increases  almost  in  proportion  to  the  ratio  CD/cf  in  the  linear  model  with  a 
relatively  larger  increase  at  the  ends  of  the  distributions .  For  the  non¬ 
linear  model,  the  dependence  of  the  friction  term  on  the  velocity  squared 
makes  the  effect  of  an  increase  in  CD/cf  somewhat  nonuniform. 

241.  The  same  simulations  as  for  the  plane-sloping  beach  were  carried 
out  for  the  barred  profile  introduced  by  Ebersole  and  Dalrymple  (1980) . 

Figures  27a  and  b  show  the  simulation  results  for  the  linear  and  nonlinear 
models,  respectively.  The  influence  of  the  wind  on  the  current  distribution 
exhibits  the  same  general  characteristics  found  for  the  plane-sloping  beach. 
However,  the  current  distribution  is  bimodal  because  the  bars  produce  multiple 
wave  breaking  and  subsequent  reformation.  The  trough  behind  the  bar  reduces 
the  wave  energy  dissipation  and  thus  the  change  in  the  alongshore -directed 
radiation  stress  driving  the  wave -induced  longshore  current.  Because  of  the 
change  in  radiation  stress  gradient,  a  minimum  in  the  current  profile  is 
located  just  shoreward  of  the  bar,  and  a  second  peak  is  produced  closer  to 
shore  as  the  energy  dissipation  again  becomes  larger.  The  relatively  large 
wave- induced  contribution  nearshore  makes  it  difficult  for  the  opposing  wind 
to  reverse  the  current  near  the  shoreline.  For  the  nonlinear  model,  the  wind 
effect  on  the  current  profile  is  relatively  small  around  the  peaks,  but 
becomes  important  in  the  trough  region  and  seaward  of  the  breaker  zone .  The 
larger  valued  drag  coefficient  giving  CD/cf  =  2.0  reverses  the  current  in  the 
trough  area  for  the  opposing  wind. 

Mean  water  surface  elevation 

242.  The  wind  not  only  enters  as  a  forcing  term  in  the  longshore 
momentum  equation,  but  also  contributes  to  the  cross -shore  momentum  equation. 
For  example,  a  wind  blowing  onshore  will  tilt  the  water  surface  shoreward  to 
balance  the  stress  at  the  air- sea  interface.  Wind  setup  is  a  menacing 
phenomenon  during  storms  as  it  allows  waves  to  attack  high  on  the  beach, 
leading  to  severe  erosion  and  inundation.  In  a  numerical  model,  the  effect  of 
the  wind  on  the  mean  water  surface  elevation  will  be  underestimated  unless  the 
grid  extends  seaward  much  farther  than  the  distance  typically  covered  in  a 
longshore  current  calculation,  and  the  model  allows  opening  of  "wetted"  cells 
on  the  previously  dry  beach  (as  allowed  in  NMLONG) . 
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a.  Linear  model 


b.  Nonlinear  model 

Figure  27 .  Effect  of  wind  on  the  longshore  current  (barred  profile) 
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243.  To  estimate  the  importance  of  wind- induced  setup,  the  simple  case 
of  an  onshore -blowing  wind  without  waves  is  examined.  If  the  wind  setup  is 
small  in  comparison  with  the  water  depth,  in  the  absence  of  waves  the  momentum 
equation  (Equation  54)  may  be  written, 
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Equation  90  expresses  a  force  balance  between  the  surface  shear  stress  and  the 
hydrostatic  pressure  caused  by  the  tilt  of  the  water  surface.  For  a  wind 
blowing  onshore,  the  setup  decreases  going  seaward,  according  to  the  defini¬ 
tion  of  sign  given  in  Figure  1.  For  simple  bottom  topographies,  Equation  90 
is  readily  integrated  to  yield  an  analytic  solution  for  rj .  In  the  case  of  a 
plane  beach  of  bottom  slope  tan) 3,  Equation  90  gives, 
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where  xQ  denotes  the  seaward  location  where  the  wind  stress  starts  to  act. 
Equation  91  can  only  be  applied  to  within  a  certain  distance  from  the  shore¬ 
line  since  at  the  shoreline  h  equals  zero,  giving  infinite  setup.  This  is 
simply  a  consequence  of  assuming  that  rj  is  negligible  in  comparison  with  h;  at 
the  shoreline  this  assumption  is  inapplicable. 

244.  In  the  previous  case  of  a  plane -sloping  beach  and  a  moderate 
onshore -blowing  wind  with  W  =  15  m/sec  and  4>  —  30  deg,  wind- induced  setup  over 
a  500-m-long  grid  was  calculated  by  Equation  91  to  be  only  about  3.0  cm  at  a 
point  1  m  from  the  shoreline.  In  the  example  CD  =  5  10’3  was  used.  To 
reproduce  realistic  wind  setup,  it  is  necessary  to  describe  wind  conditions 
over  a  longer  stretch  than  what  is  normally  used  in  calculating  wave -induced 
longshore  currents . 

Numerical  solution 

245.  Numerical  solution  of  the  momentum  equations  incorporating  wind 
driving  requires  no  special  techniques ;  only  values  of  the  drag  coefficient 
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and  wind  speed  and  direction  must  be  specified.  However,  the  boundary  condi¬ 
tions  must  be  modified  since  a  wind  stress  is  now  acting  at  these  two  points. 
Boundary  conditions  are  obtained  by  explicitly  solving  for  the  current 
velocity  without  lateral  mixing.  Thus,  the  bottom  friction  force  is  equated 
to  the  driving  force,  including  wind.  At  the  seaward  boundary,  the  wind- 
induced  current  calculated  by  the  linear  model  can  take  an  unreasonably  high 
value  if  the  grid  is  located  in  deep  water  compared  with  the  wavelength, 
because  the  current  is  inversely  proportional  to  the  wave  orbital  velocity  at 
the  bottom,  which  approaches  zero.  For  a  nonlinear  friction  term,  the  bottom 
resistance  approaches  the  familiar  velocity- squared  dependence  if  the  orbital 
velocity  is  small. 

246.  Figure  28  shows  calculated  results  for  the  aforementioned  wind  and 
beach  conditions  for  onshore  and  offshore  winds .  The  numerical  results  agree 
well  with  the  simple  analytic  model  solution  (simple  because  the  total  water 
depth  is  approximated  by  the  still -water  depth) .  The  numerical  solution  is 
asymmetric  because  the  total  water  depth  is  used. 

Wind  Setup/Setdown  (m) 


Distance  Offshore  (m) 

Figure  28.  Simulation  of  wind-induced  setup  and  setdown 
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Wave -Current  Interaction 


247.  If  the  wave-current  interaction  is  activated  in  the  model,  the  wave 
height  and  current  distributions  cannot  be  obtained  independently,  and  an 
iterative  technique  must  be  employed.  Wave-current  interaction  terms  appear 
in  the  wave  energy  flux  equation  (Equation  7)  and  wave  refraction  equation 
(Equation  19).  Selected  results  are  described  here  from  the  calculations  to 
quantify  the  importance  of  wave -current  interaction  and  how  it  influences  the 
wave  height  and  current  profiles. 

248 .  Both  a  plane  beach  and  the  barred  profile  of  Ebersole  and  Dalrymple 
(1980)  were  used  to  evaluate  the  influence  of  the  wave -current  interaction  on 
the  profile.  Empirical  parameters  in  the  wave  model  were  set  to  default 
values,  and  the  current  was  computed  with  the  nonlinear  model  with  cf  =  0.01 
and  A  =  0.5  as  in  some  of  the  previous  simulations  for  the  hypothetical  barred 
profile  (Figure  25) .  The  plane-sloping  beach  was  identical  to  the  previous 
case.  A  breaker  index  of  0.78  was  employed  in  both  examples,  following 
Ebersole  and  Dalrymple. 

249 .  The  current  velocity  was  set  to  zero  to  calculate  the  wave  height 
at  the  first  iteration,  and  subsequent  calculations  of  the  wave  height 
included  the  current  velocity  determined  at  the  previous  iteration.  This 
procedure  continued  until  two  consecutively  calculated  distributions  of  the 
wave  height  and  current  respectively  differed  by  less  than  1  percent  at  every 
grid  point.  The  calculation  normally  converged  in  two  to  three  iterations, 
the  small  number  of  iterations  indicating  that  the  wave -current  interaction 
was  weak. 

Plane -sloping  beach  (Figure  29) 

250.  The  difference  in  wave  height  incorporating  the  wave-current 
interaction  is  noticeable  only  near  the  break  point,  where  the  wave  height  is 
increased  slightly.  The  effect  of  the  interaction  on  the  longshore  current  is 
somewhat  greater,  but  still  weak,  being  most  pronounced  around  the  maximum  in 
the  current  and  extending  both  shoreward  and  seaward  of  the  maximum.  The 
wave-current  interaction  raises  the  tail  of  the  current  profile,  lowers  the 
current  velocity  at  the  maximum  while  translating  the  location  of  the  maximum 
current  seaward,  and  slightly  decreasing  the  current  magnitude  in  the  surf 
zone . 
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Barred  profile  (Figure  30) 

251.  In  this  example,  wave  breaking  occurs  on  the  bar  with  reformation 
where  the  energy  flux  falls  below  the  stable  energy  flux  in  the  deeper  part  of 
the  trough.  The  wave -current  interaction  increases  the  wave  height  slightly 
around  the  break  points,  similar  to  the  plane -sloping  beach  case,  and  most 
notably  at  the  seaward  break  point  where  the  current  velocity  is  greatest. 

The  magnitude  of  the  current  is  decreased  slightly  in  the  vicinity  of  the  two 
maxima,  and  changes  negligibly  in  the  trough  behind  the  bars. 

Separate  effects  (Figure  31) 

252.  The  wave-current  interaction  enters  in  two  ways  in  NMLONG:  through 
the  wave  energy  flux  equation  (Equation  7)  as  an  interaction  term  with  the 
radiation  stress  and  through  wave  refraction  (current  refraction) 

(Equation  19).  To  quantify  the  influence  of  these  terms,  their  individual  and 
combined  effects  were  studied  by  example  with  the  previous  plane -sloping  beach 
case.  Figure  31a  shows  that  refraction  negligibly  modifies  the  wave  height 
distribution  and  can  for  all  practical  purposes  be  neglected.  This  conclusion 
holds  only  for  the  wave -current  interaction  for  an  alongshore  system;  even 
casual  observation  shows  that  localized  strong  seaward  flow  in  rip  current 
throats  can  greatly  change  wave  properties  and  even  prematurely  break  the 
waves  before  they  become  depth  limited. 

253.  Calculations  showed  that  modifications  of  the  current  by  the 
refraction  and  radiation  stress  terms  were  of  the  same  order  of  magnitude  for 
the  test  cases  (see  Figure  31b) .  In  addition,  the  two  terms  counteract  each 
other,  the  radiation  stress  term  increasing  the  current  around  the  peak  and 
decreasing  the  current  along  the  tail,  and  the  refraction  term  acting  in  the 
opposite  way. 

Discussion 

254.  The  calculated  examples  demonstrated  that  the  wave -current  interac¬ 
tion  has  relatively  minor  Influence  on  the  wave  height  and  longshore  current 
profiles  for  the  situation  of  longshore  uniformity  with  no  strong  offshore - 
directed  flows  such  as  would  occur  at  a  rip  current.  The  examples  pertained 
to  current  velocities  reaching  slightly  more  than  1  m/sec,  which  is  consider¬ 
ably  greater  than  the  current  on  most  beaches  under  average  wave  conditions. 
Storm  waves  and  waves  approaching  at  unusually  large  incident  angles ,  as  well 
as  high  winds,  could  produce  significantly  greater  current  velocities.  For 
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b.  Longshore  current 

Figure  30.  Wave-current  interaction  on  a  barred  beach 
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Wave  Height  (m) 


a.  Wave  height 


Current  Velocity  (m/sec) 


b .  Longshore  current 

Figure  31.  Evaluation  of  the  wave-current  interaction 
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example,  longshore  current  velocities  recorded  at  the  FRF  have  reached  the 
range  of  3  to  5  m/sec  in  large  storms.  In  these  cases,  inclusion  of  the  wave- 
current  interaction  may  be  necessary.  Otherwise,  the  increased  calculation 
time  required  for  solution  with  the  wave -current  interaction  does  not  appear 
warranted  in  light  of  the  predictive  accuracy  of  a  model  driven  by  engineering 
data.  The  weakness  of  the  wave-current  interaction  terms  in  the  wave  energy 
flux  and  longshore  current  governing  equations,  combined  with  the  fact  that 
the  wave  height  decay  depends  on  two  empirical  parameters,  makes  it  doubtful 
whether  these  particular  interaction  terms  need  ever  be  activated.  The 
standard  implementation  of  NMLONG  does  not  iterate  between  the  wave  and 
current  models. 
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PART  VII:  SUMMARY  AND  CONCLUSIONS 


Modeling  System  Capabilities 

255.  A  numerical  modeling  system  called  NMLONG  was  developed  that  has 
the  capability  to  calculate  the  wave  height,  wave  direction,  mean  water  level, 
and  longshore  current  across  a  multiple  bar  and  trough  beach  profile.  The 
calculated  longshore  current  is  time -averaged  and  depth- integrated.  A  major 
practical  restriction  on  the  modeling  system  is  the  requirement  of  longshore 
uniformity,  meaning  that  the  beach  topography,  water  level,  and  wave  condi¬ 
tions  should  be  uniform  alongshore . 

256.  Major  features  of  NMLONG  may  be  summarized  as  follows: 

a.  Efficient  numerical  solution  of  both  the  wave  and  longshore 
current  equations,  while  maintaining  general  forms  of  linear- 
wave  theory  expressions  for  arbitrary  depth. 

b.  Capability  to  calculate  the  wave  height,  wave  direction,  and 
longshore  current  over  a  beach  profile  consisting  of  multiple 
longshore  bars  and  troughs . 

c.  Incorporation  of  either  linear  or  quadratic  (nonlinear)  bottom 
friction  stresses  as  options. 

d.  Incorporation  of  wave  and  wind  driving  for  both  the  longshore 
current  and  mean  water  level. 

e.  Incorporation  of  a  rational  formulation  for  lateral  mixing. 

f.  Wide-ranging  verification  and  sensitivity  testing. 

g.  Extension  of  the  regular-wave  model  to  random  wave  height  as  an 
option,  and  calculation  of  the  associated  longshore  current. 

h.  Implementation  on  a  desk- top  computer  (not  described  in  the 
present  report) . 

257.  The  three  major  calculation  components  of  the  modeling  system  were 
verified  separately.  The  wave  height  calculation  was  verified  with  data  from 
small-scale  laboratory  experiments  involving  regular  waves  incident  at  an 
angle  to  the  beach,  prototype-scale  tank  data  involving  regular  wave  transfor¬ 
mation  over  an  irregular  bottom,  and  with  high-quality  field  data  for  which  a 
Monte-Carlo  random  wave  calculation  procedure  was  successful.  The  mean  water 
level  calculation  was  verified  with  laboratory  measurements,  and  the  longshore 
current  calculation  was  verified  with  laboratory  measurements  and  two  field 
data  sets. 
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258.  Sensitivity  tests  were  performed  on  all  model  components,  partic¬ 
ularly  those  components  for  which  verification  data  were  not  available  such  as 
the  wind- induced  setup  and  setdown  and  the  wind- induced  longshore  current. 
Tests  were  also  made  of  components  operating  in  combination,  such  as  longshore 
current  generation  by  combined  wind  and  waves.  These  tests  gave  physically 
reasonable  results  and  also  verified  the  stability  of  the  numerical  solution 
algorithms  in  examination  of  general  as  well  as  limiting  cases. 

259.  Calculation  efficiency  was  an  important  design  requirement  for  the 
modeling  system,  because  a  major  goal  of  this  study  is  implementation  of 
NMLONG  on  desk- top  computers.  Calculation  speed  of  the  combined  wave  and  mean 
water  level  calculation  was  optimized  by  adopting  the  procedure  of  Dally 
(1980)  for  estimating  the  mean  water  level,  thereby  eliminating  the  necessity 
of  iteration  between  the  waves  and  water  level.  An  explicit  numerical 
solution  scheme  for  the  wave  height  and  water  level  is  more  efficient  and 
convenient  than  an  implicit  scheme,  because  the  shoreward  boundaries  cannot  be 
precalculated,  being  controlled  by  wave  breaking  and  reformation  criteria 
which  must  be  evaluated  cell  by  cell.  The  profile  of  the  mean  longshore 
current  is  obtained  by  an  implicit  double -sweep  solution  procedure.  Although 
such  algorithms  are  commonly  used  in  large  2-D  numerical  circulation  models, 
this  appears  to  be  the  first  application  to  the  simple  1-D  case.  The  double¬ 
sweep  solution  technique  is  stable,  accurate,  and  much  superior  in  efficiency 
to  time -stepping,  relaxation,  and  "shooting"  solution  procedures  that  have 
been  used  in  other  numerical  models. 

260.  Execution  time  of  the  linear  model  for  a  barred  profile  discretized 
with  150  grid  points  including  realistic  wave  and  wind  forcing  is  approxi¬ 
mately  1  sec  on  a  20 -MHz  desk- top  computer.  This  solution  yields  the  wave 
height,  water  level,  and  longshore  current  velocity  at  all  grid  points. 
Execution  time  of  the  nonlinear  friction  model  for  the  same  condition  is 
approximately  5  sec  and  is  related  to  the  product  of  the  number  of  iterations 
and  the  execution  time  of  the  linear  model  solution,  with  great  computational 
cost  savings  made  by  implementation  of  the  Nishimura  (1982)  square -wave 
approximation.  By  use  of  a  correction  factor  for  estimating  the  nonlinear 
solution  based  on  the  linear  solution,  typically  only  three  to  five  iterations 
are  required  to  obtain  a  solution  with  the  nonlinear  friction  force. 
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Fundamental  Processes 


261.  Although  the  major  objective  of  this  study  was  development  of  a 
general  pragmatic  model  of  the  longshore  current  for  engineering  use,  several 
fundamental  governing  processes  had  to  be  investigated. 

Nonlinear  bottom  friction 

262.  The  finding  of  Thornton  and  Guza  (1986)  was  confirmed  that  a 
quadratic  friction  law  acts  mainly  to  reduce  the  value  of  the  average  bottom 
friction  coefficient  compared  with  the  value  obtained  with  a  linear  friction 
law.  A  correction  factor  was  developed  to  estimate  the  nonlinear  current 
velocity  by  multiplication  of  the  linear  solution.  Alternatively,  the 
correction  factor  can  be  interpreted  as  relating  the  friction  coefficients  in 
the  nonlinear  and  linear  solutions.  An  approximation  for  this  factor  (Equa¬ 
tion  A13)  was  obtained  as  a  function  of  the  wave  angle  and  ratio  of  current 
velocity  of  the  linear  solution  and  magnitude  of  the  wave  orbital  velocity. 

263.  The  Nishimura  (1982)  square-wave  approximation  of  the  quadratic 
friction  force  was  investigated  and  found  to  be  an  accurate  and  efficient 
method  for  replacing  computation- intensive  time -integration  of  the  friction 
term.  This  technique  provides  considerable  savings  in  computation  time, 
although  iteration  is  still  required  in  obtaining  a  solution  of  the  nonlinear 
current  model . 

264.  Generally  speaking,  the  linear  current  model  provided  as  good  or 
superior  agreement  with  measurements  of  the  current  in  the  surf  zone  as  the 
nonlinear  model,  with  the  nonlinear  model  giving  better  agreement  seaward  of 
the  breaker  line  because  the  calculated  current  decreased  more  steeply.  For 
engineering  use,  at  the  present  stage  of  understanding  of  bottom  friction 
factors  and  sediment  transport  in  the  surf  zone,  the  linear  model  is  consid¬ 
ered  sufficiently  accurate. 

Wave -current  interaction 

265.  Numerical  experimentation  confirmed  previous  results  (Dalrymple 
1980,  Thornton  and  Guza  1986)  that  the  wave-current  action  is,  for  practical 
purposes,  negligible  for  an  alongshore  current  system  except,  perhaps,  for 
extreme  conditions  of  very  large  wave  angles  and  wave  heights.  The  radiation 
stress  interaction  term  and  the  current-refraction  term  are  of  the  same 
magnitude  and  act  to  cancel  each  other. 
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Lateral  mixing 

266 .  It  was  found  necessary  to  include  lateral  mixing  in  the  current 
models  to  reproduce  current  profiles  measured  in  a  carefully  controlled 
laboratory  experiment  using  monochromatic  regular  waves  as  well  as  profiles 
measured  in  the  field  produced  by  random  waves.  The  value  of  the  empirical 
mixing  coefficient  was  approximately  the  same  for  the  laboratory  and  field 
tests.  A  rational  form  of  the  eddy  viscosity  coefficient  was  introduced 
(Equation  48)  and  appears  to  be  more  applicable  to  barred  profiles  than  other 
commonly  used  expressions. 

Random  waves 

267.  A  Monte-Carlo  approach  based  on  a  Raleigh  distribution  was  develop¬ 
ed  to  give  the  wave  height  probability  distribution  of  random  waves  transform¬ 
ing  by  shoaling,  refraction,  breaking,  and  reformation  at  intervals  across  a 
surf  zone  with  a  bar  and  trough  bottom  profile.  In  this  approach,  between  100 
and  500  waves  are  randomly  selected  from  a  Rayleigh  distribution  defined  by 
the  rms  wave  height  and  then  transformed  individually.  Statistical  properties 
are  then  obtained  by  invoking  superposition.  In  this  procedure,  a  true 
breaker  index  corresponding  to  regular  monochromatic  waves  is  employed. 
Calculation  results  were  compared  with  high-quality  measurements,  and  good 
agreement  was  found.  The  Monte -Carlo  random  wave  model  drove  the  longshore 
current  model,  and  favorable  results  were  obtained  in  reproducing  current 
measurements  in  five  field  data  sets. 

268.  The  Monte-Carlo  approach  is  computation  intensive  because  on  the 
order  of  100  model  runs  must  be  made  to  obtain  a  reasonable  representation  of 
the  wave  height  probability  distribution  and  resultant  current.  The  advantage 
is  that  the  method  is  applicable  to  a  bar  and  trough  profile  without  invoking 
a  priori  assumptions  about  the  shape  of  the  transformed  probability 
distribution. 

Wind  drag  coefficient 

269.  The  order  of  magnitude  of  the  wind  drag  coefficient  was  estimated 
to  be  the  same  as  the  bottom  friction  coefficient  in  a  crude  empirical 
analysis  correlating  measured  winds  and  longshore  current  at  one  site. 

Several  formulas  for  the  drag  coefficient  were  reviewed  for  use  in  the 
longshore  current  model,  and  the  WAMDI  Group  (1988)  formula  was  selected  as 
being  most  appropriate  for  the  nearshore  zone. 


133 


Future  Work 


270.  The  next  stage  in  this  study  program  will  transfer  all  model 
developments  to  the  desk- top  computer  environment.  The  existing  prototype 
desk- top  model  will  be  completed  and  integrated  with  a  convenient  user  inter¬ 
face  that  includes  structured  data  input,  automatic  data  cleaning  and  error 
messages,  and  graphical  output.  The  results  will  be  presented  in  a  user's 
manual  to  be  distributed  with  the  model. 

271.  A  research  effort  is  required  to  understand  the  role  of  the  wind  in 
generating  the  longshore  current  and  changing  the  mean  water  surface  level . 

Of  particular  importance  is  how  the  wind  changes  the  location  of  the  break 
point  and  type  of  breaking,  processes  that  are  not  represented  in  existing 
models,  including  NMLONG. 

272.  The  availability  of  the  present  general  model  of  waves,  mean  water 
level,  and  longshore  current  in  the  surf  zone  should  stimulate  further  field 
data  collection.  The  field  program  should  include  synoptic  measurements  of 
the  profile  topography,  wave  height  and  direction,  wind  speed  and  direction, 
mean  water  level ,  and  longshore  current  for  a  wide  range  of  conditions .  One 
possible  use  of  a  well -verified  model  would  be  to  infer  information  on  the 
form  of  lateral  mixing  eddy  viscosity  coefficient,  which  is  difficult  to 
obtain  by  direct  measurement. 
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APPENDIX  A:  NISHIMURA'S  SQUARE-WAVE  APPROXIMATION 


1.  Nishimura  (1982,  1988)  introduced  an  approximate  method  for  numeri¬ 
cally  evaluating  the  time  integration  of  the  full  bottom  friction  term  in 
Equation  57  (see  main  text) .  In  this  method,  the  sinusoidally  varying 
horizontal  component  of  the  wave  orbital  velocity  u0  =  umc os (2 nt/T)  is 
replaced  by  a  square  wave  having  the  same  area  as  the  original  sinusoid.  The 
velocity  is  then  expressed  as: 


u 


O 


0  <  t  <  — ,  —T  <  t  <T 
4  4 


T 

4 


<  t  < 


(Al) 


Equation  Al  represents  the  least  squares  fit  of  a  square  wave  to  the  original 
cos (2nt/T)  term  in  Equation  57.  By  this  approximation,  the  bottom  orbital 
velocity  becomes  time  independent  in  the  different  intervals  indicated  above, 
and  the  integration  is  readily  carried  out  for  each  interval.  The  time 
integration  may  thus  be  written, 


f  =  — 

Lby  T 
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where 


Jby 

P 


cf\ju20  +  V2  +  2uoVsin0  ( V  +  uosin0) 


(A2) 


(A3) 


2.  To  evaluate  Equation  A2 ,  the  expression  for  uD  in  Equation  A3  should 
be  evaluated  for  the  respective  time  interval  according  to  Equation  Al.  Since 
the  integrand  is  constant  within  each  interval,  the  three  integrals  in 
Equation  A2  may  be  evaluated  as : 


Al 


(A4) 


—  f  4  tb  dt  =  ct  y/w2  +V2  +  2wksin0  (V  +  wsin0)  — 
p  J  o  4 


—  fT4  dfc  =  cf  y/w2  +V2  -  2wksin0  ( V  -  wsin0)  —  (A5) 

P  4  2 


—  f2T  tiy  dt  =  cf  y/w2  +V2  +  2wksin0  (V  +  wsin0)  —  (A6) 

P  J  ~T  4 


The  definition  w  =  (2/?r)um  was  introduced  in  Equations  A4  to  A6  following  the 
notation  of  Nishimura  (1982,  1988).  Thus,  Equation  A2  may  be  written, 


'by  =  ~~  \ZX  (V  +  wsin0)  +  Z2  ( V  -  wsin0)] 


(A7) 


where 

Zx  =  y/w2  +  V2  +  2wksin0  (A8) 


Z2  =  y/w2  +  V2  -  2wVsin0  (A9) 


3.  Rearrangement  of  Equation  A7  yields: 


fby  =  -y[V(Zi  +  Z2)  +  wsin0(Zx  -  Z2)] 


(A10) 


The  conjugate  expression  (Z1  +  Z2 )  is  used  to  multiply  (and  divide)  the  last 
term  within  the  bracket  of  Equation  A10,  and,  after  carrying  out  the  multi¬ 
plication,  the  following  simplified  equation  is  obtained: 


A2 


fby  =  CfV^Z  +  -^sin2ej 


(All) 


where 


Z  = 


Z1  +  ^2 
2 


(A12) 


4.  Equation  All  is  the  square -wave  approximation  for  the  time -averaged 
friction  term  as  presented  by  Nishimura  (1982,  1988)  if  no  cross-shore  current 
is  present.  It  can  be  shown  that  Equation  All  reduces  to  the  formal  expres¬ 
sions  for  the  linear  bottom  friction  term  for  weak  mean  current  (Equation  59) 
and  a  pure  quadratic  form  for  strong  mean  current  (Equation  60) . 

5.  In  the  nonlinear  current  model  developed  in  this  report,  Equation  All 
is  used  to  evaluate  the  bottom  friction  term.  Nishimura  (1988)  states  that 
Equation  All  deviates  at  most  by  10  percent  from  the  original  formulation  in 
Equation  57.  Figure  Al  displays  a  comparison  between  the  friction  term  as 
predicted  by  Equation  All  and  the  full  equation  (Equation  57)  for  various 
incident  wave  angles  and  ratios  V/um.  For  the  range  of  values  studied,  the 
approximation  is  quite  good  and  within  Nishimura' s  error  margin. 

6.  In  numerical  solution  of  the  nonlinear  longshore  momentum  equation 
for  the  current,  an  iterative  technique  is  employed.  To  evaluate  the  friction 
term  fby  the  velocity  obtained  at  the  previous  iteration  is  used  to  determine 
Z,  thereby  giving  a  linear  dependence  on  V.  Once  the  difference  between  two 
consecutively  calculated  velocity  distributions  is  smaller  than  some  predeter¬ 
mined  tolerance  at  each  grid  point,  the  calculation  is  halted.  As  an  initial 
guess,  the  linear  solution  with  a  friction  term  according  to  Equation  59  is 
used.  However,  depending  on  6  and  V/um,  this  initial  estimate  could  still  be 
far  off,  not  reducing  the  large  number  of  iterations  required.  Because 
reduction  of  execution  time  is  one  of  the  practical  objectives  of  this  study, 
improvements  in  the  initial  guess  are  of  importance. 

7.  As  discussed  in  Part  V,  if  lateral  mixing  is  neglected  it  is  possible 
to  derive  a  correction  factor  to  multiply  the  linear  solution  to  obtain  the 
nonlinear  solution.  This  factor  is  a  function  of  9  and  V/um  and  cannot  be 
expressed  analytically,  only  graphically.  To  improve  the  initial  guess  for 
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Figure  Al.  Evaluation  of  the  square-wave  approximation 

the  current  distribution,  a  correction  factor  was  derived  based  on  Equation 
All  and  Equation  59  (used  in  the  linear  solution) .  This  correction  factor  was 
then  approximated  with  the  following  empirical  equation  derived  for  the 
condition  V/lin/w  >0.5, 

a  =  (0.754  +  0.00314  |6| )  |  i^gl j  (A13) 

where  9  is  given  in  degrees  and  Flin  refers  to  the  linear  current  solution 
using  a  friction  term  according  to  Equation  59.  Thus,  if  Vlin/w  >0.5,  at  the 
first  iteration  the  linear  current  solution  obtained  is  multiplied  with  a  as 
given  from  Equation  A13;  if  Vliu/w  <0.5,  the  linear  solution  is  retained  as 
the  first  guess  without  correction.  This  correction  procedure  saves  one  to 
three  iterations  from  the  typically  four  to  eight  iterations  required  for 
convergence  of  V  to  1  percent  of  the  previous  value  at  any  point  on  the  grid. 
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APPENDIX  B:  SUMMARY  OF  DATA  USED  IN  THE  SIMULATIONS 


Laboratory  Data  from  Visser  <'1982') 

1.  The  data  for  the  current  measurements  are  taken  from  Visser  (1982). 
Because  of  observed  and  potential  disturbances  of  the  current  at  the  basin 
boundaries,  for  the  present  simulations  the  current  was  averaged  from  two 
centrally  located  transects.  The  transects  are  Cases  1,  4,  and  7,  transects  4 
and  5;  Case  3,  transects  8  and  9  (Tables  Bl,  B2 ,  B3,  and  B4) .  Mean  water 
surface  elevations  were  digitized  from  figures  appearing  in  Visser  (1984a) . 


Table  Bl 

Visser  (1982)  Case  1:  Slope  0.101 


Distance  Offshore 

m 

Wave  Height 
m 

Current 

m/sec 

Mean  Water  Surface  Elevation 

m 

0.10 

0.401 

0.30 

0.600 

0.027 

0.50 

0.029 

0.668 

0.016 

0.70 

0.048 

0.666 

0.007 

0.90 

0.057 

0.593 

-0.001 

1.10 

0.075 

0.424 

-0.003 

1.30 

0.098 

0.289 

-0.003 

1.50 

0.103 

0.184 

-0.003 

1.70 

0.096 

0.121 

-0.002 

1.90 

0.083 

0.080 

-0.002 

2.10 

0.074 

0.060 

-0.002 

2.30 

0.075 

0.042 

2.50 

0.077 

0.032 

2.70 

0.079 

0.024 

3.10 

0.082 

0.017 

-- 

Bl 


Table  B2 

Visser  ('1982')  Case  3:  Slope  0.101 


Distance  Offshore 

m 

Wave  Height 
m 

Current 

m/sec 

Mean  Water  Surface  Elevation 

m 

0.18 

-  - 

0.406 

0.020 

0.38 

0.018 

0.456 

0.014 

0.58 

0.031 

0.472 

0.007 

0.78 

0.048 

0.451 

0.002 

0.98 

0.055 

0.418 

-0.003 

1.18 

0.085 

0.354 

-0.002 

1.38 

0.094 

0.269 

-0.002 

1.58 

0.093 

0.201 

-0.002 

1.78 

0.092 

0.139 

-0.001 

1.98 

0.090 

0.095 

-0.001 

2.18 

0.092 

0.062 

2.38 

0.087 

0.039 

2.58 

0.083 

0.024 

2.98 

0.091 

0.009 
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Table  B3 

Visser  (1982)  Case  4:  Slope  0.050 


Distance  Offshore 

m 

Wave  Height 
m 

Current 

m/sec 

Mean  Water  Surface  Elevation 

m 

0.01 

0.000 

0.016 

0.21 

0.114 

0.015 

0.41 

0.275 

0.014 

0.81 

0.016 

0.332 

0.011 

1.21 

0.026 

0.389 

0.008 

1.61 

0.034 

0.375 

0.003 

2.01 

0.056 

0.301 

-0.001 

2.41 

0.087 

0.218 

-0.001 

2.81 

0.084 

0.139 

-0.001 

3.21 

0.078 

0.099 

-0.001 

3.61 

0.080 

0.063 

-0.001 

4.01 

0.079 

0.034 

-0.001 

4.41 

0.076 

0.016 

4.81 

0.076 

0.007 

5.21 

0.075 

0.000 

5.61 

0.072 

-0.006 
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Table  B4 

Visser  (1982')  Case  7:  Slope  0,050 


Offshore  Distance 

m 

Wave  Height 
m 

Current 

m/sec 

0.13 

_ 

0.000 

0.33 

0.060 

0.53 

0.102 

0.93 

0.009 

0.148 

1.33 

0.019 

0.194 

1.73 

0.029 

0.202 

2.13 

0.041 

0.172 

2.53 

0.060 

0.120 

2.93 

0.089 

0.066 

3.33 

-- 

0.040 

3.73 

0.024 

4.13 

0.011 

4.53 

0.006 

4.93 

0.001 

5.33 

-0.004 
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Field  Data  from  Kraus  and  Sasaki  (1979') 


2.  These  data  are  taken  from  Kraus  and  Sasaki  (1979),  who  list  current 
measurements  made  at  Urahama  Beach,  Japan,  at  seven  transects  alongshore.  The 
values  in  Table  B5  are  averages  from  all  transects. 


Table  B5 

Kraus  and  Sasaki  ('1979')  Field  Data 


Distance  Offshore 

m 

Depth 

m 

Current 

m/sec 

5 

0.45 

0.63 

10 

0.50 

0.68 

15 

0.52 

0.66 

20 

0.59 

0.66 

25 

0.69 

0.75 

30 

0.69 

0.77 

35 

0.76 

0.51 

40 

0.99 

0.41 

45 

1.11 

0.29 

50 

1.26 

0.28 

55 

0.18 

Field  Data  from  Thornton  and  Guza  (1986) 

3.  The  data  listed  in  Tables  B6  to  B9  are  presented  in  Thornton  and  Guza 
(1989)  and  summarize  measurements  of  wave  height  and  longshore  current  made 
during  3-6  February  1980  at  Leadbetter  Beach  in  California.  The  data  set  is 
also  discussed  in  Wu,  Thornton,  and  Guza  (1985)  and  Thornton  and  Guza  (1986). 
The  distance  offshore  is  given  with  respect  to  the  mean  shoreline  measured 
during  the  data  collection  interval.  The  wave  height  listed  is  the  root-mean- 
square  value,  whereas  the  current  velocity  represents  a  mean  value.  Wave 
characteristics  for  the  measurement  events  are  given  in  Table  1  in  the  main 


B5 


text.  An  asterisk  in  a  table  implies  that  an  average  value  is  given  of  two 
gages  located  close  together. 


Table  B6 

Thornton  and  Guza  (1989)  Field  Measurements  of  3  February 


Distance  Offshore 

m 

Depth 

m 

Wave  Height 
m 

Current 

m/sec 

6.93 

0.37 

0.40 

0.23 

9.93 

0.50 

0.31 

0.11 

12.99 

0.65 

0.50 

0.35 

15.95 

0.80 

0.48 

0.44 

18.75 

0.94 

0.54 

0.44 

21.88 

1.11 

0.60 

0.50 

24.80 

1.24 

0.59 

0.40 

31.00 

1.48 

0.67 

0.39 

35.12 

1.65 

0.68 

0.35 

39.70 

1.77 

0.71 

0.34 

57.16 

2.28 

0.68 

68.96 

3.02 

0.61 

0.15 

80.70 

3.80 

0.55* 

0.08 

128.60 

6.51 

0.49 

-0.04 
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Table  B7 


Thornton  and  Guza 

(1989)  Field 

Measurements  of  4 

February 

Distance  Offshore 

m 

Depth 

m 

Wave  Height 
m 

Current 

m/sec 

10.39 

0.52 

0.36 

0.36 

13.39 

0.65 

0.35 

0.34 

16.45 

0.78 

0.40 

0.44 

19.41 

0.89 

0.44 

0.45 

22.21 

1.01 

0.50 

0.46 

25.34 

1.16 

0.57 

0.48 

28.26 

1.24 

0.60 

0.43 

31.31 

1.34 

0.64 

0.44 

34.46 

1.42 

0.59 

0.38 

38.58 

1.56 

0.63 

0.32 

43.16 

1.70 

0.66 

0.17 

60.62 

2.37 

0.67 

72.42 

3.03 

0.60 

84.16 

3.79 

0.56* 

0.06 

132.06 

6.50 

0.49 

-0.03 
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Table  B8 

Thornton  and  Guza  ( 1989 )  Field  Measurements  of  5  February 


Distance  Offshore 

m 

Depth 

m 

Wave  Height 
in 

Current 

m/sec 

11.36 

0.65 

0.34 

0.26 

14.42 

0.78 

0.38 

0.35 

17.38 

0.88 

0.43 

0.35 

20.18 

0.97 

0.44 

0.34 

23.31 

1.07 

0.47 

0.34 

26.23 

1.12 

0.48 

0.29 

29.28 

1.18 

0.51 

0.30 

32.43 

1.26 

0.50 

0.25 

36.55 

1.37 

0.59 

0.21 

41.13 

1.52 

0.59 

0.18 

47.60 

1.78 

0.13 

58.59 

2.21 

0.55 

70.39 

2.90 

0.49 

82.13 

3.64 

0.45* 

0.04 

130.03 

6.35 

0.40 

-0.03 
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Table  B9 

Thornton  and  Guza  C 1989 Field  Measurements  of  6  February 


Distance  Offshore 

m 

Depth 

m 

Wave  Height 
m 

Current 

m/sec 

7.53 

0.50 

0.20 

0.15 

10.58 

0.65 

0.23 

0.20 

13.23 

0.70 

0.24 

0.16 

13.58 

0.71 

0.27 

0.15 

16.38 

0.77 

0.27 

0.18 

19.51 

0.85 

0.34 

0.19 

22.43 

0.90 

0.31 

0.13 

25.48 

0.95 

0.32 

0.14 

28.63 

1.03 

0.31 

0.08 

32.75 

1.18 

0.37 

0.04 

37.33 

1.38 

0.34 

0.05 

54.79 

2.22 

0.30 

66.59 

2.79 

0.28 

0.02 

78.33 

3.45 

0.26* 

0.00 

126.23 

6.16 

0.26 

-0.03 

Barred  Profile 

of  Ebersole 

and  Dalrvmnle 

61980') 

4.  Table  BIO  lists  the  data  used  to  generate  a  barred  profile  configura¬ 
tion  taken  from  Ebersole  and  Dalrymple  (1980).  The  values  were  digitized  from 
a  figure  appearing  in  that  reference . 


B9 


Table  BIO 

Barred  Profile  of  Ebersole  and  Dalrvmnle  (1980') 


Distance  Offshore 
_ m _ 


Water  Depth 
_ m _ 


-10 

0.00 

0 

-0.27 

10 

-0.50 

20 

-0.87 

30 

-1.24 

40 

-1.67 

50 

-2.02 

60 

-2.31 

70 

-2.50 

80 

-2.69 

90 

-2.89 

100 

-2.98 

110 

-3.16 

120 

-3.06 

130 

-2.81 

140 

-2.45 

150 

-2.45 

160 

-2.93 

170 

-3.41 

180 

-4.01 

190 

-4.63 

200 

-5.02 

210 

-5.31 

220 

-5.55 

230 

-5.72 

240 

-6.03 

250 

-6.28 

260 

-6.46 

270 

-6.66 

280 

-7.00 

BIO 


APPENDIX  C:  NOTATION 


cf  =  bottom  friction  coefficient 


(cf)lin  - 
(<-'f)non  = 


(4>y  )  lin 

(4  y)  non 

(4y)sq  = 

g  = 
h  = 
u  = 
Ub  = 

= 

u*  = 

V  = 

Vb  = 

w  = 


X  = 


7  = 


z  = 


friction  coefficient  associated  with  the  linearized  friction  law 

friction  coefficient  associated  with  the  nonlinear  (quadratic) 
bottom  friction  law 

total  water  depth  (m) 

formal  expression  for  the  longshore  component  of  the  bottom 
friction  stress  (m2/sec2) 

linearized  form  of  fby  (m2/sec2) 

complete  (nonlinear)  form  of  fby,  evaluated  numerically 
(m2/sec2) 

evaluation  of  fhy  by  the  square -wave  approximation  (m2/sec2) 
acceleration  produced  by  gravity  (m/sec2) 
still-water  depth  (m) 

cross -shore  component  of  the  total  current  velocity  (m/sec) 

cross -shore  component  of  the  horizontal  water  particle  velocity 
at  the  bottom  (m/s ec) 

magnitude  of  the  horizontal  water  particle  velocity  (m/sec) 
shear  velocity  (m/sec) 

longshore  component  of  the  total  current  velocity  (m/sec) 

longshore  component  of  the  horizontal  water  particle  velocity  at 
the  bottom  (m/sec) 

auxiliary  quantity  entering  the  square-wave  time-averaged  fric¬ 
tion  approximation  (m/sec) 

cross-shore  coordinate  directed  positive  seaward  (m) 
longshore  coordinate  (m) 

vertical  coordinate  directed  positive  from  still -water  surface 

(m) 


zQ  -  roughness  length  scale,  or  elevation  where  the  wind  speed  W  is 
zero  (m) 

A  =  auxiliary  quantity  used  to  simplify  notation  in  expressing  the 
longshore  current  governing  equation  and  related  to  contribu¬ 
tions  from  lateral  mixing  (m3/sec) 

AHi  =  auxiliary  quantity  used  to  simplify  notation  in  the  numerical 
solution  of  the  longshore  current  equation  and  related  to 
contributions  from  lateral  mixing  (m/sec) 


Cl 


auxiliary  quantity  used  to  simplify  notation  in  expressing  the 
longshore  current  governing  equation  and  related  to  contribu¬ 
tions  from  bottom  friction  (m/sec) 

auxiliary  quantity  used  to  simplify  notation  in  the  numerical 
solution  of  the  longshore  current  equation  and  related  to 
contributions  from  bottom  friction  (m/sec) 

wave  phase  speed  (m/sec) 

drag  coefficient  for  the  wind-water  interface 

drag  coefficient  evaluated  at  the  10-m  elevation  from  the  water 
or  land  surface 

wave  group  speed  (m/sec) 

auxiliary  quantity  used  to  simplify  notation  in  the  numerical 
solution  of  the  longshore  current  equation  and  related  to 
contributions  from  lateral  mixing  (m/sec) 

wave  phase  speed  in  deep  water  (m/sec) 

auxiliary  quantity  in  the  double -sweep  algorithm 

wave  energy  flux  (Nm/m/sec) 

cross-shore  component  of  the  bottom  friction  shear  stress 
divided  by  the  total  water  depth  (m/sec2) 

longshore  component  of  the  bottom  friction  shear  stress  divided 
by  the  total  water  depth  (m/sec2) 

wave  energy  flux  at  grid  cell  i  (Nm/m/sec) 

wave  energy  flux  corresponding  to  the  stable  wave  height 
(Nm/m/sec) 

wave  energy  flux  corresponding  to  the  stable  wave  height  at  grid 
cell  i  (Nm/m/sec) 

auxiliary  quantity  in  the  double-sweep  algorithm  (m/sec) 
wave  height  (m) 

wave  height  at  grid  boundary  i  (m) 

wave  height  threshold  having  exceedance  probability  1-p  (m) 
root-mean- square  wave  height  (m) 
wave  height  associated  with  the  stable  condition  (m) 
cross-shore  component  of  lateral  mixing  term  (m/sec2) 
longshore  component  of  the  lateral  mixing  term  (m/sec2) 
number  of  grid  cells 

auxiliary  quantity  used  to  simplify  notation  in  the  numerical 
solution  of  the  longshore  current  equation  and  related  to 
contributions  from  the  wave  and  wind  driving  terms  (m2/sec2) 

cross -shore  component  of  the  wave  driving  term  (m/sec2) 
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longshore  component  of  the  wave  driving  term  (m/sec2) 
cross -shore  component  of  the  wind  driving  term  (m/sec2) 
longshore  component  of  the  wind  driving  term  (m/sec2) 
radiation  stress  component  in  the  direction  of  the  waves  (N/m) 
shoreward-directed  radiation  stress  component  alongshore  (N/m) 
radiation  stress  component  alongshore  (N/m) 

cross -shore  component  of  the  mean  (time -averaged)  current  speed 
(m/sec) 

longshore  current  speed  (longshore  component  of  the  mean  (time- 
averaged)  current  speed)  (m/sec) 

longshore  current  speed  at  cell  i  (m/sec) 

longshore  current  speed  obtained  from  the  linear  friction  stress 
law  (m/sec) 

longshore  current  speed  obtained  from  the  nonlinear  (quadratic) 
friction  stress  law  (m/sec) 

wind  speed  (m/sec) 

wind  speed  at  the  10-m  elevation  from  the  water  surface  (m/sec) 

|W|sin^,  longshore  component  of  wind  speed  (m/sec) 

auxiliary  quantity  entering  the  square-wave  time-average  fric¬ 
tion  approximation  (m/sec) 

auxiliary  quantities  entering  the  square -wave  time -average 
friction  approximation  (m/sec) 

longshore  current  speed  calculated  with  nonlinear  model  divided 
by  the  speed  calculated  with  the  linear  model 

ratio  of  wave  height  to  water  depth  for  breaking  and  nonbreaking 
waves . 

breaker  index  (ratio  of  wave  height  and  water  depth  at  incipient 
breaking) 

calculated  breaker  index 
measured  breaker  index 

empirical  coefficient  relating  the  stable  wave  height  to  the 
water  depth 

lateral  mixing  eddy  viscosity  coefficient  (m2/sec) 

components  of  the  lateral  mixing  eddy  viscosity  coefficient 
tensor  in  which  subscripts  (i,j)  each  denote  the  coordinate  axes 
components  (x,y)  (m2/sec) 

sin# 
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t)  =  mean  displacement  of  the  water  surface  from  the  still -water 

level  as  induced  by  waves  and  wind  (called  "setup"  if  positive 
displacement  and  "setdown"  if  negative)  (m) 

r) i  =  mean  displacement  of  the  water  surface  at  grid  cell  i  (m) 

9  =  angle  of  wave  crests  to  the  bottom  contours  (m) 

9j_  =  angle  of  wave  crests  with  respect  to  the  bottom  contours  at  grid 

cell  I 

k  =  empirical  decay  coefficient  in  the  wave  dissipation  term 
k  =  von  Karman's  constant  (approximately  0.41) 

A  =  empirical  coefficient  in  lateral  mixing  eddy  viscosity 
Alin  =  value  of  A  evaluated  with  the  linear  current  model 

Anon  =  value  of  A  evaluated  with  the  nonlinear  current  model 

£  =  surf  similarity  parameter 
Jr  -  3.14159265... 

p  =  density  of  water  (kg/m3) 
pa  =  density  of  air  (kg/m3) 

rb  =  instantaneous  total  bottom  shear  stress  (N/m2) 

rbx  =  cross -shore  component  of  the  instantaneous  total  bottom  shear 
stress  (N/m2) 

rby  =  longshore  component  of  the  instantaneous  total  bottom  shear 
stress  (N/m2) 

rs  =  magnitude  of  the  total  shear  stress  induced  by  wind  at  the  water 
surface  (N/m2) 

rsx  =  cross -shore  component  of  the  shear  stress  induced  by  wind  at  the 
water  surface  (N/m2) 

rsx  =  longshore  component  of  the  shear  stress  induced  by  wind  at  the 
water  surface  (N/m2) 

4>  =  incident  angle  of  the  wind 

V>  =  ratio  of  longshore  current  speed  of  linear  model  to  the  magni¬ 
tude  of  the  wave  orbital  velocity  at  the  bottom 

<  >  =  time -averaging  operation 
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